Field Data Collection

This modeling effort makes use of the network of automated weather stations (AWS) maintained by MCM.  Additionally, it uses the glacier surface mass balance stake network for calibration and validation.  Finally, it uses the streamgauge network for validation.

Model Approach

The point observations were distributed across the valley on a 250 m resolution horizontal grid using the model MicroMet.  Additionally, MicroMet was used to generate forcing for the ~20 m high vertical cliffs that form the termini of the glaciers.

The distributed meteorological data generated by MicroMet was used to force the IceMelt model which calculates surface and below-surface melt and runoff rates for each location in the model domain.  To adequately model glacier ablation and subsurface ice temperatures required including the transmission of solar radiation into the ice and subsequent drainage of some of the melt generated below the surface.  IceMelt was run using a daily time step for the 1995-2005 time period, and then with an hourly time step for the 1996-2011 time period.  The model was calibrated using seasonal surface mass balance measurements, below-surface ice temperatures, and ice density at three locations.  The model was then run in distributed mode across all glacier ablation zones in Taylor Valley.  Model results were compared to surface mass balance stake observations and proglacial streamflow records.

Modeling code used

MicroMet: 

MicroMet is a quasi-physically based meteorological model of intermediate complexity that is suitable for generating gridded atmospheric forcings at high-resolution (~30-1000 m grid increments).  It combines spatial interpolation with known physical relationships between weather and topography for the variables of air temperature, relative humidity, wind speed and direction, incoming solar radiation, incoming longwave radiation, surface pressure, and precipitation. MicroMet includes a preprocessor for filling short gaps in the measurement time series at each station, and makes use of an autoregressive integrated moving average (ARIMA) calculation for filling gaps in diurnally varying variables. MicroMet requires a digital elevation model and the time-series of individual weather variables at various sites in or near the modeling domain. For spatial interpolation from irregularly spaced stations to a gridded dataset, MicroMet uses a Barnes objective analysis scheme.

IceMelt: 

IceMelt is a physical-based surface energy balance model for glacier ice.  It accounts for both the surface energy balance and subsurface solar heating from absorbed solar radiation, allowing it to reproduce the so-called ‘solid state greenhouse effect’ observed near the surface of Dry Valley glaciers.  The surface energy balance includes terms for net shortwave radiation, net longwave radiation, sensible and latent heat fluxes between the ice and atmosphere, heat conduction into or out of the ice, and melt at the ice surface.  Additionally, the model domain extends 15 m into the ice and the model includes heat conduction through that column, including a source term for solar radiation absorbed internally within the ice based on a two-stream approximation.  The model supports partial melting of ice within a grid cell below the surface, and it calculates drainage of water when partial melt surpasses a specified threshold water fraction.  IceMelt is a point model, but can be run in distributed mode where each location in a spatial grid is run independently.  IceMelt has been run with both daily and hourly time steps.

Model input files

Inputs to MicroMet are air temperature, relative humidity, wind speed, incoming solar radiation, and incoming longwave radiation at 3 m height at point AWS locations.  Additionally, a digital elevation model (DEM) of the model domain is required.

Inputs to IceMelt are air temperature, relative humidity, wind speed, incoming solar radiation, and incoming longwave radiation at 3 m height and surface albedo at each location to be modeled. 

Simulation output

MicroMet outputs air temperature, relative humidity, wind speed, incoming solar radiation, and incoming longwave radiation at all locations in the specified model domain.

IceMelt outputs all of the terms of the surface energy balance, and sublimation and melt at the surface.  It also outputs ice temperature, melt/freeze rate, and water drainage at all depths modeled within the ice.

Both models output all fields at every glacier ablation zone location on a 250 m grid of Taylor Valley