A methodological audit of regional hydrogeochemical literature reveals a systemic gap. Legacy environmental assessments for the Kryvyi Rih Iron Ore Basin rely on frameworks designed for fresh water, which fail when applied to highly mineralized mine waters. Standard equilibrium algorithms are inadequate for environments with extreme ionic strength.

This limitation extends to mass transport: conventional modeling overlooks density-driven convection and gravitational settling of heavy brines, leading to inaccurate spatial migration models.

Overcoming this inertia requires a return to the fundamental laws of chemical thermodynamics and multiphase flow physics, aligning predictive models with actual environmental constraints.

👇 The link to download the updated monograph (v0.80, Open Access, in Ukrainian):
https://doi.org/10.5281/zenodo.20709306

#KryvyiRihBasin #IronOre #MineWater #Hydrogeology #PHREEQC #SvystunovaGully #Geochemistry #Contamination #RStats #Groundwater #DataScience

🔥 Dear colleagues! I present to you the new version (0.80) of my monograph dedicated to the hydrogeochemical (thermodynamic) modeling of the contamination plume around the storage pond in the Svystunova Gully (in Ukrainian).
I have significantly expanded the theoretical and practical sections, adding a description of the hydrodynamic properties and behavior of highly mineralized plumes.

The comprehensive concept of the gradual advancement of the contamination plume and the resulting technogenic karst is now revealed.

❗ In the section 1.2 I critically analyze the methodologies and conclusions of previous studies on the technogenic water bodies of the Kryvbas and the propagation processes of contamination plumes from these sites.

The link to the new version of the monograph:
https://doi.org/10.5281/zenodo.20709306

#phreeqc #SvystunovaGully #geochemistry #minewater #contamination #RStats #groundwater #DataScience #IronOreMining #GIS #EnvironmentalDataScience

Updating the graphical materials for the new edition of the monograph on the thermodynamic modeling of processes in technogenic waters. The figures present pe-pH stability (Pourbaix) diagrams for the vanadium and zinc systems, calculated using the PHREEQC core engine. Unlike simplified "water–element" systems, these models compute the thermodynamic equilibrium of the comprehensive solution matrix, accounting for actual macrocomponent concentrations, ionic strength, and ligand competition at every point of the computational space. Utilizing multicomponent calculations allows for the analysis of physicochemical processes in real water samples, rather than simply reproducing idealized "textbook diagrams".

#Geochemistry #PHREEQC #RStats #MineWater #EnvironmentalScience #DataScience #SvystunovaGully #RedModRphree #FOSS

Instalando PHREEQC en Linux

Un programa que utilizo habitualmente durante el desarrollo de mi investigación es PHREEQC. Sirve para hacer modelado hidrogeoquímico para el estudio de aguas naturales o para el estudio de soluciones en laboratorio. Yo la he utilizado principalmente para preparar soluciones con unas sobresaturaciones concretas, para estudiar la evolución de formulaciones de fertilización foliar con diferentes concentraciones minerales y para estudiar la evolución de una solución en contacto con el CO2 atmosférico, por mencionar solo algunos ejemplos. Pero el software es mucho más potente y yo solo sé utilizar una pequeña parte (poco más que la punta del iceberg que es este programa). En Windows hay dos versiones, la versión con una bonita interfaz gráfica (que no me gusta) y la versión integrada en Notepad++ (que es la que he utilizado siempre). La versión de Notepad++ me permite tener un flujo de trabajo que me resulta muy cómodo, el código aparece con un código de colores para su fácil identificación y, además, tiene la función de autocompletar con el tabulador. ¿Y en Linux, qué opciones tenemos? PHREEQC no ofrece una interfaz gráfica como en Windows, pero he encontrado un modo de tener algo similar.

[...]

https://escritura.social/doclomieu/instalando-phreeqc-en-linux

PHREEQC Version 3

PHREEQC Version 3 is a computer program for speciation, batch-reaction, one-dimensional transport, and inverse geochemical calculations.

USGS

