dissolve_pb: Simulate contributions of various lead solids to total soluble lead

Description

This function takes a water data frame defined by define_water and outputs a dataframe of the controlling lead solid and total lead solubility. Lead solid solubility is calculated based on controlling solid. Total dissolved lead species (tot_dissolved_pb, M) are calculated based on lead complex calculations. Some lead solids have two k-constant options. The function will default to the EPA's default constants. The user may change the constants to hydroxypyromorphite = "Zhu" or pyromorphite = "Xie" or laurionite = "Lothenbach"

Usage

dissolve_pb(
  water,
  hydroxypyromorphite = "Schock",
  pyromorphite = "Topolska",
  laurionite = "Nasanen"
)

Value

A data frame containing only the controlling lead solid and modeled dissolved lead concentration.

Arguments

water: Source water object of class "water" created by define_water. Water must include alk and is. If po4, cl, and so4 are known, those should also be included.
hydroxypyromorphite: defaults to "Schock", the constant, K, developed by Schock et al (1996). Can also use "Zhu".
pyromorphite: defaults to "Topolska", the constant, K, developed by Topolska et al (2016). Can also use "Xie".
laurionite: defaults to "Nasanen", the constant, K, developed by Nasanen & Lindell (1976). Can also use "Lothenbach".

Details

The solid with lowest solubility will form the lead scale (controlling lead solid).

Make sure that total dissolved solids, conductivity, or ca, na, cl, so4 are used in `define_water` so that an ionic strength is calculated.

Examples

Run this code


example_pb <- define_water(
  ph = 7.5, temp = 25, alk = 93, cl = 240,
  tot_po4 = 0, so4 = 150, tds = 200
) %>%
  dissolve_pb()
example_pb <- define_water(
  ph = 7.5, temp = 25, alk = 93, cl = 240,
  tot_po4 = 0, so4 = 150, tds = 200
) %>%
  dissolve_pb(pyromorphite = "Xie")