DATABASE C:\Program Files (x86)\USGS\Phreeqc Interactive 3.0.6-7757\database\minteq.v8.dat SELECTED_OUTPUT -file Daisy model_0.045gpm_JD_10Acres_10Fe_v4.out -user_punch true #-reset true -state false -distance false -time true -reaction true -temperature true -alkalinity true USER_PUNCH -headings Al Sb As Ba Be Cd Cl Cr Cu F Fe Pb Mn N(3) N(5) Hg pH Se Ag S(6) Tl Cyanide Zn -start 10 PUNCH TOT("Al"), TOT("Sb"), TOT("As"), TOT("Ba"), TOT("Be"), TOT("Cd"), TOT("Cl"), TOT("Cr"), TOT("Cu"), TOT("F"), TOT("Fe"), TOT("Pb"), TOT("Mn"), TOT("N(3)"), TOT("N(5)"), TOT("Hg"), -LA("H+"), TOT("Se"), TOT("Ag"), TOT("S(6)"), TOT("Tl"), TOT("Cyanide"), TOT("Zn") -end END SOLUTION 0 - leachate from ET cells temp 25 pH 7.7 pe 4 O2(g) -0.68 redox pe units mg/l density 1 Alkalinity 150 as HCO3 Al 0.067 As 0.52 Ba 0.02 Be 0.0001 B 0.31 Cd 0.001 Ca 400 Cl 96 Cr 0.013 Cu 0.05 F 1.8 Fe 0.2 Pb 0.01 Mg 110 Mn 0.085 Hg 0.0089 Ni 0.012 P 0.6 K 15 Se 0.056 Ag 0.006 Na 800 charge Sr 1.9 S(6) 2400 Tl 0.0016 V 0.091 Zn 0.028 Cyanide 0.086 N 160 Sb 0.0072 -water 300 # kg ### Assumes 30% pore space end SOLUTION 1-19 rainwater temp 25 pH 5.65 pe 4 redox pe units mg/l density 1 Ca 0.74 Cl 0.38 K 0.09 Mg 0.09 N(-3) 0.41 N(5) 1.88 Na 0.8 charge S(6) 1.98 -water 300 # kg EQUILIBRIUM_PHASES 1-19 Allow to equilibriate (precipitation, no solution). Fix pe by fixing SI of O2. Al4(OH)10SO4 0 0 Alunite 0 0 Anhydrite 0 0 Ba3(AsO4)2 0 0 Barite 0 0 CO2(g) -3.5 10 Calcite 0 0 Cr2O3 0 0 Dolomite(disordered) 0 0 Dolomite(ordered) 0 0 Ferrihydrite 0 11.37 Gibbsite 0 0 Gypsum 0 0 Hydroxylapatite 0 0 Manganite 0 0 Mirabilite 0 0 O2(g) -32 10 Pyrolusite 0 0 Pyromorphite 0 0 ZnSO4:1H2O 0 0 SURFACE 1-19 #-equilibrate with solution 19 ## Groundwater chemistry #-sites DENSITY #Hfo_w 2.21325 244 1991.28 #specify a certain weight of feryhydrite, #Hfo_s 0.05675 #no link with ferriH precipitating in eq phase Hfo_sOH Ferrihydrite equilibrium_phase 0.005 64200 Hfo_wOH Ferrihydrite equilibrium_phase 0.2 -donnan 1e-008 end SOLUTION 20 Groundwater for well MW4 temp 25 pH 8.1 pe 4 O2(g) -0.68 redox pe units mg/l density 1 Alkalinity 120 as HCO3 Al 0.03 As 0.016 Ba 0.024 Be 0.002 B 0.23 Cd 0.002 Ca 11 Cl 10 Cr 0.006 Cu 0.003 F 3.2 Fe 0.02 Pb 0.001 Mg 0.24 Mn 0.0027 Hg 0.0002 Ni 0.01 K 3.4 Se 0.002 Ag 0.006 Na 73 charge S(6) 45 Tl 0.001 Zn 0.016 Cyanide 0.01 N 0.26 Sb 0.0023 -water 300 # kg ### Assumes 30% pore space END ADVECTION -cells 181 #(1m cells) -shifts 18100 #100 shift -print_cells 19 -punch_cells 19 -warnings true SAVE solution 1-19 # New solutions to go into infiltration basin USER_GRAPH 1 -headings days As Cu Hg Cd -axis_titles "Days" "Concentrations" "" -start 10 graph_x (sim_time + 4320)/(3600*24) 20 graph_y tot("As")*75000, tot("Cu")*63500, tot("Hg")*200600, tot("Cd")*112410 -end end END MIX 20 19 0.061 20 0.939 save solution 23 end #n=1 day EQUILIBRIUM_PHASES 23 Allow to equilibriate (precipitation, no solution). Fix pe by fixing SI of O2. Al4(OH)10SO4 0 0 Alunite 0 0 Anhydrite 0 0 Barite 0 0 Calcite 0 0 CO2(g) -3.5 10 Dolomite(disordered) 0 0 Dolomite(ordered) 0 0 #Fe3(OH)8 0 0 #FeAsO4:2H2O 0 0 Ferrihydrite 0 0 Gibbsite 0 0 #Goethite 0 0 Gypsum 0 0 #H-Jarosite 0 0 #K-Jarosite 0 0 Manganite 0 0 Mirabilite 0 0 #Na-Jarosite 0 0 O2(g) -32 10 Pyrolusite 0 0 ZnSO4:1H2O 0 0 SURFACE 23 Hfo_sOH Ferrihydrite equilibrium_phase 0.005 64200 Hfo_wOH Ferrihydrite equilibrium_phase 0.2 -donnan 1e-008 Use solution 23 end #PRINT # -selected_output true # -warnings 100 #SURFACE 1 Fe oxides in 1m3 alluvium - 0.14mol/kg Fe, 10% assumed Fe0H, 244m2/g, 2.27 total sites/nm2 # -sites DENSITY # Hfo_w 2.21325 244 1991.28 #amount of Hfo in 1m3 alluvium - 10% of total Fe (taking into account RMM of additonal O and H) from whole rock and assuming alluvium density of 1900 kg/m3 # Hfo_s 0.05675