olm.calcite.dissRateFromCaPCO2

olm.calcite.dissRateFromCaPCO2(Ca, PCO2, T_C, rho=2.6, method=None, impure=True, per_tol=0.001, error=False, error_num=100, Ca_err=None, PCO2_err=None, molL=False)

Calculates the calcite/limestone dissolution rate from given calcium concentration and PCO2. Optionally uses Monte Carlo error propagation to calculate uncertainty in rates.

Parameters:

Ca : float, numpy.ndarray or pandas Series

Calcium concentration, default units are mg/L. Change to mol/L by setting keyword mol_L=true.

PCO2 : float, numpy.ndarray, or pandas Series

The partial pressure of CO2 for the solution(s).

T_C : float, numpy.ndarray, or pandas Series

The temperature of the water in degrees Celcius.

rho : float

Density of rock in g/cm^3. (default=2.6)

method : string

Determines method used to calculate dissolution rates. Set to either ‘PWP’ or ‘Palmer’. This keyword is required.

impure : boolean

Used when calculating Palmer rates. Determines whether to use the table values for impure calcite (True) or pure calcite (False). Impure calcite is more representative of typical limestone. (default = True)

per_tol : float

the fractional change in H concentration between iterations upon which the iteration is terminated (see solutionFromCaPCO2). default=0.001

error: boolean

Set to true if you want to use Monte Carlo Error propagation to estimate error in dissolution rate. Requires values for Ca_err and PCO2_err. (default=False)

error_num : integer

Size of random sample used in Monte Carlo Error propagation. default=100

Ca_err: float, numpy.ndarray, or pandas Series

Percent error in calcite concentration(s) (1=100%)

PCO2_err: float, numpy.ndarray, or pandas Series

Percent error in PCO2 values (1=100%)

molL : boolean

Are Ca units in mol/L. If so, set to true. Otherwise, units assumed are mg/L. (default=False, i.e. mg/L)

Returns:

R : float, numpy.ndarray, or pandas Series

calcite dissolution rate in mm/yr

R_err : float, numpy.ndarray, or pandas Series

error in dissolution rate (returned with R if keyword error=True)