Calculate Gross Biogeochemical Flux Rates from Isotope Pool Dilution Data
Retrieve default k fractionation value for a pool
Retrieve default P fractionation value for a pool
Cost function between observed and predicted pools
Estimate initial k from heavy isotope concentration data
Optimize production and consumption parameters for pool dilution data
Optimize production and consumption parameters for pool dilution data
Predict total pool, heavy isotope pool, and atom percent
Pool dilution is a isotope tracer technique wherein a biogeochemical pool is artifically enriched with its heavy isotopologue and the gross productive and consumptive fluxes of that pool are quantified by the change in pool size and isotopic composition over time. This package calculates gross production and consumption rates from closed-system isotopic pool dilution time series data. Pool size concentrations and heavy isotope (e.g., 15N) content are measured over time and the model optimizes production rate (P) and the first order rate constant (k) by minimizing error in the model-predicted total pool size, as well as the isotopic signature. The model optimizes rates by weighting information against the signal:noise ratio of concentration and heavy- isotope signatures using measurement precision as well as the magnitude of change over time. The calculations used here are based on von Fischer and Hedin (2002) <doi:10.1029/2001GB001448> with some modifications.