Water Quality Constituents

Water quality constituents are the HEC-WQ Engine model state variables (e.g., water temperature or dissolved oxygen concentration). The WQ Engine links to three ERDC-EL libraries:

At each time step, the WQ Engine passes information to the libraries for each computational cell (reaches) or layer (reservoirs). This includes:

  • Constituent concentrations,
  • Water surface area,
  • Average depth,
  • Cell or layer volume,
  • Streamwise (reaches) or vertical (reservoir) velocities,
  • Suspended sediment concentration,
  • Water and sediment temperatures (when modeling constituents in the NSMI library),
  • Salinity, and
  • Meteorological data (as detailed in Meteorological Data)

The ERDL-EL libraries use this information to calculate the rate of change in the constituent concentration (\frac{dC}{dt}), which is passed back to the WQ Engine and used to update the cell concentrations. 

A listing of the water quality constituents available in each library is given below. Two constituents - Carbonaceous Biological Oxygen Demand and the General Constituent - support the creation of multiple instances, with each instance allowing unique calibration parameters. The modeled units for each nutrient species are typically given as mass concentrations of the parent nutrient (e.g., ammonium units are mg-N/L). 

In general, HEC-ResSim supports the ability to mix and match water quality constituents. The user selects a combination of constituents to model, and assumptions are made about dependent constituents in order to perform the simulation. For example, water temperature may be the only constituent used in a simulation. In this case, sediment temperature is assumed to be approximately equal the adjoining cell water temperature, and temperature fluxes between the water and the sediment are zero. If both water and sediment temperature are selected, sediment temperature is explicitly modeled, and sediment-water heat fluxes are calculated as described in the WTSM. Similarly, phytoplankton may be modeled without explicitly modeling all dependent nutrient groups. If phytoplankton is modeled without phosphorus groups, for example, phosphorus is assumed to be in abundance and not limit the growth of the phytoplankton. 

There are a few exceptions to this flexibility, however. Since water temperature plays a critical role in determining biochemical reaction rates and the vertical distribution of material in reservoirs, the water quality constituents in the NSMI module may not be modeled without also explicitly modeling water temperature. Sediment temperature modeling similarly requires water temperature. 

Library

Display Name

Units

WTSM

Water Temperature

°C

Sediment Temperature

°C

NSMI

Dissolved Oxygen

mg/L

Carbonaceous Biological Oxygen Demand

mg/L

Ammonium

mg/L

Nitrate

mg/L

Organic Nitrogen

mg/L

Inorganic Phosphorous

mg/L

Organic Phosphorous

mg/L

Phytoplankton

µg/L

Benthic Algae

g-dry weight/m²

Particulate Organic Carbon

mg/L

Dissolved Organic Carbon

mg/L

Dissolved Inorganic Carbon

mg/L

Suspended Particulate Organic Matter

mg-dry weight/L

Bed Sediment Particulate Organic Matter

mg-dry weight/L

Alkalinity

mg-CaCO₃/L

Pathogen

cfu/100 mL

Nitrogen Gas

mg/L

GCSM

Salinity

ppt

General Constituent

mg/L

Water Quality Parameters

Water quality parameters are the coefficients and options that are used to calculate water quality transformations in the ERDC-EL libraries. 

In HEC-ResSim any given parameter may be defined as a constant or time variable. Time variable specifications can be input with a user-defined seasonal cycle, a given DSS time series record, or dependency on a model state variable. For example, the concentration of an inflowing nutrient may be specified as a function of the inflow magnitude.

Parameters may also vary spatially throughout the water quality modeling domain; this is accomplished by defining regions within the water quality geometry. The Upper Sacramento watershed shown below is broken into many regions, colored based on their region id. In this way, individual reservoirs may have their own set of water quality parameters defined, as well as individual stretches of river on the Trinity River, Sacramento River, and Clear Creek. 


A listing of available water quality parameters is given below, organized by their associated constituents. Documentation about each these parameters is provided in Water Quality Transformation Libraries

Some parameters have units specified as "selection". In this case, the parameter is not a numerical value, but a selection from a list of options. For example, the "Hydraulic O² Reaeration Option" specifies the formulation that is used to calculate stream reaeration rates based on hydraulic parameters.

Library

Group

ParameterUnits
WTSMTemperature

Coefficient a in Wind Function

10⁻⁹ 1/mb m/s

Coefficient b in Wind Function

10⁻⁹ 1/mb

Coefficient c in Wind Function

unitless

Turbulent Diffusivity Ratio

unitless

Sediment Layer Thickness

m

Sediment Bulk Density

kg/m³

Sediment Specific Heat Capacity

J/kg/m³

Sediment Thermal Diffusivity

m²/day

Shortwave Radiation Bed Reflectivity

unitless

Background Light Attenuation

1/m

