Contaminant Plume Dispersion Interactive Calculator

Results

What is the difference between seepage velocity and Darcy velocity, and which should I use for plume calculations? +

Darcy velocity (also called specific discharge) represents the volumetric flow rate per unit cross-sectional area of the entire porous medium, including both pore space and solid grains. Seepage velocity (also called pore velocity or average linear velocity) represents the actual average speed at which water molecules travel through the interconnected pores. Because water can only flow through void spaces, not through solid grains, seepage velocity always exceeds Darcy velocity by a factor of the effective porosity: v = q/n_e, where q is Darcy velocity and n_e is effective porosity (typically 0.15-0.35 for unconsolidated sediments). For contaminant transport calculations, you must use seepage velocity because dissolved chemicals travel with the moving water through pores. Using Darcy velocity would severely underpredict transport rates and arrival times. As a practical example, if Darcy velocity is 0.10 m/day and effective porosity is 0.25, the correct seepage velocity is 0.40 m/day—four times faster than mistakenly using Darcy velocity would suggest. This distinction frequently causes confusion among practitioners transitioning from traditional groundwater flow analysis to transport modeling.

How do I determine appropriate dispersivity values for my site when no tracer test data exists? +

In the absence of site-specific tracer testing, you must rely on empirical relationships and published values from similar geological settings, while acknowledging significant uncertainty. The most widely applied approach uses the scale-dependent relationship where longitudinal dispersivity approximately equals 10% of the characteristic travel distance (α_L ≈ 0.1×L). For a plume extending 100 meters, this suggests α_L = 10 meters. However, published compilations show actual field measurements vary by orders of magnitude at any given scale. For conservative (protective) screening assessments, use lower dispersivity values which produce less spreading and higher downgradient concentrations. For unconsolidated sand and gravel aquifers, typical ranges are α_L = 1-20 m, α_T = 0.1-2 m, and α_V = 0.01-0.2 m at 100-meter scales. Fractured rock systems exhibit completely different behavior dominated by fracture network geometry rather than porous media dispersion. Conservative ratios of α_T/α_L = 0.1 and α_V/α_L = 0.01 provide reasonable starting estimates. Sensitivity analysis evaluating results across the full plausible range quantifies uncertainty and identifies whether conclusions remain robust or require site-specific testing for defensible decision-making. For critical applications like drinking water protection or litigation support, the cost of proper aquifer characterization through tracer testing represents essential investment.

When are Gaussian plume models inappropriate, and what alternatives should I consider? +

Gaussian analytical solutions break down under several common field conditions requiring numerical modeling approaches. Time-varying sources (episodic spills, seasonal loading) violate steady-state assumptions; use transient solutions or numerical codes like MT3D or MODFLOW-GWT. Strongly heterogeneous geology with preferential pathways invalidates the homogeneous medium assumption; particle tracking through calibrated flow models captures complex migration better than analytical averaging. Reactive transport involving biodegradation, radioactive decay, or chemical precipitation requires coupled reaction-transport simulators since the Gaussian model assumes conservative (non-reactive) chemicals. Dense non-aqueous phase liquids (DNAPLs) like chlorinated solvents denser than water sink below the water table creating source zones violating point-source geometry assumptions; multiphase flow and transport codes become necessary. Sorbing contaminants exhibiting retardation require modified analytical solutions incorporating retardation factors or numerical approaches handling spatially variable partitioning. Shallow aquifers with significant interaction between groundwater and surface water need coupled models. For screening assessments, regulatory compliance demonstrations with conservative assumptions, or educational purposes, Gaussian models remain invaluable. For detailed site characterization, remediation design optimization, or forensic source reconstruction, invest in appropriate numerical modeling with qualified expert implementation. The key professional judgment involves matching model complexity to decision requirements and data quality—overly complex models applied to sparse data produce illusions of precision rather than insight.

How does biodegradation affect actual plume behavior, and can I account for it with simple corrections? +

