This function runs the EPGM model for a specific simulated period. The model is based primarily upon data collected in the early 1990's along the phosphorus gradient in WCA-2A. Substantial additional data collected since then in WCA-2A and other locations indicate a need to recalibrate the model and potentially revise its structure. Recent data suggest, for example, that the relationship between cattail density and soil P needs recalibration and that actual soil P thresholds for biological impacts are probably lower than reflected in the original calibrations. There are also issues relating to interpretation of and potential anomalies in the historical soil P calibration data attributed to variations in soil core collection method and definition of the soil/water interface (inclusion vs. exclusion of floc layer). There are also indications in the recent data of biologically-mediated vertical transport and/or mixing that are not reflected in the current model structure.

EPGMProfile(
  case.no = NA,
  Start.Discharge = NA,
  STA.outflow.TPconc = NA,
  STA.outflow.vol = NA,
  FlowPath.width = NA,
  Hydroperiod = NA,
  Soil.Depth = NA,
  Soil.BulkDensity.initial = NA,
  Soil.TPConc.initial = NA,
  Vertical.SoilTPGradient.initial = NA,
  Soil.BulkDensity.final = NA,
  PSettlingRate = NA,
  P.AtmoDep = NA,
  Rainfall = NA,
  ET = NA,
  Yr.Display = 30,
  Max.Yrs = 200,
  Max.Dist = 15,
  Dist.increment.km = 0.1,
  plot.profile = TRUE,
  raw.output = FALSE,
  results.table = TRUE,
  summary.distance = c(0, 0.5, 1, 2, 4, 8, 10)
)

Arguments

case.no

Case number from the pre-loaded example data (values ranges from 1 to 12)

Start.Discharge

The year which this particular STA began discharge operations.

STA.outflow.TPconc

Outflow total phosphorus concentration (in ug L-1; micrograms per liter) for this STA.

STA.outflow.vol

Annual outflow discharge volume (in x1000 Acre-Feet Year-1) for this STA.

FlowPath.width

The width of the downstream flow path (in kilometers).

Hydroperiod

Average hydroperiod (time above ground surface) of the downstream system (in percent).

Soil.Depth

Depth of soil (in centimeters).

Soil.BulkDensity.initial

The initial bulk density prior to dicharge of the soil downstream of the system (in g cm-3).

Soil.TPConc.initial

The initial total phosphorus concentration of soil prior to discharge downstream of the system (in mg kg-1).

Vertical.SoilTPGradient.initial

The soil total phosphorus concentration gradient prior to dischage downstream of the system (in mg cm-3 cm-1).

Soil.BulkDensity.final

The final bulk density after dischage of the soil downstream of the system (in g cm-3).

PSettlingRate

The phosphorus settling rate estimated from steady-state conditions (m Year-1).

P.AtmoDep

Phosphorus atmospheric depostition loading rate (in mg m-2 Year-1).

Rainfall

Annual accumulated rainfall estimate (m Year-1).

ET

Annual evapotranspiration estimate (m Year-1).

Yr.Display

Output displays results for this time (years)

Max.Yrs

Maximum number of years simulated

Max.Dist

Maximum ditance plotted, default is 50 km

Dist.increment.km

Distance increment modeled

plot.profile

If TRUE base plot will be generate with water column distance, soil distance and cattail distance profiles.

raw.output

If TRUE a data.frame will be printed with all calculations used to estimate various parameters.Default is set to FALSE.

results.table

if TRUE summary results table will be printed in the console. Default is set to TRUE.

summary.distance

Default is c(0,0.5,1,2,4,8,10) but can be changed. Values determine what distances will be included in the summary table.

Value

This function computes and plots the distance profile along the gradient based on input values

Details

As described in the original documentation, the model is designed to simulate marsh enrichment (responses to increasing P load), not recovery (responses to decreasing in load).

Examples