A slightly tongue-in-cheek AI animation of a conceptual model showing leakage of highly mineralized mine waters into a carbonate aquifer:
https://youtu.be/WDIze8QNMYs

Yes - it’s simplified. Yes - real systems are messier.
But sometimes a "clean" visual helps to communicate density-driven flow and reactive processes more intuitively than 20 pages of equations.

The animation idea emerged while working on my monograph. Draft versions are available here:
https://zenodo.org/records/16741148

Engineering humor is allowed - even in hydrogeochemistry 😁

#Hydrogeochemistry #Groundwater #MineWater #ReactiveTransport #CarbonateAquifer #PHREEQC #AIinScience #Geoscience #SvystunovaGully

How reliable are dilution estimates in contaminated aquifers?

In my current monograph, I validated the dilution coefficient using two fully independent approaches:
- chloride concentrations as a conservative tracer
- solution density calculated in PHREEQC (Laliberté–Cooper model)

I then performed a robust statistical comparison of both dilution estimates (n ≈ 2000).
The result: no statistically significant difference between the two methods (Mann–Whitney test), with very close median values and negligible effect size.

❗What does this mean in practice?
It confirms that density-based calculations — often ignored in groundwater contamination studies — provide robust, independent validation of classical tracer-based approaches.
This is especially important for highly mineralized mine waters, where density-driven processes (e.g. buoyancy effects, Rayleigh–Taylor instability) can critically affect plume behavior.

#Groundwater #Geochemistry #PHREEQC #MineWater #Hydrogeology #RStats #SvystunovaGully

Thousands of studies describe Rayleigh–Taylor instability, and most hydrogeological models acknowledge its importance in contaminant transport.

However, in studies of the Svystunova gully mine-water storage pond, I have not seen density explicitly accounted for — despite highly mineralized waters where this factor can fundamentally change plume behaviour.
By calculating solution density within the contamination halo, I demonstrated a real risk of gravity-driven sinking of the contamination plume to the base of the aquifer. This finding adds a critical dimension to impact assessment and is now documented as a standalone section in my ongoing report.

If you work with mine water, aquifers contamination, or long-term contamination modelling, this is a parameter worth revisiting.

#Hydrogeochemistry #Groundwater #MineWater #PHREEQC #RStats #QGIS #EnvironmentalRisk #Contamination #RiskAnalysis #SvystunovaGully

Based on the metasomatic (thermodynamic) zonation of the groundwater contamination plume identified in my study, I prepared a schematic impact map.

Unlike conventional chloride concentration isolines (which are often weakly informative), this approach reflects the thermodynamic state of the aquifer. The map clearly shows that the aggressive corrosion-by-mixing zone (Subzone III) occupies a large area and, in several places, extends beyond the existing monitoring well network.

This indicates that the current monitoring system is insufficient to delineate the full zone of influence of the mine water impoundment.

The criteria for defining metasomatic subzones and the associated risks are described in detail in the current published version of my monograph:
https://zenodo.org/records/16741148

#PHREEQC #RStats #QGIS #Geochemistry #Groundwater #Contamination #Hydrogeology #EnvironmentalGeochemistry #SvystunovaGully

I’m currently open to opportunities related to environmental geochemistry, groundwater, and geochemical modeling.

My background includes PHREEQC-based modeling, reproducible workflows (R, QGIS), environmental data analysis, applied machine learning (exploratory / environmental) and applied research on mine-water impacts on carbonate aquifers.

An example of my applied research and methodology (open access):
https://zenodo.org/records/16741148

Based in Canada (Calgary). Open to discussions and professional connections.

#Hydrogeochemistry #Groundwater #EnvironmentalGeochemistry
#PHREEQC #GeochemicalModeling #ReproducibleResearch #Canada #Alberta #Calgary #YYC

Гiдрогеохiмiчне моделювання впливу фiльтрацiйних втрат високомiнералiзованих шахтних вод зi ставка-накопичувача у балцi Свистунова (версiя 0.75)

