Calculate sediment porosity

This function takes median grain size at a location and returns a measure of porosity. The default parameterisation comes from Matt Pace's thesis for fine sediments. The numerical order of parameters has been changed from Matt's to be consistent with StrathE2E.

Usage

D50_to_porosity(D50, p1 = -1.035, p2 = -0.314, p3 = -0.435, p4 = 0.302)

Arguments

D50

The median grain size at a location.

p1

A parameter in the relationship.

p2

A parameter in the relationship.

p3

A parameter in the relationship.

p4

A parameter in the relationship.

Value

A numeric vector containing estimates of sediment porosity in

See also

Other Sediment properties: mud_to_permeability()

Examples

data <- data.frame(D50 = seq(0, 1, length.out = 1000)) %>%
mutate(Porosity = D50_to_porosity(D50))

ggplot(data) +
geom_line(aes(x = D50, y = Porosity)) +
scale_x_continuous(trans = "log", breaks = c(0.001,0.01,0.1,1))
#> Warning: Transformation introduced infinite values in continuous x-axis