Shellfish Farm
Introduction
The main objective of this training course is to provide users with a basic knowledge of how-to setup a shellfish model and run a solution for Environmental Fluid Dynamics Code Plus (EFDC+). It will cover preparation of the necessary input files for the EFDC+ model and visualization of the output by using the EFDC+ Explorer Modeling System (EEMS). Kobayashi et al. (1997) is used as a case study to illustrate the general process of developing a model from basic to complex. All files for this tutorial are contained in the EEMS website. After this training, users will be able to:
1. Set up a basic water quality and shellfish model from input data using EFDC+ Explorer (EE) version 10.2
2. Visualize, extract and manage model output by EE
Model Description
The model is set up to replicate the Kobayashi et al. (1997) data for shellfish growth in interaction with water quality described in the paper "A population dynamics model for the Japanese oyster, Crassostrea gigas". The goal of the project is to better understand the impact of vegetation to the flow in open channel.
The experiment was carried out in a 300x100 m channel as described in Shellfish Farm#Figure 1 below.
Figure 1. The shellfish farm model grid.
Data Summary
All data used in the model including inputs and comparison data are collected from the paper Kobayashi et al. (1997). These data are stored in "Data" folder from the testcase model (TC-12) available on EEMS website.
Environmental Conditions
Time series of mean water temperature, salinity and chlorophyll are shown in Shellfish Farm#Table 1 and Shellfish Farm#Figure 2 below. Water temperature ranged between 8 to 30oC, while salinity during this period was relatively constant, ranging only between 28.4 to 31.9 ppt.
Table 1. Data use for shellfish model
Day | Temperature(deg. C) | Salinity(ppt) | Chl-a(mg/L) | Dry weight(g C) | Shell length(mm) |
134 | 18.3 | 30.5 | 2.3 |
| 27.14 |
162 | 22.7 | 30.5 | 5.0 | 0.34 | 39.95 |
197 | 26.8 | 28.6 | 2.9 | 0.77 | 63.41 |
226 | 30.3 | 30.8 | 6.3 | 0.61 | 73.86 |
239 | 27.8 | 31.1 | 1.0 | 0.52 | 76.27 |
254 | 28.7 | 31.4 | 1.3 | 0.69 | 76.75 |
270 | 25.5 | 28.8 | 0.3 | 0.68 | 80.23 |
288 | 21.5 | 29.2 | 2.2 | 0.82 | 82.38 |
336 | 15.5 | 29.6 | 1.1 | 1.53 | 86.53 |
393 | 8.0 | 32.0 | 7.3 | 2.81 | 93.08 |
Figure 2. Environments data for the study.
Biological Characteristics of Oyster Populations
Shellfish species (C02) | Oyster |
Initial carbon per unit of shellfish | 0.0262 (g C) |
Fraction of carbon in dry meat weight | 0.4 |
N/C ratio of shellfish | 0.167 (g N/g C) |
Caloric content of shellfish | 5210 (cal/g dry weight) |
Factor for length-weight relationship | 77.9 |
Exponent for length-weight relationship | 0.291 |
Shellfish filtration (C03) | |
Option for temperature effect on filtration | Kobayashi et al. (1997) |
Option for salinity effect on filtration | No salinity effect |
Option for sediment effect on filtration | No TSS effect |
Option for dissolved oxygen effect on filtration | No DO effect |
Shellfish respiration (C04) | |
Option for respiration model | Korean authors (2019) |
Base temperature | 20 (deg C) |
Basal metabolic rate | 0.005 (1/day) |
Temperature coefficient for respiration | 1.066 |
Energy consumption by respiration | 4.83 (cal/mL O2) |
Shellfish growth (C05) | |
Growth efficiency | 0.7 |
Assimilation efficiency | 0.75 |
Assimilation of cyanobacteria | 0.85 |
Assimilation of diatom algae | 0.85 |
Assimilation of green algae | 0 |
Assimilation of POC | 0.5 |
Shellfish mortality (C06) | |
Temperature dependent mortality flag |
|
Temperature coefficient for shellfish mortality | 1.066 |
Specific mortality rate | 0.0058 (1/day) |
Shellfish spawning (C07) | |
Option for spawning | Spawning is accounted for |
Min. dry weight for spawning | 0.3 (g C) |
Min. temperature for spawning | 25 (deg C) |
Fraction of cumulative reproductive biomass for spawning | 0.05 |
Egg dry weight | 1.3E-05 (g C) |
Caloric content of egg | 6133 (cal/g dry weight) |
Shellfish interaction with water quality (C10, C11) | |
Fraction of RPOC from dead shellfish | 0.1 |
Fraction of LPOC from dead shellfish | 0.1 |
Fraction of DOC from dead shellfish | 0.8 |
Fraction of RPOC in feces | 0.1 |
Fraction of LPOC in feces | 0.1 |
Fraction of DOC in feces | 0.7 |
Fraction of SUU in feces | 0.1 |
Fraction of ROP in urine | 0.1 |
Fraction of LOP in urine | 0.1 |
Fraction of DOP in urine | 0.1 |
Fraction of P4D in urine | 0.7 |
Fraction of RON in urine | 0.1 |
Fraction of LON in urine | 0.1 |
Fraction of DON in urine | 0.1 |
Fraction of NHX in urine | 0.7 |
Step-by-Step Guide for Building Hydrodynamics Model
This chapter will show users how to:
Generate a new hydrodynamic model
Assign initial conditions
Assign BCs.
Generate a New Model
In this exercise, we will generate a new water quality model. This exercise familiarizes users with the basic interface of EE and demonstrates how to build a model from scratch. Double click on the EE10.2 icon on your Desktop or start it from the list of programs to load the model interface. Click on the New Model button to begin generating a new model, as shown in Shellfish Farm#Figure 3.
Figure 3. Cartesian Grid Generator form.
For this demo, we will generate a Cartesian grid using the following steps. From the Cartesian Grid Generator form,
(1) Select Grid Options as Generate Uniform Grid;
(2) Set the dimension of the uniform grid as Figure 3
(3) In the example, leave the UTM Zone as 0 or any zone;
(4) After filling in the parameters, click the Generate button to generate the model grid, then click on the Remove Dry button to remove cells outside the polygon. The preview of the generated grid will show in the right panel.
Note: If the grid is not suitable for your purpose, you can modify the grid parameters and click Generate again to re-generate the grid;
(5) After the grid is generated, click the OK button to close the form and go to model settings.
After the grid is generated, the Model Control form will be displayed as Figure 4. General information about the model will be displayed that includes the model title, area, number of cells, model width, height, and model location coordinates.
Figure 4. EE Model control form.
In the model interface, click on Save to save the model.
After clicking on the Save button, a Select Directory pop-up will appear to select saving options. Before project files are saved, the grey box is empty. This is where EE project files are listed. The empty gray box shows that no files have been written yet for the study.
There are three save options:
Full Write can be selected to save all the input files,
Write All except Time Series Files can be selected to save all input files without the time series.
Save Profile File Only (EFDC.EE) can be selected for saving if users have only changed formatting options in EE.
In this project, it is recommended to use Full Write option. After choosing this option, click OK to close the form and save the model.
Figure 5. Saving the model.
Assign Model Bathymetry
The next step in the model setup is to assign the bathymetry for the model domain. Assigning bathymetry refers to assigning the elevation value to each cell in the model. In this exercise, bottom elevation is a constant value at elevation of -3m (see Shellfish Farm#Figure 6).
Figure 6. Set constant bottom elevation.
Assign Water Surface Elevation
Water surface elevation can be assigned following the import of bathymetry. In this exercise, water depth of the channel is 3m, so assign WSEL equal to 0m.
Assign Boundary Condition Time series Data
From 2DH View, assign Inflow and Outflow BCs as shown in Shellfish Farm#Figure 7 below and set constant value of inflow and outflow BCs are 10 m3/s and (-10 m3/s).
Figure 7. Assign Boundary Conditions in 2DH View.
Hydrodynamics Model
From Model Control form,
Go to the Modules tab
Expand Hydrodynamics sub-tab
RMC on Roughness to assign Bottom roughness (Z0) as shown in Figure 8. In this demo, bottom roughness of the model is 0.01
Figure 8. Hydrodynamics – Assign bottom roughness.
From the Model Control form,
Go to the Modules tab
Expand Hydrodynamics sub-tab
RMC on Shallow Water to enter the Dry Depth and the Wet Depth.
This feature is required with fluctuating water levels that will leave some of the domain dry. For an area with tidal flats, the user should probably enable the wet and dry option. Although it incurs a certain computational overhead, it is necessary to obtain accurate results as it sets the minimum water depth allowed in a point before it is taken out of calculation (Dry Depth), and also the water depth at which the point will be reentered into the calculation (Wet Depth). In this exercise shallow water parameters are set as shown in Figure 9.
Figure 9. Hydrodynamics – Shallow Water setting.
Note: The Wet Depth should always greater than the Dry Depth.
Step-by-step Guide for Building Shellfish Farm Model
From Model Control form, RMC on Modules tab under Model Control form, check on the Salinity, Temperature, Water Quality and Shellfish Farm check-boxes then click the OK button to activate the salinity, temperature, water quality and shellfish farm modules as shown in Figure 10 below.
Figure 10. Active EFDC Modules.
Salinity Module
Assign Salinity Initial Condition
From Model Control form,
RMC on Salinity sub-tab under Modules tab to open Salinity Parameters form
Click on Assign button to set salinity IC
Under Value to Set frame, set constant value of salinity equal to 30.1 ppt
Click on Apply button to assign salinity IC as shown in Figure 11 below.
Figure 11. Assign Salinity IC.
Import Salinity Data Series
After assigning salinity IC, from Salinity Parameters form, click on Edit button under Time Series Data frame to import salinity data series as shown in Figure 12 below.
Figure 12. Edit Salinity Time Series Data.
After clicking on the Edit button, Boundary Data Series form will appear for participants to add/modify data series. In Boundary Data Series form, click on the Add New button to add new data series, then select Import Data button to import Salinity series data as below as shown in Shellfish Farm#Figure 13.
In the ASCII Data Import form,
browse to the "Data_Salinity.dat" file and select the data format as EE Dat or WQ Data;
Modify Import Settings parameters; then
click OK to import data series; Import Data form appears to allow the user to select import options from this form
click OK to close the form and finish importing data .
Figure 13. Import data into data series for salinity.
Note: After importing the data series; from Boundary Data Series form, participants can set the data series name by typing in Series Name blank on the top left corner.
Temperature Module
From the Model Control form, RMC on each module sub-tab under the Modules tab; a list of options will appear for the user to select, as shown in Figure 14 below. From this pop-up menu:
(1) Go to the Settings;
(2) Assign the Initial Conditions settings