Geochemical modeling

Last updated

Geochemical modeling or theoretical geochemistry is the practice of using chemical thermodynamics, chemical kinetics, or both, to analyze the chemical reactions that affect geologic systems, commonly with the aid of a computer. It is used in high-temperature geochemistry to simulate reactions occurring deep in the Earth's interior, in magma, for instance, or to model low-temperature reactions in aqueous solutions near the Earth's surface, the subject of this article.

Contents

Applications to aqueous systems

Geochemical modeling is used in a variety of fields, including environmental protection and remediation, [1] the petroleum industry, and economic geology. [2] Models can be constructed, for example, to understand the composition of natural waters; the mobility and breakdown of contaminants in flowing groundwater or surface water; the ion speciation of plant nutrients in soil and of regulated metals in stored solid wastes; the formation and dissolution of rocks and minerals in geologic formations in response to injection of industrial wastes, steam, or carbon dioxide; the dissolution of carbon dioxide in seawater and its effect on ocean acidification; and the generation of acidic waters and leaching of metals from mine wastes.

Development of geochemical modeling

Garrels and Thompson (1962) first applied chemical modeling to geochemistry in 25 °C and one atmosphere total pressure. Their calculation, computed by hand, is now known as an equilibrium model, which predicts species distributions, mineral saturation states, and gas fugacities from measurements of bulk solution composition. By removing small aliquots of solvent water from an equilibrated spring water and repeatedly recalculating the species distribution, Garrels and Mackenzie (1967) simulated the reactions that occur as spring water evaporated. [3] This coupling of mass transfer with an equilibrium model, known as a reaction path model, enabled geochemists to simulate reaction processes.

Helgeson (1968) introduced the first computer program to solve equilibrium and reaction path models, [4] which he and coworkers used to model geological processes like weathering, sediment diagenesis, evaporation, hydrothermal alteration, and ore deposition. [5] Later developments in geochemical modeling included reformulating the governing equations, first as ordinary differential equations, then later as algebraic equations. Additionally, chemical components came to be represented in models by aqueous species, minerals, and gases, rather than by the elements and electrons which make up the species, simplifying the governing equations and their numerical solution. [2]

Recent improvements in the power of personal computers and modeling software have made geochemical models more accessible and more flexible in their implementation. [6] Geochemists are now able to construct on their laptops complex reaction path or reactive transport models which previously would have required a supercomputer. [7]

Setting up a geochemical model

An aqueous system is uniquely defined by its chemical composition, temperature, and pressure. [8] Creating geochemical models of such systems begins by choosing the basis, the set of aqueous species, minerals, and gases which are used to write chemical reactions and express composition. The number of basis entries required equals the number of components in the system, which is fixed by the phase rule of thermodynamics. Typically, the basis is composed of water, each mineral in equilibrium with the system, each gas at known fugacity, and important aqueous species. Once the basis is defined, a modeler can solve for the equilibrium state, which is described by mass action and mass balance equations for each component. [2]

In finding the equilibrium state, a geochemical modeler solves for the distribution of mass of all species, minerals, and gases which can be formed from the basis. This includes the activity, activity coefficient, and concentration of aqueous species, the saturation state of minerals, and the fugacity of gases. Minerals with a saturation index (log Q/K) equal to zero are said to be in equilibrium with the fluid. Those with positive saturation indices are termed supersaturated, indicating they are favored to precipitate from solution. A mineral is undersaturated if its saturation index is negative, indicating that it is favored to dissolve. [8]

Geochemical modelers commonly create reaction path models to understand how systems respond to changes in composition, temperature, or pressure. By configuring the manner in which mass and heat transfer are specified (i.e., open or closed systems), models can be used to represent a variety of geochemical processes. Reaction paths can assume chemical equilibrium, or they can incorporate kinetic rate laws to calculate the timing of reactions. In order to predict the distribution in space and time of the chemical reactions that occur along a flowpath, geochemical models are increasingly being coupled with hydrologic models of mass and heat transport to form reactive transport models. [2] Specialized geochemical modeling programs that are designed as cross-linkable re-entrant software objects enable construction of reactive transport models of any flow configuration. [9]

Types of reactions

Geochemical models are capable of simulating many different types of reactions. Included among them are:

