NSRA logo US-Canada NSRA Model - Version 4.3

(Manuscript version, 01/19/2021)

Use 'Ctrl + F5' to reset.

Overview and Disclaimer

Welcome to the NSRA tool. This tool accompanies the manuscript titled “Risk assessment of norovirus illness from consumption of raw oysters in the United States and in Canada” published in Risk Analysis. (DOI: 10.1111/risa.13755) (link will open in a new tab).

Manuscript Abstract

Human norovirus (NoV) is the leading cause of foodborne illness in the United States and Canada. Bivalve molluscan shellfish is one commodity commonly identified as being a vector of NoV. Bivalve molluscan shellfish are grown in waters that may be affected by contamination events, tend to bioaccumulate viruses, and are frequently eaten raw. In an effort to better assess the elements that contribute to potential risk of NoV infection and illness from consumption of bivalve molluscan shellfish, the US Department of Health and Human Services/Food and Drug Administration (FDA), Health Canada (HC), the Canadian Food Inspection Agency (CFIA), and Environment and Climate Change Canada (ECCC) collaborated to conduct a quantitative risk assessment for NoV in bivalve molluscan shellfish, notably oysters. This study describes the model and scenarios developed and results obtained to assess the risk of NoV infection and illness from consumption of raw oysters harvested from a quasi-steady-state situation. Among the many factors that influence the risk of NoV illness for raw oyster consumers, the concentrations of NoV in the influent (raw, untreated) and effluent (treated) of wastewater treatment plants (WWTP) were identified to be the most important. Thus, mitigation and control strategies that limit the influence from human waste (WWTP outfalls) in oyster growing areas have a major influence on the risk of illness from consumption of those oysters.

Access to the tool is provided to illustrate the model presented. Please refer to the manuscript for a full description of the methods, results, limitations and conclusions of the analysis.

Responsibility for the interpretation and use of the model and of the accompanying documentation lies solely with the user. Third parties' use of or acknowledgment of the system and its accompanying documentation, including through the suggested citation, does not in any way represent that US Food and Drug Administration, Health Canada, Canadian Food Inspection Agency, Environment and Climate Change Canada or JIFSAN or JIFSAN endorses such third parties or expresses any opinion with respect to their statements.

Suggested citation

Pouillot R., Smith M., Van Doren J.M., Catford A., Holtzman J., Calci K.R., Edwards R., Goblick G., Roberts C., Stobo J., White J., Woods J.,DePaola Jr. A., Buenaventura E.,Burkhardt III W. (2021) Risk assessment of norovirus illness from consumption of raw oysters in the United States and in Canada. Risk Analysis. DOI: 10.1111/risa.13755 (link will open in a new tab).

Instructions

Navigate from tab to tab to set the various parameters (Simulation setup, Influent, Waste Water Treatment Plant characteristics, Harvest Water (Estuary) characteristics, Shellfish characteristics, Meal characteristics, Dose Response). Once your set of parameters is defined:
  • if someone has already launched the simulation with the same set of parameters, “Computed Results Available” will be indicated and the results will be available immediatly. click on the “Compute/Update Results” button to get those results.
  • Otherwize, “To Be Computed” will be indicated. click on the “Compute/Update Results” button and wait until the simulation is done (see progress bars in the lower right corner).
In the “Input/Output” tab, you'll be able to get a report of the simulation, to save the current set of parameters or to upload a prviously saved set of parameters.

Various parameters

Manuscript derived using 10,000 simulations


Note: For meals including more shellfish, the meal will use a sampling with replacement. Specifying a high value will slow the simulation.
Note: a value of 0 (default) computes the number of classes using a classical statistical method (Sturges, 1926).
Note: Density graphs may be misleading near the limits.

Influent

Proportion of infective virus in influent

The following parameters will be used in a Pert distribution. If you want the proportion to be constant, set the minimum, the most probable and the maximum to the same value. Default is 1/40.
Norovirus GI
Norovirus GII

Specification of influent concentration

Norovirus GI
Norovirus GII
MSC

Variation in virus concentration over time. Default values: set to 0.


Results

Influent (log10 per l)

WWTP

Specification of effluent concentration
Norovirus GI
Norovirus GII
MSC
Note: The proportions will be normalized to 1.
Variation in virus concentration over time. Default: set to 0.

Results

Effluent (log10 per l)
Log-Reduction (log10)
Influent-Effluent (log10 per l)
Norovirus I
Norovirus II
MSC

Harvest Water parameters

Harvest Water characteristics.
Mixed semi-diurnal tides are considered Semi-Diurnal tides
You can use the "Box-Model" in the "Tools" tab to estimate the dilutions. After calculating the value in that tool, copy the numbers down and input them here.
(Set the same value in both boxes if you are considering only one estuary)
(Set the same value in both boxes if you are considering only one estuary)
Derive value uses a Pert distribution with min = 1 hour, mode = 1/2 tide period and max = 1 tide period.
Adapt the parameters to your situation

Treated Wastewater

Concentration

The virus concentration in treated wastewater will be the same as the concentration of virus in effluent, drawn from the WWTP model.

Flow

Dilution to harvest

The default dilution to harvest is specified above, in the harvest water parameters section. Change here if you expect a change in the dilution during the exceptional event (e.g. wet weather). You can use the "Box-Model" in the "Tools" tab to estimate the dilutions. After calculating the value in that tool, copy the numbers down and input them here.
(Set the same value in both boxes if you are considering only one estuary)
(Set the same value in both boxes if you are considering only one estuary)

