Shellfish Farm

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,

  1. Go to the Modules tab

  2. Expand Hydrodynamics sub-tab

  3. 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,

  1. Go to the Modules tab

  2. Expand Hydrodynamics sub-tab

  3. 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,

  1. RMC on Salinity sub-tab under Modules tab to open Salinity Parameters form

  2. Click on Assign button to set salinity IC

  3. Under Value to Set frame, set constant value of salinity equal to 30.1 ppt

  4. 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,

  1. browse to the "Data_Salinity.dat" file and select the data format as EE Dat or WQ Data;

  2. Modify Import Settings parameters; then

  3. click OK to import data series; Import Data form appears to allow the user to select import options from this form

  4. 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