Simple phase diagrams or plots are commonly used to illustrate such geochemical reactions. Eh-pH (Pourbaix) diagrams, for example, are a special type of activity diagram which represent acid-base and redox chemistry graphically.

Uncertainties in geochemical modelling

Various sources can contribute to a range of simulation results. The range of the simulation results is defined as model uncertainty. One of the most important sources not possible to quantify is the conceptual model, which is developed and defined by the modeller. Further sources are the parameterization of the model regarding the hydraulic (only when simulating transport) and mineralogical properties. [10] The parameters used for the geochemical simulations can also contribute to model uncertainty. These are the applied thermodynamic database and the parameters for the kinetic minerals dissolution. [11] Differences in the thermodynamic data (i.e. equilibrium constants, parameters for temperature correction, activity equations and coefficients) can result in large uncertainties. Furthermore, the large spans of experimentally derived rate constants for minerals dissolution rate laws can cause large variations in simulation results. Despite this is well-known, uncertainties are not frequently considered when conducting geochemical modelling. [12]

Reducing uncertainties can be achieved by comparison of simulation results with experimental data, although experimental data does not exist at every temperature-pressure condition and for every chemical system. [12] Although such a comparison or calibration can not be conducted consequently the geochemical codes and thermodynamic databases are state-of-the-art and the most useful tools for predicting geochemical processes.

Software programs in common use

The USGS website provides free access to many of the software listed above. [35]

See also

Further reading

Related Research Articles

In physics, statistical mechanics is a mathematical framework that applies statistical methods and probability theory to large assemblies of microscopic entities. Sometimes called statistical physics or statistical thermodynamics, its applications include many problems in the fields of physics, biology, chemistry, neuroscience, computer science, information theory and sociology. Its main purpose is to clarify the properties of matter in aggregate, in terms of physical laws governing atomic motion.

Geochemistry is the science that uses the tools and principles of chemistry to explain the mechanisms behind major geological systems such as the Earth's crust and its oceans. The realm of geochemistry extends beyond the Earth, encompassing the entire Solar System, and has made important contributions to the understanding of a number of processes including mantle convection, the formation of planets and the origins of granite and basalt. It is an integrated field of chemistry and geology.

<span class="mw-page-title-main">Solvation</span> Association of molecules of a solvent with molecules or ions of a solute

Solvation describes the interaction of a solvent with dissolved molecules. Both ionized and uncharged molecules interact strongly with a solvent, and the strength and nature of this interaction influence many properties of the solute, including solubility, reactivity, and color, as well as influencing the properties of the solvent such as its viscosity and density. If the attractive forces between the solvent and solute particles are greater than the attractive forces holding the solute particles together, the solvent particles pull the solute particles apart and surround them. The surrounded solute particles then move away from the solid solute and out into the solution. Ions are surrounded by a concentric shell of solvent. Solvation is the process of reorganizing solvent and solute molecules into solvation complexes and involves bond formation, hydrogen bonding, and van der Waals forces. Solvation of a solute by water is called hydration.

<span class="mw-page-title-main">Solubility</span> Capacity of a substance to dissolve in a homogeneous way

In chemistry, solubility is the ability of a substance, the solute, to form a solution with another substance, the solvent. Insolubility is the opposite property, the inability of the solute to form such a solution.

Solubility equilibrium is a type of dynamic equilibrium that exists when a chemical compound in the solid state is in chemical equilibrium with a solution of that compound. The solid may dissolve unchanged, with dissociation, or with chemical reaction with another constituent of the solution, such as acid or alkali. Each solubility equilibrium is characterized by a temperature-dependent solubility product which functions like an equilibrium constant. Solubility equilibria are important in pharmaceutical, environmental and many other scenarios.

In chemical thermodynamics, the fugacity of a real gas is an effective partial pressure which replaces the mechanical partial pressure in an accurate computation of chemical equilibrium. It is equal to the pressure of an ideal gas which has the same temperature and molar Gibbs free energy as the real gas.

Robert Minard Garrels was an American geochemist. Garrels applied experimental physical chemistry data and techniques to geology and geochemistry problems. The book Solutions, Minerals, and Equilibria co-authored in 1965 by Garrels and Charles L. Christ revolutionized aqueous geochemistry.

