Download PDF
Download page Water Quality Constituents and Parameters.
Water Quality Constituents and Parameters
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:
- Water Temperature Simulation Module (WTSM)
- Nutrient Simulation Module I (NSMI)
- General Constituent Simulation Module (GCSM)
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 | Parameter | Units |
---|---|---|---|
WTSM | Temperature | 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 | ||
NSMI | Algae | 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 | ||
NSMI | Benthic 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 | ||
NSMI | Carbon 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 | ||
NSMI | Carbonaceous Biological Oxygen Demand | CBOD Oxidation Rate | 1/day |
Half-Saturation Oxygen Attenuation Constant for CBOD Oxidation | mg-O²/L | ||
CBOD Sedimentation Rate | m/day | ||
NSMI | Nitrogen 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 | ||
NSMI | NSMI 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 | ||
NSMI | Particulate Organic Matter | POM Settling Velocity | m/day |
POM Dissolution Rate | 1/day | ||
NSMI | Pathogen | Pathogen Death Rate | 1/day |
Light Efficiency Factor for Pathogen Decay | unitless | ||
Pathogen Settling Velocity | m/day | ||
NSMI | Phosphorus Cycle | Organic P Hydrolysis Rate | 1/day |
Sediment Release Rate of DIP | g-P/m²/day | ||
GCSM | General 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 |