Time to reach harvest area

The time to reach the harvest area is specified above, in the harvest water parameters section.

Concentration

Partially treated wastewater is a mix between influent (untreated wastewater) and effluent (treated wastewater). A value of 0.4 would simulate a mix of 40% untreated wastewater, 60% treated wastewater. A value of 0 would simulate treated wastewater only, and a value of 1 would simulate untreated wastewater only.

Flow

Please input a proportion value, representing the amount of partially treated wastewater flowing during the exceptional event, relative to the design capacity.

Reduced Log Reduction

If the log reduction was altered during the exceptionnal event, provide here the effect. Example: 2 means that the log reduction will be 2 log lower than usual.

Lack of Disinfection

If the disinfection is off during the exceptionnal event: check the box.

Dilution to harvest

Same as treated wastewater.

Time to reach harvest area

Same as treated wastewater.

Concentration at emission

The virus concentration in wastewater will be the same as the concentration of virus in influent, drawn from the WWTP model. However, if this concentration has been diluted, input the factor by which is has been diluted. 1 in 100 corresponds to wastewater diluted 100 times at emission. 1 in

Dilution to harvest

The dilution to harvest for this source is different from sources coming from the WWTP outfall pipe, i.e., treated wastewater and partially treated wastewater. You can use the "Box-Model" in the "Tools" tab to estimate the dilutions
(Set the same value in both boxes if you are considering only one estuary)
(Set the same value in both boxes if you are considering only one estuary)

Time to reach harvest area

In hours.

Concentration

The concentration is expressed in log10 gc/ml (Note: the median peak amount of virus shedding from experimentally infected individuals was 1010 gc/g (Atmar et al, 2008)).

Dilution to harvest

The dilution to harvest for this source is different from sources coming from the WWTP outfall pipe. A value of 0 would be pure sewage on the animal. You can use the "Box-Model" in the "Tools" tab to estimate the dilutions. After calculating the value in that tool, copy the numbers down and input them here.
(Set the same value in both boxes if you are considering only one estuary)
(Set the same value in both boxes if you are considering only one estuary)

Time to reach harvest area

In hours.

Relative contamination at point of harvest. This output value indicates an order of magnitude of the factor by which one can multiply the nominal values, to understand the impact of the exceptional event.

Total NoV

Infective NoV, under the assumption of a 1 log10 inactivation per 24 hours


Relative contribution to contamination
Does not consider tides. Assumes 1 log10 inactivation per 24h.






Harvest Water Results

Harvest Water (log10 per 100 ml)
Concentration in Water at Harvest (Total Infective NoV (GI and GII)/ 100 ml)
(Negative time is time to reach Steady-state)

Shellfish at Harvest (Bioaccumulation an Depuration) parameters

Weight specification
Elimination specification
Bioaccumulation specification
Flannery et al (2012) does not have Single occasion among [oyster] wn occasion
Biancani et al (2012) does not have Cold Season
Only Burkhardt & Calci (2000), Single oyster single occasion, can be used with the Exponential model

Shellfish at Harvest Results


Oyster at Harvest (one oyster per iteration)

Oyster at Harvest (sum per dozen)

Prevalence (%, per animal)

Prevalence (%, per dozen)

Expected concentration in animals (Infective NoV)

(Negative time is time to reach Steady-state)

Processing and Distribution

Depuration, relay and wet storage specifications

Results

Oyster after process (one oyster per iteration)

Oyster after process (sum per dozen)

Prevalence (%, per animal)

Prevalence (%, per dozen)

Meal Characteristics

Oysters per meal
Preparation

Results

Time Temperature Profile
Log10 reduction as a function of time
In blue: virus particles, In black: infective particles

Ingested NoV (per meal)

Hazard Characterisation (Dose Response)


Results

Mean sum of infective GI and GII per meal

Expected number of infections for this number of meals

Expected number of illnesses for this number of meals

Risk of Infection

Risk of Illness

Evolution of the Risk per Day

(Negative time is time to reach Steady-state)

Risk of Illness per Day (x 100)


Save the current report

Download

Upload a set of parameters

Save the current set of parameters

Save

Parameters for this simulation

Calculation tools for Harvest water dilution parameter

Box-Model for Dilution Estimate

Nominal situation
The following tool uses a box-model to estimate dilution. This panel helps to document the fields 'Dilution at mean Tide' and 'Dilution at low tide' in the 'Harvest Water (Parameters)' tab. Document the inputs according to the situation Qo: Daily average Flow leaving WWTP through outfall
Qf: Average freshwater flow - measured or estimated (from tributary into estuary)
T: Type of tide
Ai: Area of impact
Water Depth


Wet Weather situation
Qow: Flow leaving WWTP during wet weather
Qfw: Fresh Flow during wet weather (from tributary into estuary)

Write down the outputs of this table and use them as inputs for the Dilution, or Dilution to harvest parameters. For example, fill out the required information to determine 'Dilution at mean tide', record the output and transfer to the appropriate area. The same can be done for other dilutions required such as 'Dilution at low tide' or 'Dilution to harvest' for an exceptional event.

Dilution Model Estimate

The following tool uses a near field and far field model to estimate dilution. This panel helps to document the fields 'Dilution' in the 'Harvest Water (Parameters)' tab. Document the inputs according to the situation
(ρae)/(ρe) = 0.027 (Sharp, 1994)

Initial Dilution


Final Dilution

Write down the outputs of this table and use them as inputs for the Dilution, or Dilution to harvest parameters.