PHREEQC Mixing Fraction Calculator
Calculate the endmember contribution to a mixed water sample using a conservative tracer approach commonly applied before or alongside PHREEQC modeling.
How to Calculate Mixing Fraction in PHREEQC with Technical Confidence
If you are trying to calculate a mixing fraction in PHREEQC, the key concept is simple: a measured water sample can often be represented as a combination of two or more source waters, called endmembers. In practical hydrogeochemistry, this is used to estimate how much of source A and source B are present in the mixed sample before you run deeper equilibrium, saturation, or reaction-path calculations.
In the most common two-endmember case, the fraction from endmember A is calculated using a conservative tracer concentration. The equation is: fA = (Cm – C2) / (C1 – C2), where C1 is endmember A concentration, C2 is endmember B concentration, and Cm is the mixed sample concentration. Endmember B fraction is then fB = 1 – fA. This is exactly the calculation implemented in the tool above.
Why this matters before full PHREEQC geochemical modeling
PHREEQC can do much more than linear mixing, including mineral dissolution, ion exchange, redox transformations, and gas transfer. However, if you start with a good first-pass estimate of mixing fractions from conservative tracers, your later PHREEQC model setup is cleaner and easier to calibrate. In field workflows, analysts often estimate fractions using chloride, bromide, or stable isotopes first, then build PHREEQC simulations to test whether additional reactions are needed to explain measured chemistry.
Core assumptions behind mixing fraction calculations
- The selected tracer is conservative over the travel path and time scale.
- Concentrations for both endmembers are representative and measured in the same units.
- Only two dominant endmembers control the mixed chemistry for the selected tracer.
- Analytical uncertainty is small enough not to dominate the concentration contrast between C1 and C2.
- Sampling and laboratory quality control steps are valid.
Step-by-Step Workflow for PHREEQC Projects
- Define endmembers from hydrogeologic context. Use stratigraphy, hydraulic heads, and well construction details to identify likely sources. Examples include shallow freshwater and deeper saline water, or recharge water and formation brine.
- Select a conservative tracer. Chloride is often preferred in many aquifers due to limited mineral control compared to ions such as calcium or bicarbonate that can react significantly.
- Normalize units and quality-check data. Convert all concentrations to the same unit system before calculating fractions. Check charge balance and confirm no data transcription errors.
- Calculate first-pass fractions. Use the equation in this calculator and evaluate if fA and fB remain in physically plausible ranges.
- Run PHREEQC mixing simulation. Use the MIX keyword with estimated proportions and compare modeled output against observed concentrations of reactive and conservative species.
- Iterate using independent constraints. Compare results with isotopes, electrical conductivity, and additional tracers. Adjust endmember definitions as needed.
Regulatory and Field Benchmarks You Can Use for Plausibility Checks
While mixing fraction itself is a mass-balance concept, practitioners often compare predicted concentrations with drinking-water and hydrogeologic benchmarks. The table below lists widely used values from U.S. regulatory references that can help with interpretation when your mixed sample is intended for potable use screening.
| Parameter | Benchmark Value | Type | Interpretation Use in Mixing Studies |
|---|---|---|---|
| Nitrate (as N) | 10 mg/L | EPA MCL | Checks whether mixed water may exceed health-based regulatory limits. |
| Chloride | 250 mg/L | EPA Secondary Standard | Useful taste and corrosion threshold, often paired with conservative mixing analysis. |
| Sulfate | 250 mg/L | EPA Secondary Standard | Helps evaluate aesthetic issues and broad ionic composition consistency. |
| Total Dissolved Solids (TDS) | 500 mg/L | EPA Secondary Standard | Quick aggregate indicator for salinity shifts expected from mixing fractions. |
Values shown above correspond to U.S. EPA regulatory references commonly applied during screening and interpretation.
Typical Chloride Statistics Used in Endmember Design
A second practical check is whether your chosen endmember concentrations are realistic. The ranges below are representative values used in many groundwater and salinity studies. Site-specific measurements should always override generic ranges, but these numbers are useful as sanity checks when building conceptual models.
| Water Type | Typical Chloride Concentration | Approximate Statistic Context | Mixing Interpretation |
|---|---|---|---|
| Rain-influenced recharge water | 0.2 to 5 mg/L | Low-ionic-strength meteoric signal | Useful as dilute endmember in recharge-driven systems. |
| Fresh groundwater | 5 to 50 mg/L | Common inland background range | Often represents local aquifer baseline in non-coastal settings. |
| Brackish groundwater | 250 to 5,000 mg/L | Transitional salinity interval | Typical for partial mixing between freshwater and saline sources. |
| Seawater | Approximately 19,000 mg/L | Marine chloride-dominant chemistry | Strong saline endmember in coastal intrusion modeling. |
How this connects directly to PHREEQC input files
After estimating fractions, you can convert them into PHREEQC MIX blocks. For example, if fA = 0.261 and fB = 0.739, your first trial can set solution contributions to those proportions, then compare simulated concentrations against measurements. If chloride matches but calcium and alkalinity do not, that usually indicates reactions beyond simple mixing, such as carbonate dissolution, cation exchange, or CO2 interaction.
In other words, a successful conservative tracer fraction does not guarantee all species will fit linearly. Instead, it gives a physically grounded starting point for reaction modeling. This workflow is one of the fastest ways to separate transport-driven composition changes from true geochemical transformation.
Common mistakes and how to avoid them
- Using reactive tracers as if conservative: Sodium or bicarbonate can shift due to exchange or mineral reactions.
- Ignoring denominator sensitivity: If C1 and C2 are too similar, fraction uncertainty becomes large.
- Mixing units: mg/L and mmol/L confusion can invalidate calculations immediately.
- Overlooking temporal variability: Endmembers sampled in different seasons may not represent the same process window.
- Treating out-of-range fractions as impossible by default: Values below 0 or above 1 can indicate wrong endmembers, non-conservative behavior, or additional sources.
Interpreting uncertainty like an expert
Uncertainty in mixing fractions is controlled by two dominant factors: measurement precision and contrast between endmembers. If C1 and C2 are far apart, even moderate measurement noise in Cm leads to modest fraction uncertainty. If C1 and C2 are close, the same laboratory uncertainty can produce wide fraction ranges. That is why geochemists prefer high-contrast conservative tracers whenever possible.
The calculator estimates first-order uncertainty in fA as approximately ±u/|C1 – C2|, where u is the uncertainty entered for Cm. This gives a practical confidence interval you can carry into PHREEQC scenario testing. For higher rigor, uncertainty propagation can also include C1 and C2 analytical errors and Monte Carlo simulations.
When to move from two-endmember to multi-endmember methods
If your site includes recharge, irrigation return flow, formation brine, and wastewater influence, two-endmember equations may be oversimplified. In those cases, inverse modeling with multiple tracers or isotopes is more appropriate. Still, a two-endmember estimate remains useful as an interpretive baseline and communication tool for non-specialist stakeholders.
For coastal aquifers, for example, a simple freshwater-seawater chloride fraction gives an immediate indication of intrusion extent. Then you can use PHREEQC to evaluate calcite saturation shifts, cation exchange fronts, and pH response as salinity rises.
Authoritative References for Your Workflow
- USGS PHREEQC Version 3 Software Page
- U.S. EPA Drinking Water Regulations and Contaminants
- USGS Water Science School: Salinity and Total Dissolved Solids
Practical Conclusion
To calculate mixing fraction in PHREEQC workflows, begin with a conservative-tracer mass balance, validate the fraction against field realism, then test the result in PHREEQC with MIX and reaction modules. The strongest practice is iterative: conceptual model, tracer fraction, PHREEQC simulation, independent validation, and refinement. This sequence keeps your geochemical interpretation quantitative, defensible, and decision-ready for groundwater management, contamination studies, and salinity control projects.