Biodegradation can dramatically reduce contaminant concentrations through microbial transformation, particularly for petroleum hydrocarbons undergoing aerobic respiration or chlorinated solvents experiencing reductive dechlorination under anaerobic conditions. For first-order decay kinetics (most biological processes approximate this), you can modify the Gaussian solution by multiplying by an additional exponential decay term: exp(−λt), where λ is the first-order decay constant (day⁻¹) and t is travel time (x/v). A BTEX compound with λ = 0.01 day⁻¹ traveling 100 meters at 0.5 m/day experiences travel time of 200 days, yielding decay factor exp(−0.01×200) = 0.135, reducing concentrations to 13.5% of conservative predictions—a substantial difference. However, determining appropriate decay constants requires site-specific microcosm studies or calibration against monitored natural attenuation data. Biodegradation rates vary enormously with electron acceptor availability (oxygen, nitrate, sulfate, iron), nutrient concentrations, temperature, pH, and contaminant type. Petroleum hydrocarbons often degrade rapidly (half-lives weeks to months), while some chlorinated compounds persist for decades or require specialized microorganisms. Regulatory applications typically prohibit taking credit for biodegradation without extensive documentation meeting protocols like EPA's Monitored Natural Attenuation guidelines. For conservative screening assessments, assume no degradation. For detailed fate and transport analysis supporting remediation decisions, invest in proper characterization of biodegradation capacity through multiple lines of evidence: contaminant concentration trends, daughter product presence, geochemical indicators, and laboratory rate measurements. Simple corrections provide order-of-magnitude estimates but cannot replace rigorous site investigation.

What sampling strategies best validate plume dispersion model predictions in the field? +

Effective model validation requires monitoring well networks designed to capture both centerline concentrations and lateral/vertical extent boundaries. Install monitoring wells at multiple downgradient distances (near-field, mid-field, far-field) to assess concentration attenuation trends. At each distance, position wells along the predicted plume centerline plus lateral wells offset by 1σ_y and 2σ_y to characterize transverse spreading—theoretical Gaussian profiles predict concentrations at ±1σ equal to 61% of centerline values and at ±2σ equal to 14% of centerline values. Vertical profiling through multi-level sampling systems or discrete-interval sampling reveals whether plumes remain confined to predicted depths or migrate vertically through preferential pathways. Critical validation metrics include: (1) does the plume centerline align with predicted flow direction?, (2) do concentration magnitudes match predictions within uncertainty bounds?, (3) do lateral width measurements correspond to calculated standard deviations?, and (4) do arrival times at downgradient wells match travel time calculations? Significant discrepancies indicate model assumption violations: unexpected plume width suggests dispersivity estimates require revision, concentration attenuation exceeding predictions indicates biodegradation or sorption processes, asymmetric lateral spreading reveals aquifer heterogeneity, and early arrival times expose preferential flow paths. Quarterly to annual sampling over multi-year periods establishes temporal trends distinguishing transient conditions from steady-state behavior. Geochemical parameter monitoring (dissolved oxygen, redox potential, electron acceptors/donors, degradation products) provides mechanistic understanding supporting fate predictions. This integrated approach transforms monitoring from simple compliance checking into genuine understanding of subsurface processes, enabling adaptive management as site conditions evolve.

How do I handle source zones of finite size rather than idealized point sources? +

Real contamination sources rarely approximate mathematical point sources—they extend over areas like leaking tanks, disposal trenches, or infiltration ponds. For elongated sources aligned with groundwater flow (like trenches), treat the source as multiple point sources distributed along the length and superimpose their individual Gaussian solutions at downgradient locations. For sources perpendicular to flow or covering significant areas, the initial plume width depends on source geometry rather than purely dispersive spreading. A practical approach divides the source into a grid of point sources, calculates each contribution at the receptor location, and sums the results accounting for different travel distances and lateral offsets from each source element. Alternatively, for rectangular sources of width W and height H, analytical solutions exist treating them as plane sources where initial σ_y = W/(2√3) and σ_z = H/(2√3), then allowing further dispersion downgradient using standard relationships. This approach recognizes that dispersion from finite sources produces wider plumes at near-field locations than point-source assumptions predict, but differences diminish far downgradient as dispersive spreading dominates over initial geometry. For sources larger than approximately 10% of the downgradient travel distance, finite-source corrections significantly improve prediction accuracy. Source delineation through soil boring networks or geophysical surveys provides essential input defining actual source dimensions rather than assuming point geometry. When source geometry remains uncertain, conservative practice assumes the smallest source dimensions producing highest concentrations, though this may overestimate risks at near-field locations while underestimating lateral extent.