Oyster_incr_cross_section {shelltrace}

Formula that interpolates shell increments

Description

Formula that takes the coordinates of digitized shell increments and interpolates a number of sub-increments between them to increase the time resolution of the growth model. The number of interpolated shell increments as well as the relative thickness of these increments is determined by a sinusoidal seasonal model that simulates seasonal variations in shell growth rate. de Winter, N. J. (2017) <doi:10.5194/gmd-2017-137>

Usage

Oyster_incr_cross_section(incr_matrix, cross_section, season_length,
  Tstep=1, Xstep=0.1)

Arguments

incr_matrix

Matrix listing characteristics of each digitized increment

cross_section

Digitized cross section of the shell with shell top, bottom and growth increments relative to a common X-axis

season_length

Length (in days) of the growth season of the studied bivalve

Tstep

The step size (dt) in time (days) used to interpolate the existing shell increments

Xstep

The step size (dx) in X-direction used to interpolate coordinates of shell increments

Details

Sub-increments are reconstructed by interpolating Y-values between the digitized shell increments

Value

List of two data frames:

IncG

Matrix of X- and Y-coordinates of all interpolated sub-increments

subincr_matrix

Revised version of the "incr_matrix" data frame that contains characteristics (sub-increment number, X-value of start of increment) recalculated for every sub-increment

Note

Please cite Geoscientific Model Development paper dealing with the ShellTrace model

Author(s)

Niels J. de Winter

Source

GitHub
Manuscript
Supplementary data
Author website

References

de Winter, N. J.: ShellTrace v1.0 ? A new approach for modelling growth and trace element uptake in marine bivalve shells: Model verification on pacific oyster shells (Crassostrea gigas), Geosci. Model Dev. Discuss., https://doi.org/10.5194/gmd-2017-137, in review, 2017.

Examples

  Lsub<-Oyster_incr_cross_section(incr_matrix3, cross_section, season_length=250,
  Tstep=1, Xstep=1)