<span class="mw-page-title-main">Geothermobarometry</span> History of rock pressure and temperature

Geothermobarometry is the methodology for estimating the pressure and temperature history of rocks. Geothermobarometry is a combination of geobarometry, where the pressure attained by a mineral assemblage is estimated, and geothermometry where the temperature attained by a mineral assemblage is estimated.

The word hydrolysis is applied to chemical reactions in which a substance reacts with water. In organic chemistry, the products of the reaction are usually molecular, being formed by combination with H and OH groups. In inorganic chemistry, the word most often applies to cations forming soluble hydroxide or oxide complexes with, in some cases, the formation of hydroxide and oxide precipitates.

Chemical WorkBench is a proprietary simulation software tool aimed at the reactor scale kinetic modeling of homogeneous gas-phase and heterogeneous processes and kinetic mechanism development. It can be effectively used for the modeling, optimization, and design of a wide range of industrially and environmentally important chemistry-loaded processes. Chemical WorkBench is a modeling environment based on advanced scientific approaches, complementary databases, and accurate solution methods. Chemical WorkBench is developed and distributed by Kintech Lab.

<span class="mw-page-title-main">Werner Stumm</span> Swiss geochemist (1924–1999)

Werner Stumm was a Swiss chemist. After earning his doctorate in inorganic chemistry at the University of Zürich in 1952 he moved to the U.S. where he was active as a professor at Harvard University until 1969. From 1970 until 1992 he was head of the Swiss Federal Water Resources Centre EAWAG.

The Geochemist's Workbench (GWB) is an integrated set of interactive software tools for solving a range of problems in aqueous chemistry. The graphical user interface simplifies the use of the geochemical code.

Reactive transport modeling in porous media refers to the creation of computer models integrating chemical reaction with transport of fluids through the Earth's crust. Such models predict the distribution in space and time of the chemical reactions that occur along a flowpath. Reactive transport modeling in general can refer to many other processes, including reactive flow of chemicals through tanks, reactors, or membranes; particles and species in the atmosphere; gases exiting a smokestack; and migrating magma.

Aqion is a hydrochemistry software tool. It bridges the gap between scientific software and the calculation/handling of "simple" water-related tasks in daily routine practice. The software aqion is free for private users, education and companies.

<span class="mw-page-title-main">Fred T. Mackenzie</span> American sedimentary biogeochemist (1934–2024)

Frederick T. Mackenzie was an American sedimentary and global biogeochemist. Mackenzie applied experimental and field data coupled to a sound theoretical framework to the solution of geological, geochemical, and oceanographic problems at various time and space scales.

Susan L. Brantley is an American geologist and geochemist who is the Dr. Hubert Barnes and Dr. Mary Barnes Professor at Pennsylvania State University. Her research dominantly studies interactions between fluids and minerals at low temperatures, biological reactions in water-rich fluids within soils, and the geochemical processes that convert rock into soil. However, among many other topics, she has also published work on carbon dioxide emissions from volcanoes, and the environmental impact of shale gas extraction and nuclear waste disposal. During her career, Brantley has published over 200 research papers and book chapters, has been awarded academic prizes and fellowships by many of the world's leading geoscience societies, and has been described as "one of the leading aqueous geochemists of her generation."

<span class="mw-page-title-main">Mark S. Ghiorso</span> American geochemist

Mark S. Ghiorso is an American geochemist who resides in Seattle, Washington. He is best known for creating MELTS, a software tool for thermodynamic modeling of phase equilibria in magmatic systems.

<span class="mw-page-title-main">Dimitri Sverjensky</span> Geochemistry professor

Dimitri Alexander Sverjensky is a professor in Earth and Planetary Sciences at Johns Hopkins University where his research is focused on geochemistry.

Joint Expert Speciation System (JESS) is a package of computer software and data developed collaboratively at Murdoch University and elsewhere by researchers interested in the chemical thermodynamics of water solutions with important applications in industry, biochemistry, medicine and the environment. Using information from the chemical literature, stored in databases for numerous chemical properties, JESS achieves coherence between frequently conflicting sources by automatic methods.

<span class="mw-page-title-main">Groundwater contamination by pharmaceuticals</span> Aquifer contamination by medical drugs