Ставок-накопичувач високомінералізованих шахтних вод у балці Свистунова було побудовано у 1976 р. Останні два десятиліття ставок щорічно приймає 10-13 млн/куб.м шахтних вод від південної групи шахт Криворізького залізорудного басейну. В подальшому накопичені води скидаються в ріку Інгулець. При цьому, частина вод втрачається через днище ставка, дренуючись у водоносний горизонт відкладів міоцену. За період існування ставка-накопичувача в навколишніх породах сформувався потужний ореол забруднення, який призвів до суттєвих змін гідрогеологічного режиму та хімічного складу водоносного горизонту. На південь від комплексу шламосховищ та ставка-накопичувача практично відсутні питні води четвертинного та більш захищеного неогенового водоносних горизонтів. При цьому - фактична межа контуру забруднення не встановлена. Проблема операційної діяльності ставка-накопичувача потребує багатовекторного комплексного підходу, оскільки поєднує гідрогеохімічні, екологічні, інженерні, економічні та соціальні чинники, які взаємодіють у просторі та часі. Це вимагає системного аналізу для адекватного оцінювання ризиків та прийняття рішень. Метою представленого дослідження є моделювання довготривалого геохімічного впливу фільтраційних втрат високомінералізованих шахтних вод зі ставка-накопичувача на водоносний горизонт, який залягає у вапнякових відкладах. Для досягнення поставленої мети проведено цикл математичного моделювання гідрогеохімічних процесів у горизонті підземних вод та подальший аналіз отриманої геопросторової інформації. Основним методом дослідження є обчислення термодинамічної рівноваги розчинених у горизонті підземних вод компонентів. Особливу увагу приділено рівновазі карбонатних мінералів: кальциту (Calcite), арагоніту (Aragonite), невпорядкованого доломіту (Dolomite(disordered)) та магнезиту (Magnesite). Процеси розчинення та доломітизації вапняків високомінералізованими та розбавленими водами - відомий факт, який описаний у численній науковій літературі. Ці процеси спричиняють не лише зміну пористості та проникності вапнякових товщ, але й суттєво впливають на переведення токсичних мікроелементів-домішок у міграційні форми. Дослідження базується на результатах хімічного аналізу проб води, відібраних із випускної труби, ставка-накопичувача та навколишніх спостережних свердловин упродовж березня 2012-листопада 2021 та березня 2023 року. Автор усвідомлює ймовірні методичні обмеження та недоліки у процесах відбору та хімічного аналізу проб води. Однак ці дані є єдиною доступною емпіричною базою, на якій формуються управлінські рішення та оцінюється стан довкілля. Саме це робить їх легітимним об'єктом для подальшого наукового аналізу та моделювання. Окремо автор зауважує на тому, що існуюча мережа спостережних свердловин не охоплює всю зону впливу ставка-накопичувача - навіть у найвіддаленіших свердловинах не досягається термодинамічна рівновага суміші вод із масивом водовміщуючих вапняків. На момент проведення дослідження автору не було відомо про жодні відкриті роботи, в яких використовувались би методи термодинамічних розрахунків для підземних вод, що зазнають впливу шахтних вод зі ставка-накопичувача у балці Свистунова. З огляду на дуже обмежене поширення програмного забезпечення типу PHREEQC серед українських фахівців та відсутність публікацій, присвячених рівноважному геохімічному моделюванню, можна з високою ймовірністю припустити, що подібні розрахунки раніше не виконувалися або залишалися в межах закритих відомчих звітів. Таким чином, представлена робота може розглядатися як перша спроба відкритого й прозорого аналізу умов рівноваги підземних вод із гірськими породами у зоні техногенного впливу. На відміну від класичних описових методів, представлене дослідження зосереджується на термодинамічному аналізі потенційної розчинності мінералів, змінах індексів насичення та оцінці поведінки макро- і токсичних мікроелементів. Такий підхід дозволяє виявити ймовірні процеси, які не фіксуються безпосередньо в результатах моніторингу, але можуть чинити критичний вплив на стан водоносної системи. У рамках цього дослідження автор пропонує інтерпретацію процесів забруднення водоносного горизонту як геохімічного аналога інфільтраційного метасоматозу. Такий підхід дозволяє вийти за межі описової фіксації підвищених концентрацій та розглядати інфільтрацію шахтних вод як активне заміщення мінерального складу середовища в умовах відкритої гідродинамічної системи. Також це дозволяє суттєво розширити класичну дифузійну модель, яка розглядає лише механічне <<розповзання>> забруднюючих речовин у горизонті підземних вод. Натомість ореол забруднення розглядається як складно організована гідрогеохімічна система, в межах якої формуються реакційні зони з характерними мінеральними й хімічними перетвореннями. За результатами термодинамічного моделювання та подальшого аналізу даних, в ореолі забруднення простежується зональність, наявність реакційних фронтів, процеси мобілізації та іммобілізації хімічних елементів, а також утворення вторинних мінералів у результаті локального перенасичення. Таким чином, зона впливу ставка-накопичувача може розглядатися як своєрідний метасоматичний пояс, який поступово змінює властивості водоносного горизонту та, що є найважливішим, породного масиву. Для дослідження потенційних ризиків міграції токсичних мікроелементів було здійснено гідрогеохімічне моделювання можливих форм міграції кадмію, стронцію, свинцю, цинку та ванадію. Результати моделювання вказують на високу міграційну здатність цих елементів як у воді ставка-накопичувача, так і в зоні змішування первинних та високомінералізованих вод у водоносному горизонті. Це дозволяє прогнозувати подальше їхнє потрапляння до гідрографічної мережі. Автор пропонує новий підхід до визначення межі зони техногенного впливу ставка-накопичувача на водоносний горизонт. На відміну від традиційної практики, яка базується на досягненні фонових або граничних концентрацій певних компонентів, запропонований підхід ґрунтується на показниках термодинамічної рівноваги між підземними водами та вміщуючими гірськими породами. Порушення цієї рівноваги, зокрема з кальцитом та іншими карбонатними мінералами, розглядається як фундаментальна ознака геохімічної деструкції середовища. Відповідно, справжня межа зони впливу повинна фіксуватись не там, де <<все виглядає нормально за концентраціями>>, а там, де система переходить від зміненого до збалансованого первинного стану. Особливу увагу автор приділяє якісній візуалізації результатів, що є запорукою швидкого орієнтування в багатофакторних даних і кращого розуміння просторово-геохімічної структури ореолу забруднення. Графічне подання не лише доповнює текстовий аналіз, але й виступає аналітичним інструментом, який дозволяє виявляти закономірності, аномалії, діагностично значущі кластери та межі. Візуальна перевірка даних у різних координатних проекціях, з фільтрацією або без, є важливою частиною авторського підходу до верифікації та інтерпретації отриманих результатів. Наразі, гідрогеохімічні процеси які відбуваються із високомінералізованими шахтними водами, можна поділити на декілька самостійних фаз: фаза змішування вод різного походження в системах гідрозахисту шахтних полів, фаза змішування води з різних шахтних полів та часткове відстоювання у проміжних резервуарах, фаза потрапляння та відстоювання води безпосередньо у ставку-накопичувачі, фаза змішування води із прісною водою карбонатного водоносного горизонту та фаза змішування частково <<деактивованої>> води із водою ріки Інгулець після дренування карбонатного водоносного горизонту мережею джерел у долині ріки. Окремо можна виділити фазу змішування в період скиду шахтних вод зі ставка-накопичувача. Кожна із описаних фаз має свій окремій набір термодинамічних параметрів та можливих гідрогеохімічних процесів. У даному дослідженні розглядається лише фаза перетворень втрачених зі ставка-накопичувача вод у карбонатному водоносному горизонті. Виявлені в ході дослідження закономірності вказують на те, що ореол забруднення являє собою двоступеневу зональну структуру. Внутрішня та зовнішня зони пов'язані між собою механізмом позитивного зворотного зв'язку. Така конфігурація свідчить про самопідсилення процесів поширення забруднення в поточних умовах гідрогеохімічної взаємодії. Ліквідація втрат води зі ставка-накопичувача не припинить процеси деградації масиву вапняків. Аналіз ризиків у межах цього дослідження базується передусім на вивченні геохімічних процесів, які визначають напрямки еволюції системи та потенційні сценарії розвитку подій іншого характеру - фізичного, економічного, соціального або екологічного. На переконання автора, саме геохімічні (фізико-хімічні) чинники формують початкові умови, від яких безпосередньо залежать подальші зміни в інших підсистемах. Відтак геохімія розглядається як вихідна база для міждисциплінарного аналізу та прогнозування. Формування переліку літературних джерел здійснювалося на основі персонального ознайомлення автором із матеріалами. Сучасні генеративні моделі штучного інтелекту дають можливість скласти величезний перелік дотичних до теми літературних джерел, у яких містяться підтвердження висловлених тез та гіпотез. Проте автор обрав <<класичний>> підхід до ознайомлення та обробки інформації з результатів попередніх досліджень, який полягає в індивідуальному ознайомленні з літературними джерелами та їхніми описами. Дослідження виконано на стику кількох наукових напрямів: гідрогеохімії, просторового та статистичного аналізу, чисельного моделювання та теорії метасоматичної зональності. Автор усвідомлює, що така багатошарова структура підходу може виглядати нетипово для класичних дисциплінарних рамок. Водночас, саме така інтеграція дозволяє виявити приховані закономірності та пропонувати нові інтерпретації складних природно-техногенних процесів. Основний фокус утримано не на ортодоксальності, а на пояснювальній силі моделі, побудованої на основі реальних даних. Частина запропонованих інтерпретацій є гіпотетичними та не претендує на остаточну істину. Їхня цінність полягає у спробі системного узагальнення і формулювання нових підходів, які можуть бути перевірені при подальших дослідженнях. Незважаючи на відсутність попередніх академічних публікацій автора, ця монографія є результатом багаторічної дослідницької роботи, здійсненої в умовах повної незалежності - без інституційної підтримки, грантів чи доступу до академічних баз. Структура, стиль та наукова аргументація обрані свідомо, як спроба адаптувати класичний формат монографії до умов самоорганізованого дослідження. Цей підхід не імітує академічність - він прагне зберегти її сутність у нових реаліях незалежних дослідницьких ініціатив. Після публікації цієї монографії у версії 0.72 спостерігається активізація певних організацій, які раніше не виявляли зацікавленості у системному аналізі даних стосовно проблеми ставка-накопичувача у балці Свистунова. На думку автора, це може свідчити про усвідомлення актуальності розкритих у цій роботі підходів. У цьому сенсі поява альтернативних гіпотез чи повторних досліджень слід розглядати не як критику, а як індикатор впливовості та резонансу незалежної ініціативи.

