By: Melissa Mika, P.E.

Introduction

Snowmelt can be a major driver of runoff in watersheds where seasonal snowpack accumulates and ablates during warming periods. In many operational and planning applications, the way snowmelt is modeled can strongly influence simulated hydrographs, volume timing, and peak flow magnitude. Within the Hydrologic Engineering Center's Hydrologic Modeling System (HEC-HMS), snow accumulation and ablation (melt) can be simulated using the Temperature Index snowmelt method, which is widely used due to its moderate data requirements and robust performance when properly parameterized.

A key input to the Temperature Index method is the ATI-Meltrate function, a paired-data relationship that defines melt rate as a function of Antecedent Temperature Index (ATI). The ATI-Meltrate function is often among the most sensitive parameters in Temperature Index snow modeling because it governs how melt responsiveness evolves through a season as the snowpack warms and ripens.

Estimation of the ATI-Meltrate function can be difficult, and typically relies on iterative calibration processes. To improve this workflow, the Snowmelt Parametric Aggregation Calculator (SnowPAC), originally distributed as a Jython-based utility for HEC-DSSVue and Corps Water Management System (CWMS) environments, has now been implemented directly within HEC-HMS Version 4.14. The integrated tool provides an interactive, project-based workflow to estimate ATI-Meltrate functions from observed time series data and exports results to the HMS project, substantially reducing the effort required to move from raw snow observations to a defensible initial parameter set. The new SnowPAC wizard provides numerous enhancements over the original Jython-based tool, including visualization enhancements, a new interpolation method, and the ability to view and compare multiple curves together.

Why ATI-Meltrate Functions Matter

In the Temperature Index method, melt is assumed proportional to the degree to which air temperature exceeds a base (melt) temperature. However, field conditions rarely behave as a single constant melt factor across the entire season. Early-season melt events often produce modest melt per degree-day, while later events (after the pack becomes isothermal and liquid water content increases) can generate substantially higher melt per degree-day. The ATI-Meltrate function captures this seasonal transition by linking melt rate to accumulated thawing conditions represented by ATI.

For basins where snow cover is shallow or short-lived, a nearly constant melt-rate relationship may be sufficient. In contrast, where snow persists for weeks to months and SWE becomes significant, a variable ATI-Meltrate function can meaningfully improve simulation realism and reduce calibration effort by providing a better-informed starting point.

From Observations to ATI-Meltrate

SnowPAC estimates ATI-Meltrate functions using three daily time series at a common location (or analysis point):

• Air temperature
• Precipitation
• Snow water equivalent (SWE)

The tool is designed around a practical reality of snow observations: the most informative periods for estimating melt behavior are discrete windows when temperatures remain above freezing and the snowpack is actively melting. SnowPAC formalizes these windows as melt events and provides tools to find, filter, and refine them.

Within SnowPAC, ATI is computed in a manner consistent with Temperature Index concepts, accumulating thawing degree-days when temperature exceeds the base temperature (commonly 32°F). SnowPAC then computes cumulative melt from decreases in SWE, while also checking conditions intended to isolate “clean” melt behavior – the tool uses precipitation and an adjustable rain rate limit to avoid attributing SWE changes to melt during rain-on-snow or mixed precipitation periods where SWE behavior may reflect processes beyond temperature-driven melt alone.

A defining element of SnowPAC’s event-based approach is that when temperatures fall below the base temperature, the tool resets ATI and cumulative melt to zero, thereby segmenting the record into distinct melt events that each begin at a consistent origin. This standardization makes it easier to compare events across years and to develop a single ATI-Meltrate function that represents typical basin response rather than a one-off seasonal trajectory.

Integrated HMS Workflow

The new implementation is accessed in HEC-HMS through the Snowmelt tools menu (Tools → Snowmelt → SnowPAC). Because it is integrated with the HMS project structure, SnowPAC pulls input data from Time-Series Data gages already defined in the project. This reduces file handling overhead and encourages a consistent process: perform quality control on the daily time series (prior to input to the tool), define the analysis window, identify representative melt periods, and then convert observed behavior into an HMS-ready paired data function.

After the user selects the precipitation, temperature, and SWE gages, SnowPAC determines the overlapping period of record and allows the analysis window to be restricted to a continuous period (for example, a particular water year or multi-year span). The tool then identifies candidate melt events based on user-controlled thresholds (including base temperature and minimum ATI criteria). This capability is particularly useful when working with long records that include numerous small or noisy thaw periods that are not representative of seasonal melt dynamics.

Once melt events are selected, SnowPAC provides interactive refinement to produce ATI-Meltrate functions for each event using one of two approaches:

• Linear interpolation, where the user specifies ATI ordinates and the tool constructs a piecewise relationship; the resulting melt-rate behavior is typically stepwise between ATI bands.
• Cubic spline interpolation, intended to represent more gradual changes in melt sensitivity and to better reflect natural variability where melt rate evolves smoothly with warming conditions.

The event-by-event workflow supports developing multiple candidate curves from different years or events, which can be compared before selecting one for adoption in the basin model. The resulting ATI-Meltrate function can then be exported to the HMS project for immediate use in the Temperature Index method and can also be written to DSS for archiving or external review.

Practical Benefits for Model Development and Forecasting

By moving SnowPAC into HEC-HMS, the workflow for developing ATI-Meltrate functions has become faster, more transparent, and easier to repeat across multiple analysis points. The integrated approach also supports better documentation and reproducibility: input gages, analysis windows, melt event selection criteria, and interpolation choices are all part of a consistent, user-driven process within the modeling environment.

Most importantly, SnowPAC is aimed at improving the starting point for Temperature Index snowmelt modeling. Instead of relying solely on generalized recommended meltrate values and then iteratively adjusting through calibration, users can derive a curve grounded in observed SWE response under melt conditions. This often translates into improved initial simulations and more efficient calibration, especially in operational contexts where reliable snowmelt representation is required for real-time forecasting.

Learn More

The SnowPAC tool and a step-by-step tutorial are documented on the HEC-HMS documentation space:

• Technical Reference Manual: SnowPAC: Snowmelt Parametric Aggregation Calculator
• Tutorial: Estimating ATI-Meltrate using SnowPAC
• HEC-HMS Quarterly Webinar