Groundwater contamination by pharmaceuticals, which belong to the category of contaminants of emerging concern (CEC) or emerging organic pollutants (EOP), has been receiving increasing attention in the fields of environmental engineering, hydrology and hydrogeochemistry since the last decades of the twentieth century.

References

  1. Zhu, C. and G. Anderson, 2002, Environmental Applications of Geochemical Modeling. Cambridge University Press, 300 pp.
  2. 1 2 3 4 Bethke, C.M., 2008, Geochemical and Biogeochemical Reaction Modeling. Cambridge University Press, 547 pp.
  3. Garrels, R.M. and F.T. Mackenzie, 1967, Origin of the chemical compositions of some springs and lakes. Equilibrium Concepts in Natural Waters, Advances in Chemistry Series67, American Chemical Society, Washington, DC, pp. 222-242
  4. Helgeson, H.C., 1968, Evaluation of irreversible reactions in geochemical processes involving minerals and aqueous solutions, I. Thermodynamic relations. Geochemica et Cosmochimica Acta32, 853-877
  5. Helgeson, H.C., R.M. Garrels and F.T. Mackenzie, 1969, Evaluation of irreversible reactions in geochemical processes involving minerals and aqueous solutions, II. Applications. Geochemica et Cosmochimica Acta33, 455-481
  6. Zhu, C., 2009, Geochemical Modeling of Reaction Paths and Geochemical Reaction Networks. In E.H. Oelkers and J. Schott(eds.), 2009, Thermodynamics and Kinetics of Water-Rock Interaction. Reviews in Mineralogy and Geochemistry70, 533-569
  7. Brady, P.V. and C.M. Bethke, 2000, Beyond the Kd approach. Ground Water38, 321-322
  8. 1 2 Anderson, G.M. 2009, Thermodynamics of Natural Systems. Cambridge University Press, 664 pp.
  9. Bethke, C.M., ChemPlugin User's Guide Release 15. Aqueous Solutions LLC, Champaign, IL USA https://www.chemplugin.gwb.com/documentation.php
  10. Dethlefsen, Frank; Haase, Christoph; Ebert, Markus; Dahmke, Andreas (2011-01-01). "Effects of the variances of input parameters on water-mineral interactions during CO2 sequestration modeling". Energy Procedia. 10th International Conference on Greenhouse Gas Control Technologies. 4: 3770–3777. doi: 10.1016/j.egypro.2011.02.311 .
  11. Haase, Christoph; Dethlefsen, Frank; Ebert, Markus; Dahmke, Andreas (2013-06-01). "Uncertainty in geochemical modelling of CO2 and calcite dissolution in NaCl solutions due to different modelling codes and thermodynamic databases". Applied Geochemistry. 33: 306–317. doi:10.1016/j.apgeochem.2013.03.001.
  12. 1 2 Haase, Christoph; Ebert, Markus; Dethlefsen, Frank (2016-04-01). "Uncertainties of geochemical codes and thermodynamic databases for predicting the impact of carbon dioxide on geologic formations". Applied Geochemistry. 67: 81–92. doi:10.1016/j.apgeochem.2016.01.008.
  13. Muller, B., 2004, CHEMEQL V3.0, A program to calculate chemical speciation equilibria, titrations, dissolution, precipitation, adsorption, kinetics, pX-pY diagrams, solubility diagrams. Limnological Research Center EAWAG/ETH, Kastanienbaum, Switzerland
  14. van der Lee, J., and L. De Windt, 2000, CHESS, another speciation and complexation computer code. Technical Report no. LHM/RD/93/39, Ecole des Mines de Paris, Fontainebleau
  15. Reed, M.H., 1982, Calculation of multicomponent chemical equilibria and reaction processes in systems involving minerals, gases, and aqueous phase. Geochimica et Cosmochemica Acta46, 513-528.
  16. Steefel, C.I. and A.C. Lasaga, 1994, A coupled model for transport of multiple chemical species and kinetic precipitation/dissolution reactions with application to reactive flow in single phase hydrothermal systems. American Journal of Science294, 529-592
  17. Steefel, C.I., 2001, GIMRT, Version 1.2: Software for modeling multicomponent, multidimensional reactive transport, User's Guide. Report UCRL-MA-143182, Lawrence Livermore National Laboratory, Livermore, California.
  18. Wolery, T.J., 1992a, EQ3/EQ6, a software package for geochemical modeling of aqueous systems, package overview and installation guide (version 7.0). Lawrence Livermore National Laboratory Report UCRL-MA-110662(1).
  19. Shaff, J.E., B.A. Schultz, E.J. Craft, R.T. Clark, and L.V. Kochian, 2010, GEOCHEM-EZ: a chemical speciation program with greater power and flexibility. Plant Soil330(1), 207-214
  20. Bethke, C.M., B. Farrell, and M. Sharifi, 2021, The Geochemist's Workbench Release 15 (five volumes). Aqueous Solutions LLC, Champaign, IL USA
  21. Kulik, D.A., 2002, Gibbs energy minimization approach to model sorption equilibria at the mineral-water interface: Thermodynamic relations for multi-site surface complexation. American Journal of Science302, 227-279
  22. Cheng, H.P. and G.T. Yeh, 1998, Development of a three-dimensional model of subsurface flow, heat transfer, and reactive chemical transport: 3DHYDROGEOCHEM. Journal of Contaminant Hydrology34, 47-83
  23. Westall, J.C., J.L. Zachary and F.F.M. Morel, 1976, MINEQL, a computer program for the calculation of chemical equilibrium composition of aqueous systems. Technical Note 18, R.M. Parsons Laboratory, Department of Civil and Environmental Engineering, Massachusetts Institute of Technology, Cambridge, MA.
  24. Scherer, W.D. and D.C. McAvoy, 1994, MINEQL+, A Chemical Equilibrium Program for Personal Computers, User's Manual, version 3.0. Environmental Research Software, Inc., Hallowell, ME.
  25. Allison, J.D., D.S. Brown and K.J. Novo-Gradac, 1991, MINTEQA2/ PRODEFA2, a geochemical assessment model for environmental systems, version 3.0 user's manual. US Environmental Protectiona Agency Report EPA/600/3-91/021.
  26. "ORCHESTRA | Geochemical and Transport Modelling" . Retrieved 2022-09-29.
  27. Parkhurst, D.L., 1995, User's Guide to PHREEQC, a computer model for speciation, reaction-path, advective-transport and inverse geochemical calculations. US Geological Survey Water-Resources Investigations Report 95-4227.
  28. Parkhurst, D.L. and C.A.J. Appelo, 1999, User's Guide to PHREEQC (version 2), a computer program for speciation, batch-reaction, one-dimensional transport and inverse geochemical calculations. US Geological Survey Water-Resources Investigations Report 99-4259.
  29. Leal, A.M.M., Kulik, D.A., Smith, W.R. and Saar, M.O., 2017, An overview of computational methods for chemical equilibrium and kinetic calculations for geochemical and reactive transport modeling. Pure and Applied Chemistry. 89 (5), 145–166.
  30. Perkins, E.H., 1992, Integration of intensive variable diagrams and fluid phase equilibria with SOLMINEQ.88 pc/shell. In Y.K. Kharaka and A.S. Maest (eds.), Water-Rock Interaction, Balkema, Rotterdam, p. 1079-1081.
  31. Xu, T., E.L. Sonnenthal, N. Spycher and K. Pruess, 2004, TOUGHREACT user's guide: A simulation program for non-isothermal multiphase reactive geochemical transport in variably saturated geologic media. Report LBNL-55460, Lawrence Berkeley National Laboratory, Berkeley, California.
  32. hem.bredband.net/b108693/-VisualMINTEQ_references.pdf
  33. Ball, J.W. and D.K. Nordstrom, 1991, User's manual for WATEQ4F, with revised thermodynamic data base and test cases for calculating speciation of major, trace, and redox elements in natural waters. US Geological Survey Open File Report 91-183.
  34. Tipping E., 1994, WHAM - a chemical equilibrium model and computer code for waters, sediments and soils incorporating a discrete site/electrostatic model of ion-binding by humic substances. Computers and Geosciences20, 973-1023.
  35. "Water Resources Geochemical Software". water.usgs.gov. Retrieved 2020-09-25.