{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Quebrada Sonadora Nutrient Exports\n",
"---\n",
"Fluxes were calculated using the “loadflex” program (Appling et al. 2015), which is an enhancement of the widely used USGS “LOADEST” model (Runkel et al. 2004) as implemented in R. Both loadflex and LOADEST rely primarily on relationships between concentration of a solute and instantaneous discharge at the time of sampling to estimate concentrations when measured values of river chemistry are not available. In practice, and in our case, this often means that the 15-minute record of discharge at a given station is paired with an estimate of concentration based on weekly grab samples that span a range of discharge conditions. \n",
"\n",
"First install some packages we will use. "
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"# install and load packages, you only need to do this once. \n",
"\n",
"rprofile_path = file.path(Sys.getenv(\"HOME\"), \".Rprofile\")\n",
"write('\\noptions(repos=c(getOption(\\'repos\\'),\n",
" CRAN=\\'https://cloud.r-project.org\\',\n",
" USGS=\\'https://owi.usgs.gov/R\\'))\\n',\n",
" rprofile_path, \n",
" append = TRUE)\n",
"install.packages('smwrBase')\n",
"install.packages('rloadest')\n",
"install.packages('car')\n",
"install.packages('unitted')\n",
"install.packages('dplyr')\n",
"packageurl <- \"https://github.com/USGS-R/loadflex/archive/v1.0.1.tar.gz\"\n",
"install.packages(packageurl, repos=NULL, type=\"source\")\n",
"install.packages('dplyr')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"load some r packages we will use."
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"options(warn=-1)\n",
"\n",
"library('loadflex')\n",
"library('rloadest')"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"You will need to include the data directory from this dataset in the same location you are running this notebook. \n",
"\n",
"Initially here we will load discharge and weekly grab sample data from csv files. This notebook details nutrient fluxes for the Quebrada Sonadora sampling site which has a watershed area of 261.58 hectares. "
]
},
{
"cell_type": "code",
"execution_count": 1,
"metadata": {},
"outputs": [],
"source": [
"QS_chem <- readRDS(\"data/QS_chem_UNH.rds\")\n",
"QSDischargeShort <- readRDS(\"data/QS_Discharge_USGS.rds\")\n",
"QSWatershedArea <- 261.5888\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Calculate flux totals and mean annual concentration for Na. "
]
},
{
"cell_type": "code",
"execution_count": 2,
"metadata": {},
"outputs": [
{
"ename": "ERROR",
"evalue": "",
"output_type": "error",
"traceback": [
""
]
}
],
"source": [
"QS_chemNa <- QS_chem[complete.cases(QS_chem[\"Na_mg_L\"]),]\n",
"QS_chemNa <- QS_chemNa[complete.cases(QS_chemNa[\"CFS\"]),]\n",
"\n",
"\n",
"meta <- metadata(constituent=\"Na_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(Na_mg_L ~ model(1), data=QS_chemNa,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemNa, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"# preds_lrNa <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"preds_lcNa <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcNa <- aggregateSolute(preds_lcNa,meta, format=\"flux total\", se.preds=preds_lcNa$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemNa, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcNaConc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcNaConc <- aggregateSolute(preds_lcNaConc,meta, format=\"conc\", se.preds=preds_lcNaConc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"aggs_lcNa$Na_Kg_Ha_yr <- aggs_lcNa$Flux_Total / QSWatershedArea\n",
"aggs_lcNa$Na_mg_l <- aggs_lcNaConc$Conc\n",
"FluxTotalsdf <- aggs_lcNa[, c('Na_mg_l','Na_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"#round values \n",
"years <- FluxTotalsdf$Calendar_Year\n",
"FluxTotalsdf$Calendar_Year <- NULL\n",
"FluxTotalsdf <- round(FluxTotalsdf, 2)\n",
"FluxTotalsdf$Calendar_Year <- years\n"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Calculate flux totals and mean annual concentration for Ca. "
]
},
{
"cell_type": "code",
"execution_count": 3,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(Ca_mg_L ~ model(1), data = QS_chemCa, flow = \"CFS\", dates = \"date\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
}
],
"source": [
"QS_chemCa <- QS_chem[complete.cases(QS_chem[\"Ca_mg_L\"]),]\n",
"QS_chemCa <- QS_chemCa[complete.cases(QS_chemCa[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"Ca_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"Rio Piedras, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(Ca_mg_L ~ model(1), data=QS_chemCa,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemCa, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrCa <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrCa <- aggregateSolute(preds_lrCa,meta, format=\"flux total\", se.preds=preds_lrCa$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemCa, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcCaConc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcCaConc <- aggregateSolute(preds_lcCaConc,meta, format=\"conc\", se.preds=preds_lcCaConc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"aggs_lrCa$Ca_Kg_Ha_yr <- aggs_lrCa$Flux_Total / QSWatershedArea\n",
"aggs_lrCa$Ca_mg_l <- aggs_lcCaConc$Conc\n",
"FluxTotalsdfCa <- aggs_lrCa[, c('Ca_mg_l','Ca_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfCa$Calendar_Year\n",
"FluxTotalsdfCa$Calendar_Year <- NULL\n",
"FluxTotalsdfCa <- round(FluxTotalsdfCa, 2)\n",
"FluxTotalsdfCa$Calendar_Year <- years\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfCa,FluxTotalsdf, by = \"Calendar_Year\")\n"
]
},
{
"cell_type": "code",
"execution_count": 4,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"
\n",
"Calendar_Year | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllll}\n",
" Calendar\\_Year & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year Ca_mg_l Ca_Kg_Ha_yr Na_mg_l Na_Kg_Ha_yr\n",
"1 2009 2.70 56.91 4.76 109.37 \n",
"2 2010 2.61 108.81 4.50 193.90 \n",
"3 2011 2.83 129.89 4.66 225.30 \n",
"4 2012 2.67 103.36 4.64 179.82 \n",
"5 2013 2.52 93.10 4.60 171.38 \n",
"6 2014 2.92 75.59 4.91 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"# FluxTtotalsdf <- merge(FluxTtotalsdfCa,FluxTtotalsdf, by.x='Calendar_Year', by.y='Calendar_Year')\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 5,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(Mg_mg_L ~ model(1), data = QS_chemMg, flow = \"CFS\", dates = \"date\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllll}\n",
" Calendar\\_Year & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year Mg_mg_l Mg_Kg_Ha_yr Ca_mg_l Ca_Kg_Ha_yr Na_mg_l Na_Kg_Ha_yr\n",
"1 2009 1.69 34.81 2.70 56.91 4.76 109.37 \n",
"2 2010 1.39 59.00 2.61 108.81 4.50 193.90 \n",
"3 2011 1.40 62.26 2.83 129.89 4.66 225.30 \n",
"4 2012 1.46 54.64 2.67 103.36 4.64 179.82 \n",
"5 2013 1.58 54.93 2.52 93.10 4.60 171.38 \n",
"6 2014 2.39 38.81 2.92 75.59 4.91 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"#Mg\n",
"QS_chemMg <- QS_chem[complete.cases(QS_chem[\"Mg_mg_L\"]),]\n",
"QS_chemMg <- QS_chemMg[complete.cases(QS_chemMg[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"Mg_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(Mg_mg_L ~ model(1), data=QS_chemMg,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemMg, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrMg <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrMg <- aggregateSolute(preds_lrMg,meta, format=\"flux total\", se.preds=preds_lrMg$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemMg, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcMgConc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcMgConc <- aggregateSolute(preds_lcMgConc,meta, format=\"conc\", se.preds=preds_lcMgConc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"aggs_lrMg$Mg_Kg_Ha_yr <- aggs_lrMg$Flux_Total / QSWatershedArea\n",
"aggs_lrMg$Mg_mg_l <- aggs_lcMgConc$Conc\n",
"FluxTotalsdfMg <- aggs_lrMg[, c('Mg_mg_l','Mg_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfMg$Calendar_Year\n",
"FluxTotalsdfMg$Calendar_Year <- NULL\n",
"FluxTotalsdfMg <- round(FluxTotalsdfMg, 2)\n",
"FluxTotalsdfMg$Calendar_Year <- years\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfMg,FluxTotals, by = \"Calendar_Year\")\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 8,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(NO3.N_ug_N_L ~ model(1), data = QS_chemNO3, flow = \"CFS\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | NO3_mg_l.x | NO3_Kg_Ha_yr.x | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 0.067 | 1.033 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 0.061 | 1.798 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 0.056 | 1.873 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 0.042 | 1.193 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 0.037 | 0.903 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 0.055 | 1.115 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllllll}\n",
" Calendar\\_Year & NO3\\_mg\\_l.x & NO3\\_Kg\\_Ha\\_yr.x & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l.y & NO3\\_Kg\\_Ha\\_yr.y & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 0.067 & 1.033 & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 0.061 & 1.798 & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 0.056 & 1.873 & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 0.042 & 1.193 & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 0.037 & 0.903 & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 0.055 & 1.115 & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | NO3_mg_l.x | NO3_Kg_Ha_yr.x | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 0.067 | 1.033 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 0.061 | 1.798 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 0.056 | 1.873 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 0.042 | 1.193 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 0.037 | 0.903 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 0.055 | 1.115 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year NO3_mg_l.x NO3_Kg_Ha_yr.x SO4_mg_l SO4_Kg_Ha_yr NO3_mg_l.y\n",
"1 2009 0.067 1.033 0.55 15.21 0.07 \n",
"2 2010 0.061 1.798 0.53 25.03 0.06 \n",
"3 2011 0.056 1.873 0.52 25.99 0.06 \n",
"4 2012 0.042 1.193 0.52 20.92 0.04 \n",
"5 2013 0.037 0.903 0.51 20.67 0.04 \n",
"6 2014 0.055 1.115 0.54 15.46 0.05 \n",
" NO3_Kg_Ha_yr.y Mg_mg_l Mg_Kg_Ha_yr Ca_mg_l Ca_Kg_Ha_yr Na_mg_l Na_Kg_Ha_yr\n",
"1 1.03 1.69 34.81 2.70 56.91 4.76 109.37 \n",
"2 1.80 1.39 59.00 2.61 108.81 4.50 193.90 \n",
"3 1.87 1.40 62.26 2.83 129.89 4.66 225.30 \n",
"4 1.19 1.46 54.64 2.67 103.36 4.64 179.82 \n",
"5 0.90 1.58 54.93 2.52 93.10 4.60 171.38 \n",
"6 1.11 2.39 38.81 2.92 75.59 4.91 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"QS_chemNO3 <- QS_chem[complete.cases(QS_chem[\"NO3.N_ug_N_L\"]),]\n",
"QS_chemNO3 <- QS_chemNO3[complete.cases(QS_chemNO3[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"NO3.N_ug_N_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"ug/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(NO3.N_ug_N_L ~ model(1), data=QS_chemNO3,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"ug/L\", load.units=\"kg\"))\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemNO3, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcNO3Conc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcNO3Conc <- aggregateSolute(preds_lcNO3Conc,meta, format=\"conc\", se.preds=preds_lcNO3Conc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemNO3, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrNO3 <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrNO3 <- aggregateSolute(preds_lrNO3,meta, format=\"flux total\", se.preds=preds_lrNO3$se.pred, agg.by=\"calendar year\")\n",
"aggs_lrNO3$NO3_Kg_Ha_yr <- aggs_lrNO3$Flux_Total / QSWatershedArea\n",
"aggs_lrNO3$NO3_mg_l <- aggs_lcNO3Conc$Conc / 1000.0\n",
"FluxTotalsdfNO3 <- aggs_lrNO3[, c('NO3_mg_l','NO3_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfNO3$Calendar_Year\n",
"FluxTotalsdfNO3$Calendar_Year <- NULL\n",
"FluxTotalsdfNO3 <- round(FluxTotalsdfNO3, 3)\n",
"FluxTotalsdfNO3$Calendar_Year <- years\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfNO3,FluxTotals, by = \"Calendar_Year\")\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 7,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(SO4_S_mg_L ~ model(1), data = QS_chemSO4, flow = \"CFS\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l | NO3_Kg_Ha_yr | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllll}\n",
" Calendar\\_Year & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l & NO3\\_Kg\\_Ha\\_yr & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l | NO3_Kg_Ha_yr | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year SO4_mg_l SO4_Kg_Ha_yr NO3_mg_l NO3_Kg_Ha_yr Mg_mg_l Mg_Kg_Ha_yr\n",
"1 2009 0.55 15.21 0.07 1.03 1.69 34.81 \n",
"2 2010 0.53 25.03 0.06 1.80 1.39 59.00 \n",
"3 2011 0.52 25.99 0.06 1.87 1.40 62.26 \n",
"4 2012 0.52 20.92 0.04 1.19 1.46 54.64 \n",
"5 2013 0.51 20.67 0.04 0.90 1.58 54.93 \n",
"6 2014 0.54 15.46 0.05 1.11 2.39 38.81 \n",
" Ca_mg_l Ca_Kg_Ha_yr Na_mg_l Na_Kg_Ha_yr\n",
"1 2.70 56.91 4.76 109.37 \n",
"2 2.61 108.81 4.50 193.90 \n",
"3 2.83 129.89 4.66 225.30 \n",
"4 2.67 103.36 4.64 179.82 \n",
"5 2.52 93.10 4.60 171.38 \n",
"6 2.92 75.59 4.91 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"#SO4\n",
"QS_chemSO4 <- QS_chem[complete.cases(QS_chem[\"SO4_S_mg_L\"]),]\n",
"QS_chemSO4 <- QS_chemSO4[complete.cases(QS_chemSO4[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"SO4_S_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(SO4_S_mg_L ~ model(1), data=QS_chemSO4,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemSO4, abs.or.rel.resids=\"absolute\", interp.function=triangularInterpolation)\n",
"\n",
"preds_lcSO4Conc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcSO4Conc <- aggregateSolute(preds_lcSO4Conc,meta, format=\"conc\", se.preds=preds_lcSO4Conc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemSO4, abs.or.rel.resids=\"relative\", interp.function=triangularInterpolation)\n",
"\n",
"preds_lrSO4 <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrSO4 <- aggregateSolute(preds_lrSO4,meta, format=\"flux total\", se.preds=preds_lrSO4$se.pred, agg.by=\"calendar year\")\n",
"\n",
"\n",
"aggs_lrSO4$SO4_Kg_Ha_yr <- aggs_lrSO4$Flux_Total / QSWatershedArea\n",
"aggs_lrSO4$SO4_mg_l <- aggs_lcSO4Conc$Conc\n",
"FluxTotalsdfSO4 <- aggs_lrSO4[, c('SO4_mg_l','SO4_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfSO4$Calendar_Year\n",
"FluxTotalsdfSO4$Calendar_Year <- NULL\n",
"FluxTotalsdfSO4 <- round(FluxTotalsdfSO4, 3)\n",
"FluxTotalsdfSO4$Calendar_Year <- years\n",
"\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfSO4,FluxTotals, by = \"Calendar_Year\")\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 12,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
"Calendar_Year | F_mg_l | F_Kg_Ha_yr | NO3_mg_l | NO3_Kg_Ha_yr | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2013 | 0.017 | 0.589 | 0.037 | 0.903 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 0.019 | 0.452 | 0.055 | 1.115 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\t2009 | NA | NA | 0.067 | 1.033 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | NA | NA | 0.061 | 1.798 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | NA | NA | 0.056 | 1.873 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | NA | NA | 0.042 | 1.193 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllllllll}\n",
" Calendar\\_Year & F\\_mg\\_l & F\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l & NO3\\_Kg\\_Ha\\_yr & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l.y & NO3\\_Kg\\_Ha\\_yr.y & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2013 & 0.017 & 0.589 & 0.037 & 0.903 & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 0.019 & 0.452 & 0.055 & 1.115 & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\t 2009 & NA & NA & 0.067 & 1.033 & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & NA & NA & 0.061 & 1.798 & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & NA & NA & 0.056 & 1.873 & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & NA & NA & 0.042 & 1.193 & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | F_mg_l | F_Kg_Ha_yr | NO3_mg_l | NO3_Kg_Ha_yr | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2013 | 0.017 | 0.589 | 0.037 | 0.903 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 0.019 | 0.452 | 0.055 | 1.115 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"| 2009 | NA | NA | 0.067 | 1.033 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | NA | NA | 0.061 | 1.798 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | NA | NA | 0.056 | 1.873 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | NA | NA | 0.042 | 1.193 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year F_mg_l F_Kg_Ha_yr NO3_mg_l NO3_Kg_Ha_yr SO4_mg_l SO4_Kg_Ha_yr\n",
"1 2013 0.017 0.589 0.037 0.903 0.51 20.67 \n",
"2 2014 0.019 0.452 0.055 1.115 0.54 15.46 \n",
"3 2009 NA NA 0.067 1.033 0.55 15.21 \n",
"4 2010 NA NA 0.061 1.798 0.53 25.03 \n",
"5 2011 NA NA 0.056 1.873 0.52 25.99 \n",
"6 2012 NA NA 0.042 1.193 0.52 20.92 \n",
" NO3_mg_l.y NO3_Kg_Ha_yr.y Mg_mg_l Mg_Kg_Ha_yr Ca_mg_l Ca_Kg_Ha_yr Na_mg_l\n",
"1 0.04 0.90 1.58 54.93 2.52 93.10 4.60 \n",
"2 0.05 1.11 2.39 38.81 2.92 75.59 4.91 \n",
"3 0.07 1.03 1.69 34.81 2.70 56.91 4.76 \n",
"4 0.06 1.80 1.39 59.00 2.61 108.81 4.50 \n",
"5 0.06 1.87 1.40 62.26 2.83 129.89 4.66 \n",
"6 0.04 1.19 1.46 54.64 2.67 103.36 4.64 \n",
" Na_Kg_Ha_yr\n",
"1 171.38 \n",
"2 134.12 \n",
"3 109.37 \n",
"4 193.90 \n",
"5 225.30 \n",
"6 179.82 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"#F\n",
"QS_chemF <- QS_chem[complete.cases(QS_chem[\"Fluoride_mg_L\"]),]\n",
"QS_chemF <- QS_chemF[complete.cases(QS_chemF[\"CFS\"]),]\n",
"QSDischargeShort2<- QSDischargeShort[QSDischargeShort$date >= as.POSIXct(\"2013-01-01\") \n",
" & QSDischargeShort$date <= as.POSIXct(\"2014-12-31\"),]\n",
"\n",
"meta <- metadata(constituent=\"Fluoride_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(Fluoride_mg_L ~ model(1), data=QS_chemF,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemF, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcFConc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort2, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcFConc <- aggregateSolute(preds_lcFConc,meta, format=\"conc\", se.preds=preds_lcFConc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemF, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrF <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort2, se.pred=TRUE, date=TRUE)\n",
"\n",
"#get some really high wierd values just remove them. \n",
"#preds_lrF <-preds_lrF[preds_lrF$fit < 50,]\n",
"\n",
"aggs_lrF <- aggregateSolute(preds_lrF,meta, format=\"flux total\", se.preds=preds_lrF$se.pred, agg.by=\"calendar year\")\n",
"aggs_lcFConc <- aggregateSolute(preds_lcFConc,meta, format=\"conc\", se.preds=preds_lcFConc$se.pred, agg.by=\"calendar year\")\n",
"aggs_lrF$F_Kg_Ha_yr <- aggs_lrF$Flux_Total / QSWatershedArea\n",
"aggs_lrF$F_mg_l <- aggs_lcFConc$Conc\n",
"FluxTotalsdfF <- aggs_lrF[, c('F_mg_l','F_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"# round values\n",
"years <- FluxTotalsdfF$Calendar_Year\n",
"FluxTotalsdfF$Calendar_Year <- NULL\n",
"FluxTotalsdfF <- round(FluxTotalsdfF, 3)\n",
"FluxTotalsdfF$Calendar_Year <- years\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfF,FluxTotals, by = \"Calendar_Year\", all.y=TRUE)\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": [
"QS_chemSiO2 <- QS_chem[complete.cases(QS_chem[\"SiO2_mg_L\"]),]\n",
"QS_chemSiO2 <- QS_chemSiO2[complete.cases(QS_chemSiO2[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"SiO2_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(SiO2_mg_L ~ model(1), data=QS_chemSiO2,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemSiO2, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcSiO2Conc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcSiO2Conc <- aggregateSolute(preds_lcSiO2Conc,meta, format=\"conc\", se.preds=preds_lcSiO2Conc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"\n",
"\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemSiO2, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrSiO2 <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrSiO2 <- aggregateSolute(preds_lrSiO2,meta, format=\"flux total\", se.preds=preds_lrSiO2$se.pred, agg.by=\"calendar year\")\n",
"\n",
"aggs_lrSiO2$SiO2_Kg_Ha_yr <- aggs_lrSiO2$Flux_Total / QSWatershedArea\n",
"aggs_lrSiO2$SiO2_mg_l <- aggs_lcSiO2Conc$Conc\n",
"FluxTotalsdfSiO2 <- aggs_lrSiO2[, c('SiO2_mg_l','SiO2_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfSiO2$Calendar_Year\n",
"FluxTotalsdfSiO2$Calendar_Year <- NULL\n",
"FluxTotalsdfSiO2 <- round(FluxTotalsdfSiO2, 2)\n",
"FluxTotalsdfSiO2$Calendar_Year <- years\n",
"\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfSiO2,FluxTotals, by = \"Calendar_Year\", all.y=TRUE)\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 13,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(PO4.P_ug_P_L ~ model(1), data = QS_chemPO4, flow = \"CFS\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | PO4_mg_l | PO4_Kg_Ha_yr | F_mg_l | F_Kg_Ha_yr | NO3_mg_l | NO3_Kg_Ha_yr | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 0.003 | 0.06 | NA | NA | 0.067 | 1.033 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 0.002 | 0.11 | NA | NA | 0.061 | 1.798 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 0.002 | 0.12 | NA | NA | 0.056 | 1.873 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 0.002 | 0.09 | NA | NA | 0.042 | 1.193 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 0.002 | 0.10 | 0.017 | 0.589 | 0.037 | 0.903 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 0.003 | 0.09 | 0.019 | 0.452 | 0.055 | 1.115 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllllllllll}\n",
" Calendar\\_Year & PO4\\_mg\\_l & PO4\\_Kg\\_Ha\\_yr & F\\_mg\\_l & F\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l & NO3\\_Kg\\_Ha\\_yr & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l.y & NO3\\_Kg\\_Ha\\_yr.y & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 0.003 & 0.06 & NA & NA & 0.067 & 1.033 & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 0.002 & 0.11 & NA & NA & 0.061 & 1.798 & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 0.002 & 0.12 & NA & NA & 0.056 & 1.873 & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 0.002 & 0.09 & NA & NA & 0.042 & 1.193 & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 0.002 & 0.10 & 0.017 & 0.589 & 0.037 & 0.903 & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 0.003 & 0.09 & 0.019 & 0.452 & 0.055 & 1.115 & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | PO4_mg_l | PO4_Kg_Ha_yr | F_mg_l | F_Kg_Ha_yr | NO3_mg_l | NO3_Kg_Ha_yr | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 0.003 | 0.06 | NA | NA | 0.067 | 1.033 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 0.002 | 0.11 | NA | NA | 0.061 | 1.798 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 0.002 | 0.12 | NA | NA | 0.056 | 1.873 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 0.002 | 0.09 | NA | NA | 0.042 | 1.193 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 0.002 | 0.10 | 0.017 | 0.589 | 0.037 | 0.903 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 0.003 | 0.09 | 0.019 | 0.452 | 0.055 | 1.115 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year PO4_mg_l PO4_Kg_Ha_yr F_mg_l F_Kg_Ha_yr NO3_mg_l NO3_Kg_Ha_yr\n",
"1 2009 0.003 0.06 NA NA 0.067 1.033 \n",
"2 2010 0.002 0.11 NA NA 0.061 1.798 \n",
"3 2011 0.002 0.12 NA NA 0.056 1.873 \n",
"4 2012 0.002 0.09 NA NA 0.042 1.193 \n",
"5 2013 0.002 0.10 0.017 0.589 0.037 0.903 \n",
"6 2014 0.003 0.09 0.019 0.452 0.055 1.115 \n",
" SO4_mg_l SO4_Kg_Ha_yr NO3_mg_l.y NO3_Kg_Ha_yr.y Mg_mg_l Mg_Kg_Ha_yr Ca_mg_l\n",
"1 0.55 15.21 0.07 1.03 1.69 34.81 2.70 \n",
"2 0.53 25.03 0.06 1.80 1.39 59.00 2.61 \n",
"3 0.52 25.99 0.06 1.87 1.40 62.26 2.83 \n",
"4 0.52 20.92 0.04 1.19 1.46 54.64 2.67 \n",
"5 0.51 20.67 0.04 0.90 1.58 54.93 2.52 \n",
"6 0.54 15.46 0.05 1.11 2.39 38.81 2.92 \n",
" Ca_Kg_Ha_yr Na_mg_l Na_Kg_Ha_yr\n",
"1 56.91 4.76 109.37 \n",
"2 108.81 4.50 193.90 \n",
"3 129.89 4.66 225.30 \n",
"4 103.36 4.64 179.82 \n",
"5 93.10 4.60 171.38 \n",
"6 75.59 4.91 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"#PO4\n",
"QS_chemPO4 <- QS_chem[complete.cases(QS_chem[\"PO4.P_ug_P_L\"]),]\n",
"QS_chemPO4 <- QS_chemPO4[complete.cases(QS_chemPO4[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"PO4.P_ug_P_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"ug/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(PO4.P_ug_P_L ~ model(1), data=QS_chemPO4,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"ug/L\", load.units=\"kg\"))\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemPO4, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcPO4Conc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcPO4Conc <- aggregateSolute(preds_lcPO4Conc,meta, format=\"conc\", se.preds=preds_lcPO4Conc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemPO4, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrPO4 <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrPO4 <- aggregateSolute(preds_lrPO4,meta, format=\"flux total\", se.preds=preds_lrPO4$se.pred, agg.by=\"calendar year\")\n",
"\n",
"aggs_lrPO4$PO4_Kg_Ha_yr <- aggs_lrPO4$Flux_Total / QSWatershedArea\n",
"aggs_lrPO4$PO4_mg_l <- aggs_lcPO4Conc$Conc / 1000.0\n",
"FluxTotalsdfPO4 <- aggs_lrPO4[, c('PO4_mg_l','PO4_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfPO4$Calendar_Year\n",
"FluxTotalsdfPO4$Calendar_Year <- NULL\n",
"FluxTotalsdfPO4$PO4_mg_l <- round(FluxTotalsdfPO4$PO4_mg_l, 3)\n",
"FluxTotalsdfPO4$PO4_Kg_Ha_yr <- round(FluxTotalsdfPO4$PO4_Kg_Ha_yr, 2)\n",
"\n",
"FluxTotalsdfPO4$Calendar_Year <- years\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfPO4,FluxTotals, by = \"Calendar_Year\", all.y=TRUE)\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 14,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(NH4.N_ug_L ~ model(1), data = QS_chemNH4, flow = \"CFS\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | NH4_mg_l | NH4_Kg_Ha_yr | PO4_mg_l | PO4_Kg_Ha_yr | F_mg_l | F_Kg_Ha_yr | NO3_mg_l | NO3_Kg_Ha_yr | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 5.801 | 0.17 | 0.003 | 0.06 | NA | NA | 0.067 | 1.033 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 8.741 | 0.43 | 0.002 | 0.11 | NA | NA | 0.061 | 1.798 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 4.176 | 0.22 | 0.002 | 0.12 | NA | NA | 0.056 | 1.873 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 4.791 | 0.17 | 0.002 | 0.09 | NA | NA | 0.042 | 1.193 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 5.029 | 0.19 | 0.002 | 0.10 | 0.017 | 0.589 | 0.037 | 0.903 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 5.423 | 0.19 | 0.003 | 0.09 | 0.019 | 0.452 | 0.055 | 1.115 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllllllllllll}\n",
" Calendar\\_Year & NH4\\_mg\\_l & NH4\\_Kg\\_Ha\\_yr & PO4\\_mg\\_l & PO4\\_Kg\\_Ha\\_yr & F\\_mg\\_l & F\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l & NO3\\_Kg\\_Ha\\_yr & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l.y & NO3\\_Kg\\_Ha\\_yr.y & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 5.801 & 0.17 & 0.003 & 0.06 & NA & NA & 0.067 & 1.033 & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 8.741 & 0.43 & 0.002 & 0.11 & NA & NA & 0.061 & 1.798 & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 4.176 & 0.22 & 0.002 & 0.12 & NA & NA & 0.056 & 1.873 & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 4.791 & 0.17 & 0.002 & 0.09 & NA & NA & 0.042 & 1.193 & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 5.029 & 0.19 & 0.002 & 0.10 & 0.017 & 0.589 & 0.037 & 0.903 & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 5.423 & 0.19 & 0.003 & 0.09 & 0.019 & 0.452 & 0.055 & 1.115 & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | NH4_mg_l | NH4_Kg_Ha_yr | PO4_mg_l | PO4_Kg_Ha_yr | F_mg_l | F_Kg_Ha_yr | NO3_mg_l | NO3_Kg_Ha_yr | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 5.801 | 0.17 | 0.003 | 0.06 | NA | NA | 0.067 | 1.033 | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 8.741 | 0.43 | 0.002 | 0.11 | NA | NA | 0.061 | 1.798 | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 4.176 | 0.22 | 0.002 | 0.12 | NA | NA | 0.056 | 1.873 | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 4.791 | 0.17 | 0.002 | 0.09 | NA | NA | 0.042 | 1.193 | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 5.029 | 0.19 | 0.002 | 0.10 | 0.017 | 0.589 | 0.037 | 0.903 | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 5.423 | 0.19 | 0.003 | 0.09 | 0.019 | 0.452 | 0.055 | 1.115 | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year NH4_mg_l NH4_Kg_Ha_yr PO4_mg_l PO4_Kg_Ha_yr F_mg_l F_Kg_Ha_yr\n",
"1 2009 5.801 0.17 0.003 0.06 NA NA \n",
"2 2010 8.741 0.43 0.002 0.11 NA NA \n",
"3 2011 4.176 0.22 0.002 0.12 NA NA \n",
"4 2012 4.791 0.17 0.002 0.09 NA NA \n",
"5 2013 5.029 0.19 0.002 0.10 0.017 0.589 \n",
"6 2014 5.423 0.19 0.003 0.09 0.019 0.452 \n",
" NO3_mg_l NO3_Kg_Ha_yr SO4_mg_l SO4_Kg_Ha_yr NO3_mg_l.y NO3_Kg_Ha_yr.y Mg_mg_l\n",
"1 0.067 1.033 0.55 15.21 0.07 1.03 1.69 \n",
"2 0.061 1.798 0.53 25.03 0.06 1.80 1.39 \n",
"3 0.056 1.873 0.52 25.99 0.06 1.87 1.40 \n",
"4 0.042 1.193 0.52 20.92 0.04 1.19 1.46 \n",
"5 0.037 0.903 0.51 20.67 0.04 0.90 1.58 \n",
"6 0.055 1.115 0.54 15.46 0.05 1.11 2.39 \n",
" Mg_Kg_Ha_yr Ca_mg_l Ca_Kg_Ha_yr Na_mg_l Na_Kg_Ha_yr\n",
"1 34.81 2.70 56.91 4.76 109.37 \n",
"2 59.00 2.61 108.81 4.50 193.90 \n",
"3 62.26 2.83 129.89 4.66 225.30 \n",
"4 54.64 2.67 103.36 4.64 179.82 \n",
"5 54.93 2.52 93.10 4.60 171.38 \n",
"6 38.81 2.92 75.59 4.91 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"#NH4 \n",
"\n",
"QS_chemNH4 <- QS_chem[complete.cases(QS_chem[\"NH4.N_ug_L\"]),]\n",
"QS_chemNH4 <- QS_chemNH4[complete.cases(QS_chemNH4[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"NH4.N_ug_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"ug/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(NH4.N_ug_L ~ model(1), data=QS_chemNH4,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"ug/L\", load.units=\"kg\"))\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemNH4, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcNH4Conc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcNH4Conc <- aggregateSolute(preds_lcNH4Conc,meta, format=\"conc\", se.preds=preds_lcNH4Conc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemNH4, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrNH4 <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrNH4 <- aggregateSolute(preds_lrNH4,meta, format=\"flux total\", se.preds=preds_lrNH4$se.pred, agg.by=\"calendar year\")\n",
"\n",
"aggs_lrNH4$NH4_Kg_Ha_yr <- aggs_lrNH4$Flux_Total / QSWatershedArea\n",
"aggs_lrNH4$NH4_mg_l <- aggs_lcNH4Conc$Conc\n",
"FluxTotalsdfNH4 <- aggs_lrNH4[, c('NH4_mg_l','NH4_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfNH4$Calendar_Year\n",
"FluxTotalsdfNH4$Calendar_Year <- NULL\n",
"FluxTotalsdfNH4$NH4_mg_l <- round(FluxTotalsdfNH4$NH4_mg_l, 3)\n",
"FluxTotalsdfNH4$NH4_Kg_Ha_yr <- round(FluxTotalsdfNH4$NH4_Kg_Ha_yr, 2)\n",
"FluxTotalsdfNH4$Calendar_Year <- years\n",
"\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfNH4,FluxTotals, by = \"Calendar_Year\", all.y=TRUE)\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 15,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(TDN_mg_L ~ model(1), data = QS_chemTDN, flow = \"CFS\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | TDN_mg_l | TDN_Kg_Ha_yr | NH4_mg_l | NH4_Kg_Ha_yr | PO4_mg_l | PO4_Kg_Ha_yr | F_mg_l | F_Kg_Ha_yr | NO3_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 0.14 | 3.16 | 5.801 | 0.17 | 0.003 | 0.06 | NA | NA | 0.067 | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 0.18 | 6.38 | 8.741 | 0.43 | 0.002 | 0.11 | NA | NA | 0.061 | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 0.14 | 6.53 | 4.176 | 0.22 | 0.002 | 0.12 | NA | NA | 0.056 | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 0.09 | 3.27 | 4.791 | 0.17 | 0.002 | 0.09 | NA | NA | 0.042 | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 0.10 | 3.47 | 5.029 | 0.19 | 0.002 | 0.10 | 0.017 | 0.589 | 0.037 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 0.12 | 3.36 | 5.423 | 0.19 | 0.003 | 0.09 | 0.019 | 0.452 | 0.055 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllllllllllllll}\n",
" Calendar\\_Year & TDN\\_mg\\_l & TDN\\_Kg\\_Ha\\_yr & NH4\\_mg\\_l & NH4\\_Kg\\_Ha\\_yr & PO4\\_mg\\_l & PO4\\_Kg\\_Ha\\_yr & F\\_mg\\_l & F\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l & ... & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l.y & NO3\\_Kg\\_Ha\\_yr.y & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 0.14 & 3.16 & 5.801 & 0.17 & 0.003 & 0.06 & NA & NA & 0.067 & ... & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 0.18 & 6.38 & 8.741 & 0.43 & 0.002 & 0.11 & NA & NA & 0.061 & ... & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 0.14 & 6.53 & 4.176 & 0.22 & 0.002 & 0.12 & NA & NA & 0.056 & ... & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 0.09 & 3.27 & 4.791 & 0.17 & 0.002 & 0.09 & NA & NA & 0.042 & ... & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 0.10 & 3.47 & 5.029 & 0.19 & 0.002 & 0.10 & 0.017 & 0.589 & 0.037 & ... & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 0.12 & 3.36 & 5.423 & 0.19 & 0.003 & 0.09 & 0.019 & 0.452 & 0.055 & ... & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | TDN_mg_l | TDN_Kg_Ha_yr | NH4_mg_l | NH4_Kg_Ha_yr | PO4_mg_l | PO4_Kg_Ha_yr | F_mg_l | F_Kg_Ha_yr | NO3_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 0.14 | 3.16 | 5.801 | 0.17 | 0.003 | 0.06 | NA | NA | 0.067 | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 0.18 | 6.38 | 8.741 | 0.43 | 0.002 | 0.11 | NA | NA | 0.061 | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 0.14 | 6.53 | 4.176 | 0.22 | 0.002 | 0.12 | NA | NA | 0.056 | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 0.09 | 3.27 | 4.791 | 0.17 | 0.002 | 0.09 | NA | NA | 0.042 | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 0.10 | 3.47 | 5.029 | 0.19 | 0.002 | 0.10 | 0.017 | 0.589 | 0.037 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 0.12 | 3.36 | 5.423 | 0.19 | 0.003 | 0.09 | 0.019 | 0.452 | 0.055 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year TDN_mg_l TDN_Kg_Ha_yr NH4_mg_l NH4_Kg_Ha_yr PO4_mg_l\n",
"1 2009 0.14 3.16 5.801 0.17 0.003 \n",
"2 2010 0.18 6.38 8.741 0.43 0.002 \n",
"3 2011 0.14 6.53 4.176 0.22 0.002 \n",
"4 2012 0.09 3.27 4.791 0.17 0.002 \n",
"5 2013 0.10 3.47 5.029 0.19 0.002 \n",
"6 2014 0.12 3.36 5.423 0.19 0.003 \n",
" PO4_Kg_Ha_yr F_mg_l F_Kg_Ha_yr NO3_mg_l ... SO4_mg_l SO4_Kg_Ha_yr NO3_mg_l.y\n",
"1 0.06 NA NA 0.067 ... 0.55 15.21 0.07 \n",
"2 0.11 NA NA 0.061 ... 0.53 25.03 0.06 \n",
"3 0.12 NA NA 0.056 ... 0.52 25.99 0.06 \n",
"4 0.09 NA NA 0.042 ... 0.52 20.92 0.04 \n",
"5 0.10 0.017 0.589 0.037 ... 0.51 20.67 0.04 \n",
"6 0.09 0.019 0.452 0.055 ... 0.54 15.46 0.05 \n",
" NO3_Kg_Ha_yr.y Mg_mg_l Mg_Kg_Ha_yr Ca_mg_l Ca_Kg_Ha_yr Na_mg_l Na_Kg_Ha_yr\n",
"1 1.03 1.69 34.81 2.70 56.91 4.76 109.37 \n",
"2 1.80 1.39 59.00 2.61 108.81 4.50 193.90 \n",
"3 1.87 1.40 62.26 2.83 129.89 4.66 225.30 \n",
"4 1.19 1.46 54.64 2.67 103.36 4.64 179.82 \n",
"5 0.90 1.58 54.93 2.52 93.10 4.60 171.38 \n",
"6 1.11 2.39 38.81 2.92 75.59 4.91 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"QS_chemTDN <- QS_chem[complete.cases(QS_chem[\"TDN_mg_L\"]),]\n",
"QS_chemTDN <- QS_chemTDN[complete.cases(QS_chemTDN[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"TDN_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(TDN_mg_L ~ model(1), data=QS_chemTDN,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemTDN, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrTDN <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrTDN <- aggregateSolute(preds_lrTDN,meta, format=\"flux total\", se.preds=preds_lrTDN$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemTDN, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcTDNConc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcTDNConc <- aggregateSolute(preds_lcTDNConc,meta, format=\"conc\", se.preds=preds_lcTDNConc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"aggs_lrTDN$TDN_Kg_Ha_yr <- aggs_lrTDN$Flux_Total / QSWatershedArea\n",
"aggs_lrTDN$TDN_mg_l <- aggs_lcTDNConc$Conc\n",
"FluxTotalsdfTDN <- aggs_lrTDN[, c('TDN_mg_l','TDN_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfTDN$Calendar_Year\n",
"FluxTotalsdfTDN$Calendar_Year <- NULL\n",
"FluxTotalsdfTDN <- round(FluxTotalsdfTDN, 2)\n",
"FluxTotalsdfTDN$Calendar_Year <- years\n",
"\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfTDN,FluxTotals, by = \"Calendar_Year\", all.y=TRUE)\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 16,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(DON_mg_L ~ model(1), data = QS_chemDON, flow = \"CFS\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | DON_mg_l | DON_Kg_Ha_yr | TDN_mg_l | TDN_Kg_Ha_yr | NH4_mg_l | NH4_Kg_Ha_yr | PO4_mg_l | PO4_Kg_Ha_yr | F_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 0.07 | 1.99 | 0.14 | 3.16 | 5.801 | 0.17 | 0.003 | 0.06 | NA | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 0.06 | 2.93 | 0.18 | 6.38 | 8.741 | 0.43 | 0.002 | 0.11 | NA | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 0.07 | 4.47 | 0.14 | 6.53 | 4.176 | 0.22 | 0.002 | 0.12 | NA | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 0.04 | 1.91 | 0.09 | 3.27 | 4.791 | 0.17 | 0.002 | 0.09 | NA | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 0.06 | 2.50 | 0.10 | 3.47 | 5.029 | 0.19 | 0.002 | 0.10 | 0.017 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 0.06 | 1.91 | 0.12 | 3.36 | 5.423 | 0.19 | 0.003 | 0.09 | 0.019 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllllllllllllllll}\n",
" Calendar\\_Year & DON\\_mg\\_l & DON\\_Kg\\_Ha\\_yr & TDN\\_mg\\_l & TDN\\_Kg\\_Ha\\_yr & NH4\\_mg\\_l & NH4\\_Kg\\_Ha\\_yr & PO4\\_mg\\_l & PO4\\_Kg\\_Ha\\_yr & F\\_mg\\_l & ... & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l.y & NO3\\_Kg\\_Ha\\_yr.y & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 0.07 & 1.99 & 0.14 & 3.16 & 5.801 & 0.17 & 0.003 & 0.06 & NA & ... & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 0.06 & 2.93 & 0.18 & 6.38 & 8.741 & 0.43 & 0.002 & 0.11 & NA & ... & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 0.07 & 4.47 & 0.14 & 6.53 & 4.176 & 0.22 & 0.002 & 0.12 & NA & ... & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 0.04 & 1.91 & 0.09 & 3.27 & 4.791 & 0.17 & 0.002 & 0.09 & NA & ... & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 0.06 & 2.50 & 0.10 & 3.47 & 5.029 & 0.19 & 0.002 & 0.10 & 0.017 & ... & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 0.06 & 1.91 & 0.12 & 3.36 & 5.423 & 0.19 & 0.003 & 0.09 & 0.019 & ... & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | DON_mg_l | DON_Kg_Ha_yr | TDN_mg_l | TDN_Kg_Ha_yr | NH4_mg_l | NH4_Kg_Ha_yr | PO4_mg_l | PO4_Kg_Ha_yr | F_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 0.07 | 1.99 | 0.14 | 3.16 | 5.801 | 0.17 | 0.003 | 0.06 | NA | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 0.06 | 2.93 | 0.18 | 6.38 | 8.741 | 0.43 | 0.002 | 0.11 | NA | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 0.07 | 4.47 | 0.14 | 6.53 | 4.176 | 0.22 | 0.002 | 0.12 | NA | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 0.04 | 1.91 | 0.09 | 3.27 | 4.791 | 0.17 | 0.002 | 0.09 | NA | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 0.06 | 2.50 | 0.10 | 3.47 | 5.029 | 0.19 | 0.002 | 0.10 | 0.017 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 0.06 | 1.91 | 0.12 | 3.36 | 5.423 | 0.19 | 0.003 | 0.09 | 0.019 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year DON_mg_l DON_Kg_Ha_yr TDN_mg_l TDN_Kg_Ha_yr NH4_mg_l\n",
"1 2009 0.07 1.99 0.14 3.16 5.801 \n",
"2 2010 0.06 2.93 0.18 6.38 8.741 \n",
"3 2011 0.07 4.47 0.14 6.53 4.176 \n",
"4 2012 0.04 1.91 0.09 3.27 4.791 \n",
"5 2013 0.06 2.50 0.10 3.47 5.029 \n",
"6 2014 0.06 1.91 0.12 3.36 5.423 \n",
" NH4_Kg_Ha_yr PO4_mg_l PO4_Kg_Ha_yr F_mg_l ... SO4_mg_l SO4_Kg_Ha_yr\n",
"1 0.17 0.003 0.06 NA ... 0.55 15.21 \n",
"2 0.43 0.002 0.11 NA ... 0.53 25.03 \n",
"3 0.22 0.002 0.12 NA ... 0.52 25.99 \n",
"4 0.17 0.002 0.09 NA ... 0.52 20.92 \n",
"5 0.19 0.002 0.10 0.017 ... 0.51 20.67 \n",
"6 0.19 0.003 0.09 0.019 ... 0.54 15.46 \n",
" NO3_mg_l.y NO3_Kg_Ha_yr.y Mg_mg_l Mg_Kg_Ha_yr Ca_mg_l Ca_Kg_Ha_yr Na_mg_l\n",
"1 0.07 1.03 1.69 34.81 2.70 56.91 4.76 \n",
"2 0.06 1.80 1.39 59.00 2.61 108.81 4.50 \n",
"3 0.06 1.87 1.40 62.26 2.83 129.89 4.66 \n",
"4 0.04 1.19 1.46 54.64 2.67 103.36 4.64 \n",
"5 0.04 0.90 1.58 54.93 2.52 93.10 4.60 \n",
"6 0.05 1.11 2.39 38.81 2.92 75.59 4.91 \n",
" Na_Kg_Ha_yr\n",
"1 109.37 \n",
"2 193.90 \n",
"3 225.30 \n",
"4 179.82 \n",
"5 171.38 \n",
"6 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"#DON\n",
"QS_chemDON <- QS_chem[complete.cases(QS_chem[\"DON_mg_L\"]),]\n",
"QS_chemDON <- QS_chemDON[complete.cases(QS_chemDON[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"DON_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(DON_mg_L ~ model(1), data=QS_chemDON,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemDON, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrDON <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrDON <- aggregateSolute(preds_lrDON,meta, format=\"flux total\", se.preds=preds_lrDON$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemDON, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcDONConc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcDONConc <- aggregateSolute(preds_lcDONConc,meta, format=\"conc\", se.preds=preds_lcDONConc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"\n",
"aggs_lrDON$DON_Kg_Ha_yr <- aggs_lrDON$Flux_Total / QSWatershedArea\n",
"aggs_lrDON$DON_mg_l <- aggs_lcDONConc$Conc\n",
"FluxTotalsdfDON <- aggs_lrDON[, c('DON_mg_l','DON_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfDON$Calendar_Year\n",
"FluxTotalsdfDON$Calendar_Year <- NULL\n",
"FluxTotalsdfDON <- round(FluxTotalsdfDON, 2)\n",
"FluxTotalsdfDON$Calendar_Year <- years\n",
"\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfDON,FluxTotals, by = \"Calendar_Year\", all.y=TRUE)\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 17,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(DOC_mg_L ~ model(1), data = QS_chemDOC, flow = \"CFS\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | DOC_mg_l | DOC_Kg_Ha_yr | DON_mg_l | DON_Kg_Ha_yr | TDN_mg_l | TDN_Kg_Ha_yr | NH4_mg_l | NH4_Kg_Ha_yr | PO4_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 2.22 | 98.32 | 0.07 | 1.99 | 0.14 | 3.16 | 5.801 | 0.17 | 0.003 | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 2.05 | 159.90 | 0.06 | 2.93 | 0.18 | 6.38 | 8.741 | 0.43 | 0.002 | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 1.95 | 152.43 | 0.07 | 4.47 | 0.14 | 6.53 | 4.176 | 0.22 | 0.002 | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 1.74 | 97.33 | 0.04 | 1.91 | 0.09 | 3.27 | 4.791 | 0.17 | 0.002 | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 1.76 | 143.95 | 0.06 | 2.50 | 0.10 | 3.47 | 5.029 | 0.19 | 0.002 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 1.61 | 62.29 | 0.06 | 1.91 | 0.12 | 3.36 | 5.423 | 0.19 | 0.003 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllllllllllllllllll}\n",
" Calendar\\_Year & DOC\\_mg\\_l & DOC\\_Kg\\_Ha\\_yr & DON\\_mg\\_l & DON\\_Kg\\_Ha\\_yr & TDN\\_mg\\_l & TDN\\_Kg\\_Ha\\_yr & NH4\\_mg\\_l & NH4\\_Kg\\_Ha\\_yr & PO4\\_mg\\_l & ... & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l.y & NO3\\_Kg\\_Ha\\_yr.y & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 2.22 & 98.32 & 0.07 & 1.99 & 0.14 & 3.16 & 5.801 & 0.17 & 0.003 & ... & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 2.05 & 159.90 & 0.06 & 2.93 & 0.18 & 6.38 & 8.741 & 0.43 & 0.002 & ... & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 1.95 & 152.43 & 0.07 & 4.47 & 0.14 & 6.53 & 4.176 & 0.22 & 0.002 & ... & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 1.74 & 97.33 & 0.04 & 1.91 & 0.09 & 3.27 & 4.791 & 0.17 & 0.002 & ... & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 1.76 & 143.95 & 0.06 & 2.50 & 0.10 & 3.47 & 5.029 & 0.19 & 0.002 & ... & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 1.61 & 62.29 & 0.06 & 1.91 & 0.12 & 3.36 & 5.423 & 0.19 & 0.003 & ... & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | DOC_mg_l | DOC_Kg_Ha_yr | DON_mg_l | DON_Kg_Ha_yr | TDN_mg_l | TDN_Kg_Ha_yr | NH4_mg_l | NH4_Kg_Ha_yr | PO4_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 2.22 | 98.32 | 0.07 | 1.99 | 0.14 | 3.16 | 5.801 | 0.17 | 0.003 | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 2.05 | 159.90 | 0.06 | 2.93 | 0.18 | 6.38 | 8.741 | 0.43 | 0.002 | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 1.95 | 152.43 | 0.07 | 4.47 | 0.14 | 6.53 | 4.176 | 0.22 | 0.002 | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 1.74 | 97.33 | 0.04 | 1.91 | 0.09 | 3.27 | 4.791 | 0.17 | 0.002 | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 1.76 | 143.95 | 0.06 | 2.50 | 0.10 | 3.47 | 5.029 | 0.19 | 0.002 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 1.61 | 62.29 | 0.06 | 1.91 | 0.12 | 3.36 | 5.423 | 0.19 | 0.003 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year DOC_mg_l DOC_Kg_Ha_yr DON_mg_l DON_Kg_Ha_yr TDN_mg_l\n",
"1 2009 2.22 98.32 0.07 1.99 0.14 \n",
"2 2010 2.05 159.90 0.06 2.93 0.18 \n",
"3 2011 1.95 152.43 0.07 4.47 0.14 \n",
"4 2012 1.74 97.33 0.04 1.91 0.09 \n",
"5 2013 1.76 143.95 0.06 2.50 0.10 \n",
"6 2014 1.61 62.29 0.06 1.91 0.12 \n",
" TDN_Kg_Ha_yr NH4_mg_l NH4_Kg_Ha_yr PO4_mg_l ... SO4_mg_l SO4_Kg_Ha_yr\n",
"1 3.16 5.801 0.17 0.003 ... 0.55 15.21 \n",
"2 6.38 8.741 0.43 0.002 ... 0.53 25.03 \n",
"3 6.53 4.176 0.22 0.002 ... 0.52 25.99 \n",
"4 3.27 4.791 0.17 0.002 ... 0.52 20.92 \n",
"5 3.47 5.029 0.19 0.002 ... 0.51 20.67 \n",
"6 3.36 5.423 0.19 0.003 ... 0.54 15.46 \n",
" NO3_mg_l.y NO3_Kg_Ha_yr.y Mg_mg_l Mg_Kg_Ha_yr Ca_mg_l Ca_Kg_Ha_yr Na_mg_l\n",
"1 0.07 1.03 1.69 34.81 2.70 56.91 4.76 \n",
"2 0.06 1.80 1.39 59.00 2.61 108.81 4.50 \n",
"3 0.06 1.87 1.40 62.26 2.83 129.89 4.66 \n",
"4 0.04 1.19 1.46 54.64 2.67 103.36 4.64 \n",
"5 0.04 0.90 1.58 54.93 2.52 93.10 4.60 \n",
"6 0.05 1.11 2.39 38.81 2.92 75.59 4.91 \n",
" Na_Kg_Ha_yr\n",
"1 109.37 \n",
"2 193.90 \n",
"3 225.30 \n",
"4 179.82 \n",
"5 171.38 \n",
"6 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"#DOC\n",
"QS_chemDOC <- QS_chem[complete.cases(QS_chem[\"DOC_mg_L\"]),]\n",
"QS_chemDOC <- QS_chemDOC[complete.cases(QS_chemDOC[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"DOC_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(DOC_mg_L ~ model(1), data=QS_chemDOC,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemDOC, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrDOC <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrDOC <- aggregateSolute(preds_lrDOC,meta, format=\"flux total\", se.preds=preds_lrDOC$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemDOC, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcDOCConc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcDOCConc <- aggregateSolute(preds_lcDOCConc,meta, format=\"conc\", se.preds=preds_lcDOCConc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"aggs_lrDOC$DOC_Kg_Ha_yr <- aggs_lrDOC$Flux_Total / QSWatershedArea\n",
"aggs_lrDOC$DOC_mg_l <- aggs_lcDOCConc$Conc\n",
"FluxTotalsdfDOC <- aggs_lrDOC[, c('DOC_mg_l','DOC_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfDOC$Calendar_Year\n",
"FluxTotalsdfDOC$Calendar_Year <- NULL\n",
"FluxTotalsdfDOC <- round(FluxTotalsdfDOC, 2)\n",
"FluxTotalsdfDOC$Calendar_Year <- years\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfDOC,FluxTotals, by = \"Calendar_Year\", all.y=TRUE)\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 18,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(Cl_mg_L ~ model(1), data = QS_chemCl, flow = \"CFS\", dates = \"date\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | Cl_mg_l | Cl_Kg_Ha_yr | DOC_mg_l | DOC_Kg_Ha_yr | DON_mg_l | DON_Kg_Ha_yr | TDN_mg_l | TDN_Kg_Ha_yr | NH4_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 7.22 | 181.60 | 2.22 | 98.32 | 0.07 | 1.99 | 0.14 | 3.16 | 5.801 | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 6.69 | 306.03 | 2.05 | 159.90 | 0.06 | 2.93 | 0.18 | 6.38 | 8.741 | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 6.72 | 327.98 | 1.95 | 152.43 | 0.07 | 4.47 | 0.14 | 6.53 | 4.176 | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 7.02 | 281.04 | 1.74 | 97.33 | 0.04 | 1.91 | 0.09 | 3.27 | 4.791 | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 7.25 | 290.83 | 1.76 | 143.95 | 0.06 | 2.50 | 0.10 | 3.47 | 5.029 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 7.38 | 213.48 | 1.61 | 62.29 | 0.06 | 1.91 | 0.12 | 3.36 | 5.423 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllllllllllllllllllll}\n",
" Calendar\\_Year & Cl\\_mg\\_l & Cl\\_Kg\\_Ha\\_yr & DOC\\_mg\\_l & DOC\\_Kg\\_Ha\\_yr & DON\\_mg\\_l & DON\\_Kg\\_Ha\\_yr & TDN\\_mg\\_l & TDN\\_Kg\\_Ha\\_yr & NH4\\_mg\\_l & ... & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l.y & NO3\\_Kg\\_Ha\\_yr.y & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 7.22 & 181.60 & 2.22 & 98.32 & 0.07 & 1.99 & 0.14 & 3.16 & 5.801 & ... & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 6.69 & 306.03 & 2.05 & 159.90 & 0.06 & 2.93 & 0.18 & 6.38 & 8.741 & ... & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 6.72 & 327.98 & 1.95 & 152.43 & 0.07 & 4.47 & 0.14 & 6.53 & 4.176 & ... & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 7.02 & 281.04 & 1.74 & 97.33 & 0.04 & 1.91 & 0.09 & 3.27 & 4.791 & ... & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 7.25 & 290.83 & 1.76 & 143.95 & 0.06 & 2.50 & 0.10 & 3.47 & 5.029 & ... & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 7.38 & 213.48 & 1.61 & 62.29 & 0.06 & 1.91 & 0.12 & 3.36 & 5.423 & ... & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | Cl_mg_l | Cl_Kg_Ha_yr | DOC_mg_l | DOC_Kg_Ha_yr | DON_mg_l | DON_Kg_Ha_yr | TDN_mg_l | TDN_Kg_Ha_yr | NH4_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 7.22 | 181.60 | 2.22 | 98.32 | 0.07 | 1.99 | 0.14 | 3.16 | 5.801 | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 6.69 | 306.03 | 2.05 | 159.90 | 0.06 | 2.93 | 0.18 | 6.38 | 8.741 | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 6.72 | 327.98 | 1.95 | 152.43 | 0.07 | 4.47 | 0.14 | 6.53 | 4.176 | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 7.02 | 281.04 | 1.74 | 97.33 | 0.04 | 1.91 | 0.09 | 3.27 | 4.791 | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 7.25 | 290.83 | 1.76 | 143.95 | 0.06 | 2.50 | 0.10 | 3.47 | 5.029 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 7.38 | 213.48 | 1.61 | 62.29 | 0.06 | 1.91 | 0.12 | 3.36 | 5.423 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year Cl_mg_l Cl_Kg_Ha_yr DOC_mg_l DOC_Kg_Ha_yr DON_mg_l DON_Kg_Ha_yr\n",
"1 2009 7.22 181.60 2.22 98.32 0.07 1.99 \n",
"2 2010 6.69 306.03 2.05 159.90 0.06 2.93 \n",
"3 2011 6.72 327.98 1.95 152.43 0.07 4.47 \n",
"4 2012 7.02 281.04 1.74 97.33 0.04 1.91 \n",
"5 2013 7.25 290.83 1.76 143.95 0.06 2.50 \n",
"6 2014 7.38 213.48 1.61 62.29 0.06 1.91 \n",
" TDN_mg_l TDN_Kg_Ha_yr NH4_mg_l ... SO4_mg_l SO4_Kg_Ha_yr NO3_mg_l.y\n",
"1 0.14 3.16 5.801 ... 0.55 15.21 0.07 \n",
"2 0.18 6.38 8.741 ... 0.53 25.03 0.06 \n",
"3 0.14 6.53 4.176 ... 0.52 25.99 0.06 \n",
"4 0.09 3.27 4.791 ... 0.52 20.92 0.04 \n",
"5 0.10 3.47 5.029 ... 0.51 20.67 0.04 \n",
"6 0.12 3.36 5.423 ... 0.54 15.46 0.05 \n",
" NO3_Kg_Ha_yr.y Mg_mg_l Mg_Kg_Ha_yr Ca_mg_l Ca_Kg_Ha_yr Na_mg_l Na_Kg_Ha_yr\n",
"1 1.03 1.69 34.81 2.70 56.91 4.76 109.37 \n",
"2 1.80 1.39 59.00 2.61 108.81 4.50 193.90 \n",
"3 1.87 1.40 62.26 2.83 129.89 4.66 225.30 \n",
"4 1.19 1.46 54.64 2.67 103.36 4.64 179.82 \n",
"5 0.90 1.58 54.93 2.52 93.10 4.60 171.38 \n",
"6 1.11 2.39 38.81 2.92 75.59 4.91 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"\n",
"#Cl\n",
"QS_chemCl <- QS_chem[complete.cases(QS_chem[\"Cl_mg_L\"]),]\n",
"QS_chemCl <- QS_chemCl[complete.cases(QS_chemCl[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"Cl_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(Cl_mg_L ~ model(1), data=QS_chemCl,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemCl, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrCl <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrCl <- aggregateSolute(preds_lrCl,meta, format=\"flux total\", se.preds=preds_lrCl$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemCl, abs.or.rel.resids=\"relative\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcClConc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lcClConc <- aggregateSolute(preds_lcClConc,meta, format=\"conc\", se.preds=preds_lcClConc$se.pred, agg.by=\"calendar year\")\n",
"\n",
"\n",
"aggs_lrCl$Cl_Kg_Ha_yr <- aggs_lrCl$Flux_Total / QSWatershedArea\n",
"aggs_lrCl$Cl_mg_l <- aggs_lcClConc$Conc\n",
"FluxTotalsdfCl <- aggs_lrCl[, c('Cl_mg_l','Cl_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfCl$Calendar_Year\n",
"FluxTotalsdfCl$Calendar_Year <- NULL\n",
"FluxTotalsdfCl <- round(FluxTotalsdfCl, 2)\n",
"FluxTotalsdfCl$Calendar_Year <- years\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfCl,FluxTotals, by = \"Calendar_Year\", all.y=TRUE)\n",
"FluxTotals"
]
},
{
"cell_type": "code",
"execution_count": 19,
"metadata": {},
"outputs": [
{
"name": "stderr",
"output_type": "stream",
"text": [
"Warning message in loadReg(K_mg_L ~ model(1), data = QS_chemK, flow = \"CFS\", dates = \"date\", :\n",
"\"The minimum spacing between observed loads is 6 days. The time between observations should be at least 7 days to avoid autocorrelation issues.\""
]
},
{
"data": {
"text/html": [
"\n",
"Calendar_Year | K_mg_l | K_Kg_Ha_yr | Cl_mg_l | Cl_Kg_Ha_yr | DOC_mg_l | DOC_Kg_Ha_yr | DON_mg_l | DON_Kg_Ha_yr | TDN_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr |
\n",
"\n",
"\t2009 | 0.29 | 6.99 | 7.22 | 181.60 | 2.22 | 98.32 | 0.07 | 1.99 | 0.14 | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 |
\n",
"\t2010 | 0.27 | 11.45 | 6.69 | 306.03 | 2.05 | 159.90 | 0.06 | 2.93 | 0.18 | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 |
\n",
"\t2011 | 0.21 | 10.49 | 6.72 | 327.98 | 1.95 | 152.43 | 0.07 | 4.47 | 0.14 | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 |
\n",
"\t2012 | 0.24 | 9.55 | 7.02 | 281.04 | 1.74 | 97.33 | 0.04 | 1.91 | 0.09 | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 |
\n",
"\t2013 | 0.23 | 8.80 | 7.25 | 290.83 | 1.76 | 143.95 | 0.06 | 2.50 | 0.10 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 |
\n",
"\t2014 | 0.30 | 8.55 | 7.38 | 213.48 | 1.61 | 62.29 | 0.06 | 1.91 | 0.12 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|lllllllllllllllllllllllllllll}\n",
" Calendar\\_Year & K\\_mg\\_l & K\\_Kg\\_Ha\\_yr & Cl\\_mg\\_l & Cl\\_Kg\\_Ha\\_yr & DOC\\_mg\\_l & DOC\\_Kg\\_Ha\\_yr & DON\\_mg\\_l & DON\\_Kg\\_Ha\\_yr & TDN\\_mg\\_l & ... & SO4\\_mg\\_l & SO4\\_Kg\\_Ha\\_yr & NO3\\_mg\\_l.y & NO3\\_Kg\\_Ha\\_yr.y & Mg\\_mg\\_l & Mg\\_Kg\\_Ha\\_yr & Ca\\_mg\\_l & Ca\\_Kg\\_Ha\\_yr & Na\\_mg\\_l & Na\\_Kg\\_Ha\\_yr\\\\\n",
"\\hline\n",
"\t 2009 & 0.29 & 6.99 & 7.22 & 181.60 & 2.22 & 98.32 & 0.07 & 1.99 & 0.14 & ... & 0.55 & 15.21 & 0.07 & 1.03 & 1.69 & 34.81 & 2.70 & 56.91 & 4.76 & 109.37\\\\\n",
"\t 2010 & 0.27 & 11.45 & 6.69 & 306.03 & 2.05 & 159.90 & 0.06 & 2.93 & 0.18 & ... & 0.53 & 25.03 & 0.06 & 1.80 & 1.39 & 59.00 & 2.61 & 108.81 & 4.50 & 193.90\\\\\n",
"\t 2011 & 0.21 & 10.49 & 6.72 & 327.98 & 1.95 & 152.43 & 0.07 & 4.47 & 0.14 & ... & 0.52 & 25.99 & 0.06 & 1.87 & 1.40 & 62.26 & 2.83 & 129.89 & 4.66 & 225.30\\\\\n",
"\t 2012 & 0.24 & 9.55 & 7.02 & 281.04 & 1.74 & 97.33 & 0.04 & 1.91 & 0.09 & ... & 0.52 & 20.92 & 0.04 & 1.19 & 1.46 & 54.64 & 2.67 & 103.36 & 4.64 & 179.82\\\\\n",
"\t 2013 & 0.23 & 8.80 & 7.25 & 290.83 & 1.76 & 143.95 & 0.06 & 2.50 & 0.10 & ... & 0.51 & 20.67 & 0.04 & 0.90 & 1.58 & 54.93 & 2.52 & 93.10 & 4.60 & 171.38\\\\\n",
"\t 2014 & 0.30 & 8.55 & 7.38 & 213.48 & 1.61 & 62.29 & 0.06 & 1.91 & 0.12 & ... & 0.54 & 15.46 & 0.05 & 1.11 & 2.39 & 38.81 & 2.92 & 75.59 & 4.91 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"Calendar_Year | K_mg_l | K_Kg_Ha_yr | Cl_mg_l | Cl_Kg_Ha_yr | DOC_mg_l | DOC_Kg_Ha_yr | DON_mg_l | DON_Kg_Ha_yr | TDN_mg_l | ... | SO4_mg_l | SO4_Kg_Ha_yr | NO3_mg_l.y | NO3_Kg_Ha_yr.y | Mg_mg_l | Mg_Kg_Ha_yr | Ca_mg_l | Ca_Kg_Ha_yr | Na_mg_l | Na_Kg_Ha_yr | \n",
"|---|---|---|---|---|---|\n",
"| 2009 | 0.29 | 6.99 | 7.22 | 181.60 | 2.22 | 98.32 | 0.07 | 1.99 | 0.14 | ... | 0.55 | 15.21 | 0.07 | 1.03 | 1.69 | 34.81 | 2.70 | 56.91 | 4.76 | 109.37 | \n",
"| 2010 | 0.27 | 11.45 | 6.69 | 306.03 | 2.05 | 159.90 | 0.06 | 2.93 | 0.18 | ... | 0.53 | 25.03 | 0.06 | 1.80 | 1.39 | 59.00 | 2.61 | 108.81 | 4.50 | 193.90 | \n",
"| 2011 | 0.21 | 10.49 | 6.72 | 327.98 | 1.95 | 152.43 | 0.07 | 4.47 | 0.14 | ... | 0.52 | 25.99 | 0.06 | 1.87 | 1.40 | 62.26 | 2.83 | 129.89 | 4.66 | 225.30 | \n",
"| 2012 | 0.24 | 9.55 | 7.02 | 281.04 | 1.74 | 97.33 | 0.04 | 1.91 | 0.09 | ... | 0.52 | 20.92 | 0.04 | 1.19 | 1.46 | 54.64 | 2.67 | 103.36 | 4.64 | 179.82 | \n",
"| 2013 | 0.23 | 8.80 | 7.25 | 290.83 | 1.76 | 143.95 | 0.06 | 2.50 | 0.10 | ... | 0.51 | 20.67 | 0.04 | 0.90 | 1.58 | 54.93 | 2.52 | 93.10 | 4.60 | 171.38 | \n",
"| 2014 | 0.30 | 8.55 | 7.38 | 213.48 | 1.61 | 62.29 | 0.06 | 1.91 | 0.12 | ... | 0.54 | 15.46 | 0.05 | 1.11 | 2.39 | 38.81 | 2.92 | 75.59 | 4.91 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" Calendar_Year K_mg_l K_Kg_Ha_yr Cl_mg_l Cl_Kg_Ha_yr DOC_mg_l DOC_Kg_Ha_yr\n",
"1 2009 0.29 6.99 7.22 181.60 2.22 98.32 \n",
"2 2010 0.27 11.45 6.69 306.03 2.05 159.90 \n",
"3 2011 0.21 10.49 6.72 327.98 1.95 152.43 \n",
"4 2012 0.24 9.55 7.02 281.04 1.74 97.33 \n",
"5 2013 0.23 8.80 7.25 290.83 1.76 143.95 \n",
"6 2014 0.30 8.55 7.38 213.48 1.61 62.29 \n",
" DON_mg_l DON_Kg_Ha_yr TDN_mg_l ... SO4_mg_l SO4_Kg_Ha_yr NO3_mg_l.y\n",
"1 0.07 1.99 0.14 ... 0.55 15.21 0.07 \n",
"2 0.06 2.93 0.18 ... 0.53 25.03 0.06 \n",
"3 0.07 4.47 0.14 ... 0.52 25.99 0.06 \n",
"4 0.04 1.91 0.09 ... 0.52 20.92 0.04 \n",
"5 0.06 2.50 0.10 ... 0.51 20.67 0.04 \n",
"6 0.06 1.91 0.12 ... 0.54 15.46 0.05 \n",
" NO3_Kg_Ha_yr.y Mg_mg_l Mg_Kg_Ha_yr Ca_mg_l Ca_Kg_Ha_yr Na_mg_l Na_Kg_Ha_yr\n",
"1 1.03 1.69 34.81 2.70 56.91 4.76 109.37 \n",
"2 1.80 1.39 59.00 2.61 108.81 4.50 193.90 \n",
"3 1.87 1.40 62.26 2.83 129.89 4.66 225.30 \n",
"4 1.19 1.46 54.64 2.67 103.36 4.64 179.82 \n",
"5 0.90 1.58 54.93 2.52 93.10 4.60 171.38 \n",
"6 1.11 2.39 38.81 2.92 75.59 4.91 134.12 "
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"QS_chemK <- QS_chem[complete.cases(QS_chem[\"K_mg_L\"]),]\n",
"QS_chemK <- QS_chemK[complete.cases(QS_chemK[\"CFS\"]),]\n",
"\n",
"meta <- metadata(constituent=\"K_mg_L\", flow=\"CFS\", \n",
" dates=\"date\", conc.units=\"mg/L\", flow.units=\"cfs\", load.units=\"kg\", \n",
" load.rate.units=\"kg d^-1\", station=\"QS, PR\")\n",
"\n",
"\n",
"lr <- loadReg2(loadReg(K_mg_L ~ model(1), data=QS_chemK,\n",
" flow=\"CFS\", dates=\"date\", time.step=\"instantaneous\", \n",
" flow.units=\"cfs\", conc.units=\"mg/L\", load.units=\"kg\"))\n",
"\n",
"lc <- loadComp(reg.model=lr, interp.format=\"flux\", \n",
" interp.data=QS_chemK, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lrK <- predictSolute(lc, \"flux\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"aggs_lrK <- aggregateSolute(preds_lrK,meta, format=\"flux total\", se.preds=preds_lrK$se.pred, agg.by=\"calendar year\")\n",
"\n",
"lcConc <- loadComp(reg.model=lr, interp.format=\"conc\", \n",
" interp.data=QS_chemK, abs.or.rel.resids=\"absolute\", interp.function=rectangularInterpolation)\n",
"\n",
"preds_lcKConc <- predictSolute(lcConc, \"conc\",interval=\"prediction\", QSDischargeShort, se.pred=TRUE, date=TRUE)\n",
"\n",
"\n",
"aggs_lrK$K_Kg_Ha_yr <- aggs_lrK$Flux_Total / QSWatershedArea\n",
"aggs_lcKConc <- aggregateSolute(preds_lcKConc,meta, format=\"conc\", se.preds=preds_lcKConc$se.pred, agg.by=\"calendar year\")\n",
"aggs_lrK$K_mg_l <- aggs_lcKConc$Conc\n",
"FluxTotalsdfK <- aggs_lrK[, c('K_mg_l','K_Kg_Ha_yr','Calendar_Year')]\n",
"\n",
"#round values \n",
"years <- FluxTotalsdfK$Calendar_Year\n",
"FluxTotalsdfK$Calendar_Year <- NULL\n",
"FluxTotalsdfK <- round(FluxTotalsdfK, 2)\n",
"FluxTotalsdfK$Calendar_Year <- years\n",
"\n",
"FluxTotals <- merge(FluxTotalsdfK,FluxTotals, by = \"Calendar_Year\", all.y=TRUE)\n",
"FluxTotals"
]
},
{
"cell_type": "markdown",
"metadata": {},
"source": [
"Format the table of fluxes and concentrations. "
]
},
{
"cell_type": "code",
"execution_count": 20,
"metadata": {},
"outputs": [
{
"data": {
"text/html": [
"\n",
" | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 |
\n",
"\n",
"\tCalendar_Year | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 |
\n",
"\tK_mg_l | 0.29 | 0.27 | 0.21 | 0.24 | 0.23 | 0.30 |
\n",
"\tK_Kg_Ha_yr | 6.99 | 11.45 | 10.49 | 9.55 | 8.80 | 8.55 |
\n",
"\tCl_mg_l | 7.22 | 6.69 | 6.72 | 7.02 | 7.25 | 7.38 |
\n",
"\tCl_Kg_Ha_yr | 181.60 | 306.03 | 327.98 | 281.04 | 290.83 | 213.48 |
\n",
"\tDOC_mg_l | 2.22 | 2.05 | 1.95 | 1.74 | 1.76 | 1.61 |
\n",
"\tDOC_Kg_Ha_yr | 98.32 | 159.90 | 152.43 | 97.33 | 143.95 | 62.29 |
\n",
"\tDON_mg_l | 0.07 | 0.06 | 0.07 | 0.04 | 0.06 | 0.06 |
\n",
"\tDON_Kg_Ha_yr | 1.99 | 2.93 | 4.47 | 1.91 | 2.50 | 1.91 |
\n",
"\tTDN_mg_l | 0.14 | 0.18 | 0.14 | 0.09 | 0.10 | 0.12 |
\n",
"\tTDN_Kg_Ha_yr | 3.16 | 6.38 | 6.53 | 3.27 | 3.47 | 3.36 |
\n",
"\tNH4_mg_l | 5.801 | 8.741 | 4.176 | 4.791 | 5.029 | 5.423 |
\n",
"\tNH4_Kg_Ha_yr | 0.17 | 0.43 | 0.22 | 0.17 | 0.19 | 0.19 |
\n",
"\tPO4_mg_l | 0.003 | 0.002 | 0.002 | 0.002 | 0.002 | 0.003 |
\n",
"\tPO4_Kg_Ha_yr | 0.06 | 0.11 | 0.12 | 0.09 | 0.10 | 0.09 |
\n",
"\tF_mg_l | NA | NA | NA | NA | 0.017 | 0.019 |
\n",
"\tF_Kg_Ha_yr | NA | NA | NA | NA | 0.589 | 0.452 |
\n",
"\tNO3_mg_l | 0.067 | 0.061 | 0.056 | 0.042 | 0.037 | 0.055 |
\n",
"\tNO3_Kg_Ha_yr | 1.033 | 1.798 | 1.873 | 1.193 | 0.903 | 1.115 |
\n",
"\tSO4_mg_l | 0.55 | 0.53 | 0.52 | 0.52 | 0.51 | 0.54 |
\n",
"\tSO4_Kg_Ha_yr | 15.21 | 25.03 | 25.99 | 20.92 | 20.67 | 15.46 |
\n",
"\tNO3_mg_l.y | 0.07 | 0.06 | 0.06 | 0.04 | 0.04 | 0.05 |
\n",
"\tNO3_Kg_Ha_yr.y | 1.03 | 1.80 | 1.87 | 1.19 | 0.90 | 1.11 |
\n",
"\tMg_mg_l | 1.69 | 1.39 | 1.40 | 1.46 | 1.58 | 2.39 |
\n",
"\tMg_Kg_Ha_yr | 34.81 | 59.00 | 62.26 | 54.64 | 54.93 | 38.81 |
\n",
"\tCa_mg_l | 2.70 | 2.61 | 2.83 | 2.67 | 2.52 | 2.92 |
\n",
"\tCa_Kg_Ha_yr | 56.91 | 108.81 | 129.89 | 103.36 | 93.10 | 75.59 |
\n",
"\tNa_mg_l | 4.76 | 4.50 | 4.66 | 4.64 | 4.60 | 4.91 |
\n",
"\tNa_Kg_Ha_yr | 109.37 | 193.90 | 225.30 | 179.82 | 171.38 | 134.12 |
\n",
"\n",
"
\n"
],
"text/latex": [
"\\begin{tabular}{r|llllll}\n",
" & 2009 & 2010 & 2011 & 2012 & 2013 & 2014\\\\\n",
"\\hline\n",
"\tCalendar\\_Year & 2009 & 2010 & 2011 & 2012 & 2013 & 2014 \\\\\n",
"\tK\\_mg\\_l & 0.29 & 0.27 & 0.21 & 0.24 & 0.23 & 0.30 \\\\\n",
"\tK\\_Kg\\_Ha\\_yr & 6.99 & 11.45 & 10.49 & 9.55 & 8.80 & 8.55 \\\\\n",
"\tCl\\_mg\\_l & 7.22 & 6.69 & 6.72 & 7.02 & 7.25 & 7.38 \\\\\n",
"\tCl\\_Kg\\_Ha\\_yr & 181.60 & 306.03 & 327.98 & 281.04 & 290.83 & 213.48\\\\\n",
"\tDOC\\_mg\\_l & 2.22 & 2.05 & 1.95 & 1.74 & 1.76 & 1.61 \\\\\n",
"\tDOC\\_Kg\\_Ha\\_yr & 98.32 & 159.90 & 152.43 & 97.33 & 143.95 & 62.29\\\\\n",
"\tDON\\_mg\\_l & 0.07 & 0.06 & 0.07 & 0.04 & 0.06 & 0.06 \\\\\n",
"\tDON\\_Kg\\_Ha\\_yr & 1.99 & 2.93 & 4.47 & 1.91 & 2.50 & 1.91 \\\\\n",
"\tTDN\\_mg\\_l & 0.14 & 0.18 & 0.14 & 0.09 & 0.10 & 0.12 \\\\\n",
"\tTDN\\_Kg\\_Ha\\_yr & 3.16 & 6.38 & 6.53 & 3.27 & 3.47 & 3.36 \\\\\n",
"\tNH4\\_mg\\_l & 5.801 & 8.741 & 4.176 & 4.791 & 5.029 & 5.423 \\\\\n",
"\tNH4\\_Kg\\_Ha\\_yr & 0.17 & 0.43 & 0.22 & 0.17 & 0.19 & 0.19 \\\\\n",
"\tPO4\\_mg\\_l & 0.003 & 0.002 & 0.002 & 0.002 & 0.002 & 0.003 \\\\\n",
"\tPO4\\_Kg\\_Ha\\_yr & 0.06 & 0.11 & 0.12 & 0.09 & 0.10 & 0.09 \\\\\n",
"\tF\\_mg\\_l & NA & NA & NA & NA & 0.017 & 0.019 \\\\\n",
"\tF\\_Kg\\_Ha\\_yr & NA & NA & NA & NA & 0.589 & 0.452 \\\\\n",
"\tNO3\\_mg\\_l & 0.067 & 0.061 & 0.056 & 0.042 & 0.037 & 0.055 \\\\\n",
"\tNO3\\_Kg\\_Ha\\_yr & 1.033 & 1.798 & 1.873 & 1.193 & 0.903 & 1.115 \\\\\n",
"\tSO4\\_mg\\_l & 0.55 & 0.53 & 0.52 & 0.52 & 0.51 & 0.54 \\\\\n",
"\tSO4\\_Kg\\_Ha\\_yr & 15.21 & 25.03 & 25.99 & 20.92 & 20.67 & 15.46 \\\\\n",
"\tNO3\\_mg\\_l.y & 0.07 & 0.06 & 0.06 & 0.04 & 0.04 & 0.05 \\\\\n",
"\tNO3\\_Kg\\_Ha\\_yr.y & 1.03 & 1.80 & 1.87 & 1.19 & 0.90 & 1.11 \\\\\n",
"\tMg\\_mg\\_l & 1.69 & 1.39 & 1.40 & 1.46 & 1.58 & 2.39 \\\\\n",
"\tMg\\_Kg\\_Ha\\_yr & 34.81 & 59.00 & 62.26 & 54.64 & 54.93 & 38.81 \\\\\n",
"\tCa\\_mg\\_l & 2.70 & 2.61 & 2.83 & 2.67 & 2.52 & 2.92 \\\\\n",
"\tCa\\_Kg\\_Ha\\_yr & 56.91 & 108.81 & 129.89 & 103.36 & 93.10 & 75.59\\\\\n",
"\tNa\\_mg\\_l & 4.76 & 4.50 & 4.66 & 4.64 & 4.60 & 4.91 \\\\\n",
"\tNa\\_Kg\\_Ha\\_yr & 109.37 & 193.90 & 225.30 & 179.82 & 171.38 & 134.12\\\\\n",
"\\end{tabular}\n"
],
"text/markdown": [
"\n",
"| | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 | \n",
"|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|\n",
"| Calendar_Year | 2009 | 2010 | 2011 | 2012 | 2013 | 2014 | \n",
"| K_mg_l | 0.29 | 0.27 | 0.21 | 0.24 | 0.23 | 0.30 | \n",
"| K_Kg_Ha_yr | 6.99 | 11.45 | 10.49 | 9.55 | 8.80 | 8.55 | \n",
"| Cl_mg_l | 7.22 | 6.69 | 6.72 | 7.02 | 7.25 | 7.38 | \n",
"| Cl_Kg_Ha_yr | 181.60 | 306.03 | 327.98 | 281.04 | 290.83 | 213.48 | \n",
"| DOC_mg_l | 2.22 | 2.05 | 1.95 | 1.74 | 1.76 | 1.61 | \n",
"| DOC_Kg_Ha_yr | 98.32 | 159.90 | 152.43 | 97.33 | 143.95 | 62.29 | \n",
"| DON_mg_l | 0.07 | 0.06 | 0.07 | 0.04 | 0.06 | 0.06 | \n",
"| DON_Kg_Ha_yr | 1.99 | 2.93 | 4.47 | 1.91 | 2.50 | 1.91 | \n",
"| TDN_mg_l | 0.14 | 0.18 | 0.14 | 0.09 | 0.10 | 0.12 | \n",
"| TDN_Kg_Ha_yr | 3.16 | 6.38 | 6.53 | 3.27 | 3.47 | 3.36 | \n",
"| NH4_mg_l | 5.801 | 8.741 | 4.176 | 4.791 | 5.029 | 5.423 | \n",
"| NH4_Kg_Ha_yr | 0.17 | 0.43 | 0.22 | 0.17 | 0.19 | 0.19 | \n",
"| PO4_mg_l | 0.003 | 0.002 | 0.002 | 0.002 | 0.002 | 0.003 | \n",
"| PO4_Kg_Ha_yr | 0.06 | 0.11 | 0.12 | 0.09 | 0.10 | 0.09 | \n",
"| F_mg_l | NA | NA | NA | NA | 0.017 | 0.019 | \n",
"| F_Kg_Ha_yr | NA | NA | NA | NA | 0.589 | 0.452 | \n",
"| NO3_mg_l | 0.067 | 0.061 | 0.056 | 0.042 | 0.037 | 0.055 | \n",
"| NO3_Kg_Ha_yr | 1.033 | 1.798 | 1.873 | 1.193 | 0.903 | 1.115 | \n",
"| SO4_mg_l | 0.55 | 0.53 | 0.52 | 0.52 | 0.51 | 0.54 | \n",
"| SO4_Kg_Ha_yr | 15.21 | 25.03 | 25.99 | 20.92 | 20.67 | 15.46 | \n",
"| NO3_mg_l.y | 0.07 | 0.06 | 0.06 | 0.04 | 0.04 | 0.05 | \n",
"| NO3_Kg_Ha_yr.y | 1.03 | 1.80 | 1.87 | 1.19 | 0.90 | 1.11 | \n",
"| Mg_mg_l | 1.69 | 1.39 | 1.40 | 1.46 | 1.58 | 2.39 | \n",
"| Mg_Kg_Ha_yr | 34.81 | 59.00 | 62.26 | 54.64 | 54.93 | 38.81 | \n",
"| Ca_mg_l | 2.70 | 2.61 | 2.83 | 2.67 | 2.52 | 2.92 | \n",
"| Ca_Kg_Ha_yr | 56.91 | 108.81 | 129.89 | 103.36 | 93.10 | 75.59 | \n",
"| Na_mg_l | 4.76 | 4.50 | 4.66 | 4.64 | 4.60 | 4.91 | \n",
"| Na_Kg_Ha_yr | 109.37 | 193.90 | 225.30 | 179.82 | 171.38 | 134.12 | \n",
"\n",
"\n"
],
"text/plain": [
" 2009 2010 2011 2012 2013 2014 \n",
"Calendar_Year 2009 2010 2011 2012 2013 2014 \n",
"K_mg_l 0.29 0.27 0.21 0.24 0.23 0.30 \n",
"K_Kg_Ha_yr 6.99 11.45 10.49 9.55 8.80 8.55 \n",
"Cl_mg_l 7.22 6.69 6.72 7.02 7.25 7.38 \n",
"Cl_Kg_Ha_yr 181.60 306.03 327.98 281.04 290.83 213.48\n",
"DOC_mg_l 2.22 2.05 1.95 1.74 1.76 1.61 \n",
"DOC_Kg_Ha_yr 98.32 159.90 152.43 97.33 143.95 62.29\n",
"DON_mg_l 0.07 0.06 0.07 0.04 0.06 0.06 \n",
"DON_Kg_Ha_yr 1.99 2.93 4.47 1.91 2.50 1.91 \n",
"TDN_mg_l 0.14 0.18 0.14 0.09 0.10 0.12 \n",
"TDN_Kg_Ha_yr 3.16 6.38 6.53 3.27 3.47 3.36 \n",
"NH4_mg_l 5.801 8.741 4.176 4.791 5.029 5.423 \n",
"NH4_Kg_Ha_yr 0.17 0.43 0.22 0.17 0.19 0.19 \n",
"PO4_mg_l 0.003 0.002 0.002 0.002 0.002 0.003 \n",
"PO4_Kg_Ha_yr 0.06 0.11 0.12 0.09 0.10 0.09 \n",
"F_mg_l NA NA NA NA 0.017 0.019 \n",
"F_Kg_Ha_yr NA NA NA NA 0.589 0.452 \n",
"NO3_mg_l 0.067 0.061 0.056 0.042 0.037 0.055 \n",
"NO3_Kg_Ha_yr 1.033 1.798 1.873 1.193 0.903 1.115 \n",
"SO4_mg_l 0.55 0.53 0.52 0.52 0.51 0.54 \n",
"SO4_Kg_Ha_yr 15.21 25.03 25.99 20.92 20.67 15.46 \n",
"NO3_mg_l.y 0.07 0.06 0.06 0.04 0.04 0.05 \n",
"NO3_Kg_Ha_yr.y 1.03 1.80 1.87 1.19 0.90 1.11 \n",
"Mg_mg_l 1.69 1.39 1.40 1.46 1.58 2.39 \n",
"Mg_Kg_Ha_yr 34.81 59.00 62.26 54.64 54.93 38.81 \n",
"Ca_mg_l 2.70 2.61 2.83 2.67 2.52 2.92 \n",
"Ca_Kg_Ha_yr 56.91 108.81 129.89 103.36 93.10 75.59\n",
"Na_mg_l 4.76 4.50 4.66 4.64 4.60 4.91 \n",
"Na_Kg_Ha_yr 109.37 193.90 225.30 179.82 171.38 134.12"
]
},
"metadata": {},
"output_type": "display_data"
}
],
"source": [
"rownames(FluxTotals) <- FluxTotals$Calendar_Year\n",
"FluxTotals2 <- as.data.frame(t(FluxTotals))\n",
"FluxTotals2$Calendar_Year <- NULL\n",
"FluxTotals2"
]
},
{
"cell_type": "code",
"execution_count": null,
"metadata": {},
"outputs": [],
"source": []
}
],
"metadata": {
"kernelspec": {
"display_name": "R",
"language": "R",
"name": "ir"
},
"language_info": {
"codemirror_mode": "r",
"file_extension": ".r",
"mimetype": "text/x-r-source",
"name": "R",
"pygments_lexer": "r",
"version": "3.5.1"
}
},
"nbformat": 4,
"nbformat_minor": 2
}