Title: Estuarine
1Estuarine Coastal Modelling Course
MODULE IX WATER QUALITY MODELLING PROCESS
Michael Hartnett, Environmental Change Institute,
NUI Galway
2- WATER QUALITY MODELLING PROCESS
OBJECTIVE Guidance on the major aspects in
developing hydrodynamic and water quality models
3- WATER QUALITY MODELLING PROCESS
CONTENT- I Boundary Conditions II Initial
conditions III Hydrodynamics IV Pollutant
loads V Water quality VI Heavy
metals VII Pathogens
4I BOUNDARY CONDITIONS
5I BOUNDARY CONDITIONS
- Boundary conditions (BCs)?
-
- Dominant processes along the boundaries
(edges) of a model. -
- BCs propagate throughout the interior of the
model (Domain) - resulting in some physical, chemical or
biological effects.
6I BOUNDARY CONDITIONS
- Example
- Specification of a tidal boundary induces
water level changes - and currents within a model domain
- Driving forces behind a model
- Poor specification is a major source of
modelling error
7I BOUNDARY CONDITIONS
BOUNDARY LOCATIONS
Boundaries located far enough away from region of
interest (ROI) Balance between accuracy and
speed Always ask Are the boundaries located
correctly?
ROI
8I BOUNDARY CONDITIONS
- Hydrodynamic BCs
- Tidal boundary elevations
- flows
- Streamline !
- River flows
- Wind
9I BOUNDARY CONDITIONS
- Hydrodynamic BCs
- Tidal elevation common
- Relatively easily obtained tide tables, surveys
10I BOUNDARY CONDITIONS
- Hydrodynamic BCs
- Tidal flows less common
- Difficult to obtain surveys, output from other
model - May be more accurate
11I BOUNDARY CONDITIONS
- Streamline BCs
-
- No friction/shear with water outside domain
- Sea boundary
- Free condtion
12I BOUNDARY CONDITIONS
- River flows
- Velocity V(m/s)
- Time varying
- Contributes volume and energy
- Note Q(m3/s) V x Area of flow
13II INITIAL CONDITIONS
14II INITIAL CONDITIONS
- Initial conditions?
-
- Models are run from time zero to simulation
time. - Initial conditions are the values specified
for all model - variables at time zero.
- Realistic examples
- Domain water level set to low water
- Water currents set to zero
- Dissolved Oxygen set to 3mg/l
15II INITIAL CONDITIONS
Example Orthophosphate initial conditions
Based on salinity Specified at each
model grid
16III HYDRODYNAMICS
17III HYDRODYNAMICS
- Examples of BCs
- Tides constant repeating tide e.g. spring tide
- repeating spring/neap 14 day cycle
- realistic tides from tidal
constituents - River annual/monthly averages
- flow recorded data
- Wind annual/monthly averages
- recorded data
- spatially constant (usually)
18III HYDRODYNAMICS
- Actual examples of BCs
- Dublin Bay
- Cork Harbour
19III HYDRODYNAMICS
Dublin Bay
MODEL BOUNDARIES _____ Streamline _____ Tidal
elevations _____ Flow
20III HYDRODYNAMICS
River flow
River flow
River flow
Tidal elevation
Streamline
21IV POLLUTANT LOADS
22IV POLLUTANT LOADS
- Each pollutant source separately defined
- Number of pollutants
- Location
- Hydraulic flow rate Q(m3/s)
- Pollutant concentration C(mg/l)
- (C1, C2, C3Cn)
- Note Mass loading M(mg/s) Q x C
23IV POLLUTANT LOADS
- Load specification
- Constant
- Time varying
- Model inputs depends on available data
- Ideally time varying flow and load!
24IV POLLUTANT LOADS
Temporal variations
Recorded flow
Recorded nitrate nitrogen
25IV POLLUTANT LOADS
- Example Cork Harbour
- 20 discharges
- Sewage outfalls
- Industrial outfalls
- Rivers
- Marine source
26IV POLLUTANT LOADS
- Outfall loads
- Often based on IPC
- Sometimes limited measured data
- Inputs rarely completely known
- Ensure scenarios are conservative
27IV POLLUTANT LOADS
28IV POLLUTANT LOADS
- River inputs
- Significant nutrient loads
- Fewer sources than outfall
- Highly variable seasonally
- Flow varies
- Nutrients vary
- Usually limited data
29IV POLLUTANT LOADS
- Location of rivers/diffuse
30IV POLLUTANT LOADS
Flow data not always readily available
31IV POLLUTANT LOADS
Nutrient data not always readily available
32IV POLLUTANT LOADS
33V WATER QUALITY
34V WATER QUALITY
COMPLEX
35V WATER QUALITY
- Components
- Hydrodynamic boundary/initial conditions
- Pollutant boundary/initial conditions
- Phytoplankton cycle
- Reaction rates/constants
?
?
36V WATER QUALITY
- Simplified phytoplankton cycle
37V WATER QUALITY
- Phytoplankton Kinetics
- Species present
- Temperature
- Light
- Nutrients
- Transport processes
38V WATER QUALITY
Numerical Representation
CP phytoplankton population (mg Chl-a L-1)
GPI growth rate constant (d-1)
DPI death plus respiration rate constant
VS4 settling velocity (m d-1) D
total water depth
39V WATER QUALITY
Numerical Representation
GRTS Arrhenius equation GRNU Monod
equation Limiting Nutrient
40Nitrogen Modelling
41V WATER QUALITY
- Reaction rates/constants
- Determine the rate at which a process progresses
- Must specify rate values for various processes
- Example Phytoplankton growth
CP phytoplankton population (mg Chl-a L-1)
GPI growth rate constant (d-1)
DPI death plus respiration rate constant
VS4 settling velocity (m d-1)
42V WATER QUALITY
- Reaction rates/constants determination?
- Sources
- Experiments
- Texts e.g. Surface Water quality by Chapra
- Specialist literature e.g. Rates, Constants
Kinetics, USEPA - Journal papers
- Modelling manuals
43V WATER QUALITY
- Experiments
- Effects of light on growth rate
44V WATER QUALITY
- Reaction rates/constants
- Mostly from literature
- Often used to tune Water Quality model
- Can vary from report to report
- Sensitivity analysis to assess impacts
- Select conservative values !
45Reaction rates constants
46V WATER QUALITY
References
(1) Brown, L.C. and Barnwell, T.O., (1987), The
Enhanced Stream Water Quality Models QUAL2E and
QUAL2E-UNCAS Documentation and User Manual, U.S.
Environmental Protection Agency. (2) Bowie,
G.L., Mills, W.B., Porcella, D.B., Campbell,
C.L., Pagenkopf, J.R., Rupp, G.L., Johnson, K.M.,
Chan, P.W.H. and Gherini, S.A., (1985), Rates,
Constants, and Kinetic Formulations in Surface
Water Quality Modelling, Environmental Research
Laboratory, Office of Research and Development,
U.S. Environmental Protection Agency. (3)
Chapra, S.C., (1997), Surface Water Quality
Modelling, McGraw-Hill International Editions,
Singapore. (4) Falconer, R.A. and Liu, S.Q.,
(1988), Modelling Solute Transport Using the
QUICK Scheme, Journal of Environmental
Engineering, 114 (1) 3-20.
47V WATER QUALITY
48VI HEAVY METALS
49VI HEAVY METALS
- Toxic metals
- lead, mercury, zinc etc
- Linked to sediments and salinity
50VI HEAVY METALS
- Sediment dependency
- Heavy metal mass fractioned in water
- Some behave as solute
- Some attach to fine, cohesive sediments
- Thus must first model sediment transport
- Very difficult !
- Sediment transport module
- Validated/calibrated etc.
51VI HEAVY METALS
- Salinity dependency
- Heavy metal mass fractioned in water
- P (Particulate bound mass)/(dissolved mass)
- Partitioning coefficient (P)
- Function of water chemistry salinity
- Model salinity
- Determine P
52VI HEAVY METALS
- Salinity dependency
- Very poorly understood research stage
- P varies from estuary to estuary
- Models
- P constant
- P exponential decrease with salinity
- Locally derived function from measurements
53VI HEAVY METALS
Salinity dependency - exponential
54VI HEAVY METALS
Salinity dependency locally derived
55VI HEAVY METALS
Model predictions
56VI HEAVY METALS
- Stages
- Hydrodynamics
- Salinity
- Sediments
- Metal inputs
- Metal transport
- Very difficult !
- Considerable measurements required
57VII PATHOGENS
58VII PATHOGENS
- Stages
- Hydrodynamics
- Pathogen inputs
- Pathogen (solute) transport
- Die-off
- Modelling difficulties
- Loading
- Die-off rate
59VII PATHOGENS
- Die-off rate
- T90, time for 90 pathogens to die
- Typically around 15hours
- T90 range X X hours
- T90 is not a constant
- Temperature
- Sunlight
- Turbidity
- Sediments
?
60- WATER QUALITY MODELLING PROCESS
CONTENT- I Boundary conditions II Initial
conditions III Hydrodynamics IV Discharges/pollu
tant loads V Water quality VI Heavy
metals VII Pathogens
61Estuarine Coastal Modelling Course
MODULE IX WATER QUALITY MODELLING PROCESS
Michael Hartnett, Environmental Change Institute,
NUI Galway