EPGMProfile(case.no=11)
#> Warning: conversion failure on 'Concentration (μ g L⁻¹)' in 'mbcsToSbcs': dot substituted for <e2>
#> Warning: conversion failure on 'Concentration (μ g L⁻¹)' in 'mbcsToSbcs': dot substituted for <81>
#> Warning: conversion failure on 'Concentration (μ g L⁻¹)' in 'mbcsToSbcs': dot substituted for <bb>
#> Warning: conversion failure on 'Concentration (mg kg⁻¹)' in 'mbcsToSbcs': dot substituted for <e2>
#> Warning: conversion failure on 'Concentration (mg kg⁻¹)' in 'mbcsToSbcs': dot substituted for <81>
#> Warning: conversion failure on 'Concentration (mg kg⁻¹)' in 'mbcsToSbcs': dot substituted for <bb>
#> $Time.yrs #> [1] 30 #> #> $Simulated.Zone #> Parameter Value #> 1 Distance.km 15.0 #> 2 Width.km 10.5 #> 3 Area.km2 157.5 #> 4 STA.outflow.volume.kAcftyr 281.3 #> 5 Hydroperiod.pct 91.4 #> 6 Soil.Depth.cm 10.0 #> 7 P.Settle.Rate.myr 10.2 #> 8 STA.outflow.Conc.ugL 122.0 #> 9 STA.outflow.Load.mtyr 42.4 #> #> $DistanceProfile #> 0 0.5 1 2 4 8 #> WaterCol.Pconc.ugL 122.00 106.90 93.70 72.10 43.10 16.90 #> SteadyState.WC.Conc.ugL 122.00 106.90 93.70 72.10 43.10 16.90 #> SteadyState.Soil.Conc.mgkg 2131.00 1925.00 1744.00 1449.00 1053.00 694.00 #> Time.to.Steady.State.yrs 15.00 15.50 16.00 17.20 20.90 35.20 #> NewSoil.Depth.cm 10.00 10.00 10.00 10.00 10.00 8.50 #> Soil.Mass.Accret.kgm2yr 0.67 0.65 0.63 0.58 0.48 0.28 #> Cattail.Density.pct 100.00 100.00 99.00 95.00 53.00 5.00 #> SteadyState.Cattail.Density.pct 100.00 100.00 99.00 95.00 53.00 9.00 #> 10 #> WaterCol.Pconc.ugL 11.50 #> SteadyState.WC.Conc.ugL 11.50 #> SteadyState.Soil.Conc.mgkg 620.00 #> Time.to.Steady.State.yrs 46.30 #> NewSoil.Depth.cm 6.50 #> Soil.Mass.Accret.kgm2yr 0.22 #> Cattail.Density.pct 2.00 #> SteadyState.Cattail.Density.pct 5.00 #> #> $Water.Budget #> Total.Flow.m Total.Flow.hm3 Sim.Avg.Flow.myr #> Inflow 66.14 10417 2.20 #> Rainfall 34.80 5481 1.16 #> ET 41.40 6520 1.38 #> Outflow 59.54 9378 1.98 #> #> $P.MassBalance #> PMass.mgm2 PMass.mtons Sim.Avg.Load.mgm2yr #> Inflow 8069 1270.9 269.0 #> Rainfall 1493 235.1 49.8 #> Removal 8986 1415.3 299.5 #> Outflow 576 90.7 19.2 #> #> $Soils #> SoilMass.kgm2 PMass.mgm2 PConc.mgkg BulkDensity.gcm3 PVol.mgcm3 #> Initial Storage 10.20 2020 198 0.102 0.020 #> Current Storage 8.46 6761 800 0.085 0.068 #> Accretion 8.22 8986 1093 0.080 0.087 #> Burial 9.97 4245 426 0.097 0.041 #>
EPGMProfile(NA,1991,38,526,15.3,50,10,0.05,257,-0.004,0.04,15.2,45,1.3,1.4)
#> Warning: conversion failure on 'Concentration (μ g L⁻¹)' in 'mbcsToSbcs': dot substituted for <e2>
#> Warning: conversion failure on 'Concentration (μ g L⁻¹)' in 'mbcsToSbcs': dot substituted for <81>
#> Warning: conversion failure on 'Concentration (μ g L⁻¹)' in 'mbcsToSbcs': dot substituted for <bb>
#> Warning: conversion failure on 'Concentration (mg kg⁻¹)' in 'mbcsToSbcs': dot substituted for <e2>
#> Warning: conversion failure on 'Concentration (mg kg⁻¹)' in 'mbcsToSbcs': dot substituted for <81>
#> Warning: conversion failure on 'Concentration (mg kg⁻¹)' in 'mbcsToSbcs': dot substituted for <bb>
#> $Time.yrs #> [1] 30 #> #> $Simulated.Zone #> Parameter Value #> 1 Distance.km 15.0 #> 2 Width.km 15.3 #> 3 Area.km2 229.5 #> 4 STA.outflow.volume.kAcftyr 526.0 #> 5 Hydroperiod.pct 50.0 #> 6 Soil.Depth.cm 10.0 #> 7 P.Settle.Rate.myr 15.2 #> 8 STA.outflow.Conc.ugL 38.0 #> 9 STA.outflow.Load.mtyr 24.7 #> #> $DistanceProfile #> 0 0.5 1 2 4 8 #> WaterCol.Pconc.ugL 38.00 35.30 32.80 28.50 21.70 13.70 #> SteadyState.WC.Conc.ugL 38.00 35.30 32.80 28.50 21.70 13.70 #> SteadyState.Soil.Conc.mgkg 887.00 856.00 829.00 780.00 705.00 615.00 #> Time.to.Steady.State.yrs 12.30 12.80 13.30 14.40 17.10 23.70 #> NewSoil.Depth.cm 10.00 10.00 10.00 10.00 10.00 10.00 #> Soil.Mass.Accret.kgm2yr 0.81 0.78 0.75 0.69 0.59 0.42 #> Cattail.Density.pct 26.00 23.00 19.00 15.00 9.00 5.00 #> SteadyState.Cattail.Density.pct 26.00 23.00 19.00 15.00 9.00 5.00 #> 10 #> WaterCol.Pconc.ugL 11.40 #> SteadyState.WC.Conc.ugL 11.40 #> SteadyState.Soil.Conc.mgkg 589.00 #> Time.to.Steady.State.yrs 27.30 #> NewSoil.Depth.cm 10.00 #> Soil.Mass.Accret.kgm2yr 0.37 #> Cattail.Density.pct 4.00 #> SteadyState.Cattail.Density.pct 4.00 #> #> $Water.Budget #> Total.Flow.m Total.Flow.hm3 Sim.Avg.Flow.myr #> Inflow 84.88 19480 2.83 #> Rainfall 39.00 8950 1.30 #> ET 42.00 9639 1.40 #> Outflow 81.88 18791 2.73 #> #> $P.MassBalance #> PMass.mgm2 PMass.mtons Sim.Avg.Load.mgm2yr #> Inflow 3225 740.2 107.5 #> Rainfall 1755 402.8 58.5 #> Removal 3908 896.9 130.3 #> Outflow 1073 246.2 35.8 #> #> $Soils #> SoilMass.kgm2 PMass.mgm2 PConc.mgkg BulkDensity.gcm3 PVol.mgcm3 #> Initial Storage 5.00 1285 257 0.050 0.013 #> Current Storage 4.02 2633 655 0.040 0.026 #> Accretion 5.72 3908 684 0.040 0.027 #> Burial 6.70 2560 382 0.047 0.018 #>