three dimensional modelling of hydrodynamics and dissolved...
Post on 22-Sep-2019
7 Views
Preview:
TRANSCRIPT
Three dimensional modelling of hydrodynamics and dissolved oxygen
in Lake Erie using ELCOM-CAEDYM
Luis F. Leon, Bob Hecky, Ralph Smith, Stephanie Guildford, Kevin Lamb (U. Waterloo)Jose Romero, Jason Antenucci (U. West. Australia)
CAEDYM (Computational Aquatic Ecosystem Dynamics Model)
• designed to be linked to hydrodynamic models...
– ELCOM as the 3D hydrodynamic model driving CAEDYM,– model includes comprehensive process representation of the C, N, P, Si and DO cycles,– several size classes of inorganic suspended solids, and phytoplankton dynamics,– operates on any sub-daily time step to resolve algal processes (diurnal photosynthesis and nocturnal respiration),– generally runs at the same time interval as the hydrodynamic model.
Developed, supported, and distributed by Centre for Water Research, University of Western Australia
‘Dead Zone’ in Lake Erie...
Source: US EPA http://www.epa.gov/glnpo/glindicators/water/oxygenb.html
Initial Testing: 1994• Including main tributaries and tracers setup...
– river tracers (flow weighted)– area tracers (i.e. basin exchanges)
TP Loads (concentrations)
Next step: Setup initial conditions & output sites…
With assistance and collaboration of Dave Schwab (GLERL) and Dave Dolan (UWGB).
Sampling Sites...
Initial Conditions...Average values for post-invasion years (1989-1993) in SpringParameter Unit West Central East
TP µgP/l 20.3 10.7 10.3SRP µgP/l 2.8 2.8 4.6NO3 mgN/l 0.8 0.28 0.28TKN mgN/l 0.18 0.16 0.14NH3 µgN/l 34.9 6.6 4.3
Silicate mg/l 1.4 0.4 0.7Chla µg/l 4.1 3.8 1.3
• IC (constant values by basin)
• IC (spatial variation)
Source: Makarewicz et al. 2000 JGLR 26:82
Source: Cruise values from NWRI Star Database
Config & Datablock Output...
Sheets:Surface Layer
Vertical Averages
Bottom Layer
Curtains:Longitudinal
West-Central
Central-East
Profiles:Profile (6 sites)
Profile WQ (27 sites)
TS @ 3m depth
TS@ -3m
Schematic of oxygen cycle in CAEDYM
Results (Animations)...
1994_85d_Erie_5sht.movSheets: Average5 Parameters (i.e. Temp, DO, TP, TN & Chla)
1994_85d_Erie_3Params.movSheets: Average & Curtains: Longitudinal
3 Parameters (i.e. Temp, DO & TP)
1994_85d_Erie_TEMPsvb.mov1994_85d_Erie_DOsvb.mov
Sheets: Surface, Average & Bottom
Individual Parameters (i.e. Temp & DO)
1994_85d_Erie_Temp.mov1994_85d_Erie_DO.mov
Sheet: Average & Curtain: Longitudinal
Individual Parameters (i.e. Temp & DO)
MoviesDescriptionOutput
Simulation Details:Starting Date: 12-May-1994 Time Step: 900s 8,160 iterations (85 days) Running Time: 27hrs Ratio: 80 (with CAEDYM, 350 ELCOM only)
Results (Profiles)...
Verification (DO spatial)...
Verification (DO profiles)...MSN 879 (23,151) : Sound Depth 62.5 m
0
10
20
30
40
50
600 5 10 15
DO (mg/l)
Dep
th (m
)
MSN 940 (29,154) : Sound Depth 44 m
0
10
20
30
40
50
600 5 10 15
DO (mg/l)
Dep
th (m
)
MSN 880 (57-77) : Sound Depth 25 m
0
10
20
30
40
50
600 5 10 15
DO (mg/l)
Dep
th (m
)
MSN 965 (79-42) : Sound Depth 13 m
0
10
20
30
40
50
600 5 10 15
DO (mg/l)
Dep
th (m
)
MSN 972 (58_15) : Sound Depth 8 m
0
10
20
30
40
50
600 5 10 15
DO (mg/l)
Dep
th (m
)
MSN 968 (66_32) : Sound Depth 10 m
0
10
20
30
40
50
600 5 10 15
DO (mg/l)
Dep
th (m
)
West Central East
Verification (TP spatial)...
TS - Main Variables...
TS - Dissolved Oxygen...
Details…
Schematic of P cycle in CAEDYM
P components (PO4)...
Back…
P components (POP)...
Back…
P components (DOP)...
Back…
TS - Nitrogen...
Details…
Details…
Details…
Details…
TS – Phytoplankton (Cryptophytes & Diatoms)...Details…
TS – Phytoplankton (Cyanobacteria & Nodularia)...
Next steps...
• Fine tune parameters, initial data and loads...– Revise algal group parameterization (done)– Gather Grand River daily flows (done).– Extend simulation to October (done)
• Additional setup for...– 2002 (loading and cruise data). LETS– 2004 & eventually 2005.– Detailed (higher resolution) nearshore East Basin.
N components (NO3)...
Back…
N components (NH4)...
Back…
N components (DON)...
Back…
N components (PON)...
Back…
Phytoplankton TS-Output...
Back…
}
Objectives...
• Learn how to use the coupled ELCOM-CAEDYM & preliminary setup-testing in Lake Erie...
– Program setup for version 2.2 and running executables.– Input file generation (loads and initial conditions, output definition, parameter configuration file, etc.).– Output transformation to netCD files (web and in-house scripts).– CWR Modeler as a quick output browser.– MatLab as the output visualization tool (plots, contours, time series, animations, etc.).
TS - Phosphorus...
Details…
Details…
Details…
Main Processes...
Programs Schematic...
top related