Characterization of unknown contaminant sources is an essential first step for effective remediation and management of contaminated groundwater resources. Linked optimization methodologies are widely used for characterization of unknown groundwater contaminated source, in terms of flux magnitude, location and activity history in contaminated aquifers   . Transport of contaminants comprising of multiple chemically reactive species, in contaminated groundwater systems, is complex and highly non-linear process. Therefore, robust numerical simulation models incorporating chemical reactions of the contaminant species in the transport process are necessary in linked simulation-optimization methodology. The purpose of the present study is to propose and evaluate a linked simulation optimization methodology, linking a multiple species’ reactive transport simulation model to an optimization model for unknown groundwater source characterization. The methodology is evaluated for a former mine site in Queensland, Australia involving multiple reactive contaminant species undergoing complex geochemical reactions.
PHREEQE  , a reactive, chemical transport model, is a widely used computer code using the two-step procedure for solution. Tebes-Stevensa, Valocchia  developed a reactive transport simulation code (FEREACT) to investigate the chemical reaction-based transport for cobalt and EDTA (Ethylene Diamine Tetra-acetic Acid) through an iron oxide-coated sand column. Prommer, Barry  presented numerical computation results for both reactive and non-reactive chemical species transport in an aquifer site contaminated by hydrocarbons. They used the PHT3D  model, coupling the non-reactive part of MT3DMS  model with the reactive geochemical model PHREEQC-2  . Moreover, MT3DMS is widely used in various chemical reactive transport simulators such as RT3D  , SEAM3D  , BioRedox  , DART  , PHAST  and PH- WAT  . These models have been coupled with saturated media flow simulation models. On the other hand, one-dimensional movement simulators such as HP1  and RICH-PHREEQ  have been used in recent years as unsaturated porous media flow models incorporating detailed geochemistry. The numerical model HYDROGEOCHEM 5.0   for simulating transport of reactive chemical species in an aquifers, has also gained popularity in recent years     . Wissmeier and Barry  presented a software tool for groundwater flow and transport movement in combination with geochemistry in two and three-dimensional systems. Their study coupled the reaction-based transport equations used in IPHREEQC  , to the ﬁnite element method solver COMSOL  for the simulation of ﬂow and transport processes in the form of a computer code. In this study, the PHREEQC/MT3DMS-based reactive multi- component transport model PHT3D  is used to simulate the reactive transport process and it is linked to Simulated Annealing (SA) based optimization algorithm to identify the optimal source characteristics.
Initial attempts for identifying the flux magnitude, location and activity duration of unknown groundwater pollution sources utilized linear optimization model based on linear response matrix approach  and statistical pattern recognition  . Other methodologies proposed latter include: determining optimal estimates of the unknown model parameters and source characteristics using non- linear maximum likelihood estimation based inverse models  ; source identification using embedded nonlinear optimization technique    ; classical optimization based approach  ; Genetic Algorithm (GA) based optimization approach   ; Artificial Neural Network (ANN) approach; heuristic harmony for source identification  ; Simulated Annealing (SA); and Adaptive Simulated Annealing (ASA) based linked optimization algorithm for source identification      ; and differential evolution algorithm for ground- water source identification   . Optimization techniques for solving source identification problems are reviewed in Chadalavada, Datta  , Ketabchi and Ataie-Ashtiani  and Sreekanth and Datta  .
A linked simulation optimization approach is implemented in the present stu- dy. This approach links multiple species reactive transport simulation model with SA based optimum decision model. Performance of its efficiency and reliability is evaluated for a contaminated aquifer site hydro-geologically representing an abandoned mine site in Queensland, Australia.
In this study, chemically reactive transport process is simulated using PHT3D, and is compared with non-reactive transport process in complex study area, such as mine site area. The proposed simulation model is then utilized in the linked simulation-optimization based methodology for simulating the transport process of the contaminant species in the aquifer for unknown groundwater source characterization in the contaminated aquifer. The following two main steps explain the proposed methodology development and its application to a mine site in Queensland, Australia. The first step comprises of implementing the flow and transport simulation model based on the estimated hydrological and geochemical properties of the mine site area. MODFLOW is used as a flow simulation model. PHT3D as a numerical coupled physical and chemical transport process simulator is utilized to simulate the reactive chemical transport processes in the contaminated aquifer. The flow simulation results are compared with field measurement data  to obtain acceptable calibration results based on measured flow data. In the second step, the simulation model is linked to Simulated Annealing (SA) optimization algorithm within a linked simulation-op- timization framework to obtain the optimal source characteristics in terms of location, flux magnitude and time of activity.
2.1. Simulation of Groundwater Flow and Physical Transport
The widely used three-dimensional numerical flow simulation model MODFLOW  is used for simulation of the flow in the selected aquifer study area. The MODFLOW can simulate steady and transient flow. As several variables affect the groundwater flow, it is generally modelled by the spatial and temporal governing partial differential equations.
MT3DMS  is a modular three-dimensional transport model used to simulate the transport process. It is capable of incorporating dispersion, advection, and chemical reactions of dissolved chemical species in groundwater systems. MT3DMS utilizes the flow simulation MODFLOW to compute the pollutants transport processes. In this study, MT3DMS is used for physical transport pro- cess. Equation (1)  represents the three-dimensional reactive transport pro- cess in groundwater systems:
C is the pollutant concentration (M∙L−3);
t is time (T);
xi is the Cartesian coordinate axis i (L);
Dij is the hydrodynamic dispersion coefficient tensor (L2∙T−1);
vi is the linear pore water or seepage velocity (L∙T−1);
qs is volumetric flux of water per unit volume of aquifer (positive for source and negative for sinks (T−1));
Cs is the source or sinks concentration (M∙L−3);
θ is the porosity of the porous medium (dimension less);
Rk is the chemical reaction term representing each of the N species (M∙L−3∙T−1).
The multiple component multiple species reactive transport model PHT3D  is utilized in this study or linking with the optimal source characterization model. PHT3D combines MT3DMS and PHREEQC  . The PHREEQC component allows for a variety of low temperature aqueous geochemical reactions. The PHT3D is used in this study to simulate both physical and chemical transport processes in the contaminant aquifer. In PHT3D several species can be integrated. The number of species and their initial concentrations on basic transport package need to be changed compared the MT3DMS module. It is also required to add the concentrations for each new species in the source/sink package. The chemical equations  which are solved by PHREEQC are:
rreac is a source/sink chemical reaction rate term; c is the concentration (molar) of the uncomplexed aqueous component, si is the concentration (molar) of the ith complexed species, n is the number of dissolved chemical species t complexed with the aqueous component, is the stoichiometric coefficient of the aqueous component in the ith complexed species.
2.2. Source Characterization Model Using Simulated Annealing Algorithm
The evolutionary optimization algorithm Simulated Annealing (SA)  is used as the optimization algorithm for optimal solution of the contaminant source characterization model linked with PHT3D as the reactive transport simulation model. The SA is based on analogy with thermodynamics. The analogy is with the process of slowly cooling (annealing) molted metals to a state of low energy. The important features of SA which make it a potentially powerful optimization algorithm are: it can handle multiple optima and the algorithm can move uphill, thus escaping from local optima to find global optima. The capability of SA for convergence to a global optimal solution in complex, multivariate problems involving higher degree non-linear functions makes it an ideal choice for solving unknown groundwater pollution source characterization problems.
In the search process of the SA algorithm, a random nearby solution is replaced utilizing the probability depending on the difference between the corresponding function values and algorithm control parameters (initial temperature, temperature reduction rate, etc.). Here, SIMANN a FORTRAN code developed by  is utilized as the optimization algorithm.
A new linked simulation optimization model linking SIMANN, the SA based optimization algorithm with reactive transport simulation model (PHT3D) is developed for this study. Figure 1 shows the schematic linking of the reactive transport simulation model (PHT3D) with the SA based optimization model in a linked simulation-optimization framework. The reactive transport simulation model PHT3D utilizes the flow field simulated by MODFLOW.
In order to achieve computational efficiency, the SA algorithm parameters need to be suitably specified. The SA parameters needed to be specified include the following:
The initial temperature (T);
The error tolerance for termination (called EPS);
The temperature reduction factor (RT).
The objective function of the optimal source characterization model can be described as   :
Figure 1. Representative flowchart of the linked simulation-optimization model using si- mulated annealing algorithm.
Subject to the implicit constraints defining the simulation model of the flow and reactive transport processes  ,
Abs is the absolute difference;
is the estimated concentration simulated by the transport model at monitoring point iob and at end of the time period k (mole/L);
nk is the number of time steps for monitoring;
nob is the total number of monitoring points;
is the observed concentration at monitoring location iob and at the end of the period k (mole/L).
3. Performance Evaluation
Performance of the proposed source characterization methodology is evaluated for a polluted aquifer study area hydrogeologically resembling the Mount Morgan, a former mine site in Queensland, Australia. The hydrogeological and geochemical parameters for the study area is obtained from previous studies [Esfahani and Datta 22; Wels and Findlater 45].
3.1. Polluted Aquifer Site Description
Mount Morgan deposit is situated in the Calliope Block, which is an important tectonic zone, and occurs along the eastern margin of Australia  . The geology of Mount Morgan mainly constitute of volcanic and metamorphic rocks (quartz, feldspar, basalt etc.). Mount Morgan mine site was operational for a period close to 100 years starting 1889 till 1990. It was once a large mine located in Queensland, Australia (Figure 2). It is surrounded by Dee River on the downstream of the mine site which flows into the Fitzroy River. The ground topography generally slopes from north east towards the river in south. Gold, silver, and copper ores were primarily mined during its operational life time.
The climate at the site show seasonal variation with average maximum monthly temperatures ranging from 32˚C in January to 23.1˚C in July. The average annual rainfall is approximately 815.5 mm of which a large portion is received during the wet season from November until March. During the winter season the average rainfall is in between 22 and 43 mm, while, from October to March, the rainfall may exceed 100 mm  .
Hazardous contaminants which mainly remain buried deep in the earth’s crust are exposed to the atmosphere due to open cuts, waste deposits from the ore and the created tailing dams. These hazardous compounds travel through the drainage path and contaminate the groundwater as well as surface water in the vicinity of the mine site. Once these hazardous contaminants and their components interact with water (in from of rainfall or get collected in the tailings dam as a result of overland flow), generate Acid Rock Drainage (ARD) and Acid Mine Drainage (AMD). Although several remediation measures were undertaken for this area, accurate contaminant sources characterization of ARD and AMD, physical and chemical behaviors of the contaminants in the groundwater system all along the transport pathway are required to achieve effective and efficient remediation strategies, and for optimal groundwater management.
During the early 1900s there was little or no environmental regulation. Heigh- tened mining activities with no environmental regulations have left a legacy of impacts on the Dee River and the site itself. The current environmental pollution
Figure 2. Plan view of the contaminated aquifer area in Queensland, Australia.
can possibly be attributed to acid and metalliferous drainage also known as acid mine drainage. Such pollution occurs when sulfidic rocks, such as pyrite, are exposed to atmospheric water and oxygen rather than being isolated as crystalline rock underground  . This exposure leads to the formation of sulfuric acid, which in turn dissolves high concentrations of salts and hazardous materials, such as copper, arsenic, nickel, cadmium, zinc, aluminum, iron and many more. Acid mine drainage is a very toxic concoction of hazardous chemicals threatening aquatic ecosystems and biodiversity  . Starting in 1982, tailings were dredged from Sandstone Gully and treated using the Carbon-In-Pulp (CIP) process for recovery of gold, before being backfilled into the open cut. After final closure in 1990, the partially backfilled open cut (and Sandstone Gully) were allowed to flood by natural inflows (surface runoff and groundwater inflow)  .
A major part of the waste rock deposit from the open cut is acid-forming as it contains up to 10% sulphur with sulphide minerals such as pyrite (FeS2, also called fool's gold). The term “Mundic” describes a copper ore that begins to melt. Historically, the Mundic tailings were placed into Mundic Creek drainage channel and the other tailings were placed into tailing dams. Geochemical tests carried out  show that the Mundic Red tailings were unreactive whereas the Mundic Grey tailings are highly reactive releasing sulphate (), iron (Fe), aluminum (Al) and copper (Cu). If tailings were deposited without proper containment, this would be a real environmental issue   . Thus, it is evident that the mining activity at Mount Morgan throughout its period of operation produced many chemically reactive hazardous pollutants in large quantities leading to environmental pollution issues.
The major ore at the site is Devonian Mount Morgan Tonalite which is almost impermeable. Devonian rhyolitic volcanic depositshave provided Au-Cu mineralization. Permeability of the porous medium comprise of primary permeability resulting from the formation of rocks (mostly considered negligible in rock formations in Australia) andsecondary permeability controlled by the presence of faults and fractures in the rock strata. Compartmented dykes in the area causes groundwater discharge deeper into the mine site  . Mining activities for extracting Pyrite ore mineral has resulted in large waste deposit comprising mainly of quartz?pyrite waste rocks. Sulphide oxidation, pyrite waste rock exposition, and also tailings at the site generate hazardous acidic and metal-rich drainage. The absence of carbonate rocks to naturalize groundwater pH has exacerbated this situation  . The geology of the site has been described in detail in    .
The high pH values as associated with AMD result in dissolution of metals (e.g., copper, iron, aluminum) in the water. As and when the pH value rises and become closer to neutrality, the dissolved metals precipitates leaving behind a yellow-orange substance or a light-blue substance as in Figure 3.
The groundwater flow occurs mainly in the permeable mine waste dumps and in shallow bedrocks that may have been fractured by mining and have become permeable. The area consists of 4 layers  : waste rock dumps and tailing; highly
Figure 3. Acid mine drainage consequences in the dee river.
weathered bedrock; partially weathered bedrock; and tight bedrock. Groundwater flow is simulated with specified hydrogeological parameters as given in Table 1. An aerial view of the mine site is shown in Figure 4.
The flow model is calibrated against the available head data. The flow models were calibrated with the aim to assign proper head boundary conditions, and estimate the spatially and temporally averaged recharge. The calibration was partially based on the earlier flow model calibration for the same study area as reported in  . The flow model was calibrated against limited and sparse measured hydraulic head field data available for the site. Hydraulic head measurement data recorded between 2002 and 2003  from 11 observation locations spread across the impacted area were used for model calibration. A portion of the average rainfall intensity per year is specified as recharge to the groundwater aquifer for calibration of the flow simulation model. Calibration targets were set to be within two meter intervals of the observed hydraulic head value with a confidence level of 90 percent. The model boundary conditions were manually adjusted based on best judgment of the existing boundary conditions, to achieve the calibration targets. Figure 5 shows the calibration results in form of box plots. More rigorous calibration can be undertaken if a large number of spatiotemporal field measurement data is available.
This study only serves as a demonstration, which includes the implementation of the flow and multispecies reactive transport simulation model in a hydrogeologic contaminated site with multiple sources of contamination. Copper and sulfate are considered as two representative chemically reactive species, although more species and more complex chemical activities are present in this site. This illustrative application is used to demonstrate the limited application of the linked simulation optimization based methodology for unknown source charac-
Figure 4. Open cut and dee river illustration at mount morgan site.
Figure 5. Flow simulation and calibration results  .
Table 1. Hydrogeologic parameters used in flow and transport models of the study area  .
terization in such a complex aquifer system, while incorporating chemically re- active transport processes. Table 2 shows the chemical equations related to the contaminant transport processes in the contaminated aquifer.
3.2. Pollutant Transport Simulation in the Mine Site Area
To demonstrate the role of chemically reactive transport process with chemical reactive multi-species in the contaminated aquifer, the reactive contaminant transport simulation results are compared with the simulation model results incorporating only the physical and chemical transport processes (MT3DMS). PHT3D is utilized to simulate the chemically reactive transport process in the contaminated aquifer. For the purpose of implementation, and to synthetically generate simulated concentration measurement data at specified monitoring locations, specified pollutant concentrations are specified along the boundary of the mining pit. The transport model uses the flow field generated by the flow simulation model to predict the movement of the pollutants in the impacted area of the aquifer over a specified five year period.
The Dispersion Package in MT3DMS and PHT3Dis utilized to simulate the effect of advection, longitudinal dispersion, transverse dispersion and the effective molecular diffusion. The relevant parameter values are given in Table 1. Advection and dispersion occur in every transport process for both non-reactive contaminants and chemical reactive contaminants. But for transport of the reactive species, chemical reactions also impact the spreading and movement of the pollutant plume. These reactions are handled by using PHT3D reactive species transport simulation model.
Since advection and dispersion are not impacted by reactions, the input parameters for these two aspects of the transport are the same as in advective and dispersive transport process. PHT3D handle the reactions that the user need to consider. In the developed linked simulation optimization model developed incorporating PHT3D, the chemical reactions need to be specified in a database file, which will be used by the model during the simulation process.
There are certain norms for creation of the reactions input data base for PHT3D. This reaction database enables consideration of the specified reactions.
Table 2. Copper reactions and constants of reaction  .
In this study, reaction databases were utilized to simulate copper related contamination. Once the database is created, it has to be linked to PHT3D in order to launch the simulation. In the Basic Transport Package for PHT3D, defining species also implies choosing a reactions database. This is where the database can be linked to the simulation program. The different active species existing in the databases is then considered in the simulation (Table 2).
Since every species (even the products of the reaction) are considered, the initial concentration must be defined for every species in all the areas and in the pit (a potential source of contamination). As only copper is considered in this study, the species that are considered for the reactions are Cu2+, and water for the reagents. Water does not have to be defined as a species; indeed it is part of the flow model. Then the products of the reactions between copper and water and copper and sulphate are: Cu(OH)+ (named Cu(3) in the database), Cu(OH)2 (named Cu(4)), (named Cu(5)), (named Cu(6)) and CuSO4 (named Cu(7)). There is also Cu+, which is called Cu(1) in the reaction database, which is formed due to the presence of electrons.
Initial value for these species in the study area is defined in Figure 6. A concentration of 1.0 × 10−30 Mole/L is set in order to avoid numerical error because of a zero concentration for all the compounds. pH and pe also have to be defined. An arbitrary value of 6 is set for the pH because the media is rather acidic than neutral with heterogeneous values. The pe is set at 13.5. The Source/Sink Mixing Package also need to be change.
3.3. Comparison of Advection Dispersion and Reaction
The concentrations of all the reaction products in the pit are the same as in all the rest of the study area. Indeed, only copper (II) is assumed to originate from the pit for this simulation. When these parameters are defined, the simulation is
Figure 6. Screen shot of the Source/Sink Mixing Package concentrations input for the reaction part.
ready to be launched. The plume contours in the study area are obtained based on advection, dispersion, and chemical reaction in the transport process. To show the contaminant transport contours with better resolution, adjust the area close to the pit it showed in Figures 7-9. Figure 7 shows the contaminants spreading with time from North to South. Indeed, with advection, the pollutants spread following the direction of flow, thus from North to South.
Figure 7. Evolution of the contaminant plume from the Open Pit between T1 = 36.5 days and T2 = 1825 days with only advection.
Figure 8. Evolution of the contaminant plume from the Open Pit between T1 = 36.5 days and T2 = 1825 days with advection and dispersion.
Figure 9. Evolution of the contaminant plume from the Open Pit between T1 = 36.5 days and T2 = 1825 days with advection, dispersion and chemical reaction.
fect of dispersion as well. The spread of the pollutant plume is similar to that observed with advection, meaning that the spread occurs generally from North to South, following flow direction. However, with the added effect of dispersion the plume is more diffuse compared to that with only advection. Indeed, with dispersion, the contaminant spreading follows the general direction of flow as well as the longitudinal and transverse directions. Thus the dispersivity allows the contaminant to reach locations in the site further than that with advection.
Subsequently, the reaction process is added to the simulation process, in order to see the effect of reactive chemical species transport with Cu2+. Figure 9 show the simulation results at the beginning and at the end of the simulation time span. Figure 9 shows that the contaminant spread is smaller compared to that for advection and dispersion without reaction. The contaminant concerned in this study is the Cu2+. However, Cu2+ chemically reacts with water and sulphate in the site. Thus it is because copper (Cu2+) is transformed to other forms of copper or products, and there is an associated loss of mass resulting in smaller concentrations when the reactive transport process is incorporated.
To illustrate a better comparison of different transport processes (Advection, Dispersion, and Chemical Reaction), the Cu2+ concentration in particular at different times is considered using each variation in the transport process. For this purpose, the point A (layer (L) 2, row (R) 53 and column (C) 30) is selected. This point is chosen because it is situated east of the open pit, where the contaminant spread is evident and the distinction between the advection, dispersion and reaction is prominent (Figure 10).
As shows in Figure 10, for the same point, the Cu2+ concentration is larger with advection and dispersion and smaller with copper reactions. Indeed, the advection allows the contaminants to be spread following the flow direction in the site. However, due to dispersion in addition to advection, at certain locations and time, larger total quantity of contaminant species may be transported, compared to advection alone. Thus, at the point A, the dispersion advection process
Figure 10. Graph representing the copper concentration (Mole/L) as a function of time at the point A, with advection, dispersion or reaction.
results in more contaminant quantity than the advection process alone. Moreover, the addition of reaction package to the dispersion transport model incorporates the copper product formation following the reactions presented in Table 3. A very small amount of Cu is produced from the chemical reaction equations. There are some reasons for it such as all equations are equilibrium and without any changes in the ion balance the Cu production is very low. Only the Cu+2 transport is shown here.
3.4. Copper Products Repartition
The different forms of copper are studied at the location A (L: 2, R: 53, C: 30). This point is located on the right hand side (east) of the open pit, in the area where the copper is prominently present and most spread out of the source. Therefore at this location, the Cu2+ reactive transport is likely. Figure 11 and Figure 12 show the different copper forms as a function of time. As shows in Figure 11, the majority of copper at the point A is in the form of Cu+ at a concentration of 2.34E−6 Mole/L at 1825 days. The reaction allowing the formation
Figure 11. Graph representing the concentrations of the different forms of copper (Mole/L) as a function of time at the point A.
Table 3. Results of the SA algorithm and comparison with real concentrations in the open pit.
Figure 12. Graph representing the concentrations of the different forms of copper (Mole/L) as a function of time at the Point A without Cu+.
of Cu+ from Cu2+ (Table 2) has a reaction constant of 2.72 that means a strong tendency to form Cu+ at the detriment of the Cu2+. The Cu+ is presented with the most concentration on this point but other products are presented in smaller concentrations (Figure 12).
As shown on Figure 12, two other species are major forms of copper at the point A are: Cu(OH)2 and the CuSO4 at concentrations of 6.87E−15 Mole/L and 1.22E−15 Mole/L respectively at 1825 days. The formation constant for the CuSO4 is 2.36, that is higher than that for Cu(OH)2 which is −16.19. But the water is in excess that could explain why the Cu(OH)2 concentration is more important than CuSO4 concentration. The other copper forms such as CuOH+, Cu(OH)3− and Cu(OH)42− are not in large enough concentration to be considered as major copper forms. Indeed, their formation constants are lower than the Cu+, Cu(OH)2 and CuSO4. That is why they are almost not traceable. Moreover, Cu(OH)2 is the product of the reaction between water and Cu2+ the most stable. Thus if formed, CuOH+ is transformed to Cu(OH)2. Very small amount if any, of Cu(OH)3− and Cu(OH)42− are formed. The four major forms of copper in the site are Cu2+, Cu+, Cu(OH)2 and CuSO4. These four species are considered as the relevant contaminants for the groundwater pollution source identification using the linked SA algorithm and PHT3D reactive transport simulation model
4. Source Characterization Using PHT3D as the Reactive Chemical Species Transport Model
PHT3D model is used as simulators of the flow and transport processes in the linked simulation-optimization model for optimal source characterization of distributed sources in the contaminated mine site area. For evaluation purpose, concentration measurements are simulated using the calibrated numerical flow and transport simulation model to generate synthetic concentration measurement data at specified monitoring locations. These concentration measurement data are generated for the study area for specified source characteristics. Then these synthetic concentration measurement data are utilized together with the linked simulation-optimization models, to evaluate the potential applicability, accuracy and feasibility of the developed methodology. The Cp it is considered as a distributed source of reactive contaminants. For simplification, it is assumed that the unknown contamination source characterization problem can be repre- sented by using unknown concentrations on the periphery of the open pit. Therefore discharge into the aquifer from the pit with specified water surface elevation (estimated by the calibrated flow model) and the unknown concentrations to be estimated as the source concentration by using the linked optimization simulation model, will determine the unknown contaminant mass flux. Therefore, the liked source characterization model essentially needs to estimate the unknown concentrations along the periphery of the open pit in order to characterize the unknown contaminant sources.
5. Results and Discussion
The performance evaluation results of source characterization using PHT3D model linked SA optimization algorithm are presented in Table 3. The contaminant source concentration magnitude is in mole per liter. In Table 3, the linked optimization model solution results for reactive transport process are compared with synthetically simulated actual concentrations at the observation wells for specified pollution sources. The errors associated with source characterization using the PHT3D model linked optimization model are very small. However, source magnitude results is comparatively large for some Cu components due to site complexity with associated chemically reactive multiple species. These performance evaluation results demonstrate that the developed methodology with PHT3D models is capable of optimally identifying the actual source magnitudes and locations in the mine site area with small estimation errors. Figure 13 illustrates the number of iteration in source characterization optimization model and the corresponding objective function value. The objective function value should converge to nearly a zero value for an optimal solution. The objective function value trend line, by decreasing, indicates that the optimal function value converges towards zero that should mean that the estimated sulphate and copper concentrations converge towards their respective actual values related to the open pit source.
In order to verify if the optimal concentration values given by SA algorithm correspond to the real initial concentrations integrated in the Open Pit, these values are used as input to the PHT3D simulation model to simulate the resulting concentrations at the monitoring locations. A typical solution result obtained for the Copper (II) at the point A (L = 2, R = 57, C = 28) is presented on Figure 14.
Figure 14 shows that the Cu2+ concentrations obtained for different times at location A (L = 2, R = 57, C = 28) are very close to the concentrations measured at this location. Indeed, the difference between the actual concentrations and the
Figure 13. Graph representing the objective function values as a function of the number of iterations for one typical trial.
Figure 14. Graphs representing the concentrations of Cu2+ obtained from the con- trol source and the concentrations of Cu2+ obtained when the source concentra- tions correspond to the concentrations given by the SA simulation, as a function of times in the point 2.
concentrations obtained by using the estimated source does not exceed 16%. Thus, this trial shows that the linked simulation optimization model incorporating reactive chemical species transport is potentially useful for estimating the unknown contaminant source accurately. No doubt more rigorous performance evaluations are necessary to establish this beyond doubt, or to determine the limitations.
Mount Morgan, a former gold, copper and silver mine in Queensland, Australia, has an impact on groundwater surface water quality. Indeed, because of the mixture of pyrite, water and oxygen in the site, Acid Rock Drainage is formed and, due to the hydrogeology of the site, is transported by groundwater through the aquifer to the Dee River, situated downstream of the mine site. This phenomenon is characterized by acid water rich in metals which polluted the ecosystem of the water systems present in the mine site.
Following this problem of contamination, the aim of this study was to establish a transport model of the contaminants: the copper (II); and to validate this model in order to be able to use it for the identification of unknown pollution sources present in the site. For that, a three-dimensional reactive transport simulation model PHT3D, along with the groundwater flow simulation model MODFLOW (also available within the processing software GMS) has been used. This reactive species transport simulation model allows numerical three-dimen- sional transient simulation of reactive chemical species.
Once the transport models for different transport process scenarios are implemented, one aim was to verify the validity of these models. For that, an unknown pollution sources characterization software is used. The dispersion and reactions models, which are the more representative of the phenomenon present on the site, are verified by testing the simulation model for a specified source: the Open Pit, and considering copper related reactions. Very limited evaluations show that the methodology can identify the source concentrations with the largest error of 16%. These errors are comparatively larger than those obtained when using MT3DMS with advection and dispersion only. However, those results do not represent realistic description of all the hydrogeochemical processes in this site. Including the reactive processes represents more realistic description of the actual transport processes in the site and these results are within acceptable range of accuracy. Therefore, the methodology using PHT3D as a simulation model linked to the optimization algorithm (SA) can be considered as potentially applicable for characterization of unknown contaminant sources in a complex contaminated aquifer with chemically reactive contaminant species. However, more rigorous evaluations are necessary to establish widespread applicability of the proposed computational code and the methodology. It is also noted that the SA optimization algorithm may have few inherent limitations in terms of efficiency in reaching an optimal solution. Application of Adoptive Simulated Annealing Algorithm (ASA)  could be a better choice to ensure more efficient convergence to a global optimal solution.
This study is a continuation of studies previously undertaken, with financial su- pport from CRC-CARE, Australia through Project No. 126.96.36.199.09/10 (2.6.03), CRC- CARE-Bithin Datta, at James Cook University, Australia.