Zenodo

#SvystunovaGully
I’ve published a short methodological paper on Zenodo, derived from my ongoing monograph on mine-water impact on carbonate aquifers.

The note introduces the pre-dissolution effect — a positive feedback mechanism where aggressive saline waters accelerate carbonate dissolution, increasing permeability and prolonging contamination even after the source is isolated.

The focus is not on a single case study, but on a transferable geochemical interpretation framework (PHREEQC-based, reproducible, field-data driven).

📄 Open access: https://zenodo.org/records/18142106

#Hydrogeochemistry #EnvironmentalGeochemistry #PHREEQC #Groundwater #ReproducibleResearch #RStats #QGIS #FOSS

Pre-dissolution effect: accelerated contamination of a carbonate aquifer by mine waters

High-salinity mine water infiltration into carbonate aquifers may enhance contaminant migration due to the pre-dissolution of carbonate minerals. This short article presents the results of PHREEQC thermodynamic modeling of the Svystunova Gully mine water impoundment (Ukraine). The calculations demonstrate that aggressive mine waters, initially undersaturated with respect to carbonates, induce local dissolution, thereby potentially increasing permeability. This, in turn, further accelerates the spread of contamination. Pre-dissolution is therefore identified as a key mechanism explaining both the intensity and extent of groundwater contamination in carbonate environments.

Zenodo