Suspended Sediment Light Attenuation

L/mg 1/m

Suspended Sediment Concentration

mg/L

Wind Sheltering Coefficient

unitless

Shortwave Shading Coefficient

unitless

Longwave Insulation Coefficient

unitless

Constant Sediment Temperature

°C

Turbulent Heat Fluxes Option

selection

NSMIAlgae

Algal Biomass (Dry Weight) Stoichiometry

mg-Dry weight

Algal Carbon Stoichiometry

mg-C

Algal Nitrogen Stoichiometry

mg-N

Algal Phosphorus Stoichiometry

mg-P

Algal Chlorophyll-a Stoichiometry

µg-Chl-a

Maximum Algal Growth Rate

1/day

Algal Respiration Rate

1/day

Algal Mortality Rate

1/day

Algal Settling Velocity

m/day

Light Limiting Constant for Algal Growth

W/m²

Half-saturation N Limiting Constant for Algal Growth

mg-N/L²

Half-saturation P Limiting Constant for Algal Growth

mg-P/L²

NH4 Preference Factor for Algal Growth

unitless

Fraction of Algal Mortality into POC

unitless

Algal Growth Equation Option

selection

Algal Light Limitation Function Option

selection

NSMIBenthic Algae

Benthic Algae Biomass (Dry Weight) Stoichiometry

mg-Dry weight

Benthic Algae Carbon Stoichiometry

mg-C

Benthic Algae Nitrogen Stoichiometry

mg-N

Benthic Algae Phosphorus Stoichiometry

mg-P

Benthic Algae Chlorophyll-a Stoichiometry

µg-Chl-a

Maximum Benthic Algal Growth Rate

1/day

Benthic Algae Respiration Rate

1/day

Benthic Algae Mortality Rate

1/day

Light Limiting Constant for Benthic Algae Growth

W/m²

Half-saturation N Limiting Constant for Benthic Algae Growth

mg-N/L

Half-saturation P Limiting Constant for Benthic Algae Growth

mg-P/L

Half-saturation Density Constant for Benthic Algae Growth

g-Dry mass/m²

NH4 Preference Factor for Benthic Algae Growth

unitless

Fraction of Benthic Algae Mortality into POC

unitless

Fraction of Bottom Area for Benthic Algae

unitless

Fraction of Benthic Algae Mortality into Water

unitless

Benthic Algae Growth Equation Option

selection

Benthic Algae Light Limitation Function Option

selection

NSMICarbon Cycle

POC Hydrolysis Rate

1/day

DOC Oxidation Rate

1/day

Half-saturation Oxygen Attenuation Constant for DOC Oxidation

mg-O²/L

Fraction of CO2 in Total Inorganic Carbon

unitless

Partial Pressure of CO2

ppm

NSMICarbonaceous Biological Oxygen Demand

CBOD Oxidation Rate

1/day

Half-Saturation Oxygen Attenuation Constant for CBOD Oxidation

mg-O²/L

CBOD Sedimentation Rate

m/day

NSMINitrogen Cycle

Organic N Hydrolysis Rate

1/day

Nitrification Rate

1/day

Sediment Release Rate of NH4

g-N/m²/day

Denitrification Rate

1/day

Half-saturation Oxygen Attenuation Constant for Denitrification

mg-O²/L

Sediment Denitrification Velocity

m/day

NSMINSMI Global Parameters

Particulate Organic Matter Light Attenuation

L/mg/m

Linear Algal Light Attenuation

1/m 1/(µg-Chla/L)

Nonlinear Algal Light Attenuation

1/m 1/(µg-Chla/L)⅔

Partitioning Coefficient of Inorganic P

L/kg

Solids Settling Velocity

m/day

Organic N Settling Velocity

m/day

Organic P Settling Velocity

m/day

POC Settling Velocity

m/day

Sediment Oxygen Demand

g-O²/m²/day

Half saturation oxygen Attenuation Constant for SOD

mg-O²/L

Hydraulic O² Reaeration Rate

1/day

Wind O² Reaeration Velocity

m/day

Hydraulic O² Reaeration Option

selection

Wind O² Reaeration Option

selection

pH Numerical Solution Option

selection

NSMIParticulate Organic Matter

POM Settling Velocity

m/day

POM Dissolution Rate

1/day

NSMIPathogen

Pathogen Death Rate

1/day

Light Efficiency Factor for Pathogen Decay

unitless

Pathogen Settling Velocity

m/day

NSMIPhosphorus Cycle

Organic P Hydrolysis Rate

1/day

Sediment Release Rate of DIP

g-P/m²/day

GCSMGeneral Constituent

0-order Reaction Decay Rate

mg/L/day

1st-order Reaction Decay Rate

1/day

Sediment Release Rate

g/m²/day

Settling Velocity

m/day