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, confidence=0.0, return_samples=False)[source]

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

Parameters:
Cafloat, numpy.ndarray or pandas Series

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

PCO2float, numpy.ndarray, or pandas Series

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

T_Cfloat, numpy.ndarray, or pandas Series

The temperature of the water in degrees Celcius.

rhofloat

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

methodstring

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

impureboolean

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_tolfloat

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_numinteger

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%)

molLboolean

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

confidencefloat

If non-zero then confidence intervals will be used in error estimation (e.g. 90 = 90% confidence). Default is 0.

return_samplesboolean

If true (default is false), then return entire random samples within error arrays.

Returns:
Rfloat, numpy.ndarray, or pandas Series

calcite dissolution rate in mm/yr

R_errfloat, numpy.ndarray, or pandas Series

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