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

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

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

Seven months after publishing the draft of my monograph on Zenodo, it has passed 2,000 downloads.

I take this as a signal of ongoing interest in the topic and in the proposed methodological approach. It also suggests that the framework may be applicable beyond the original case study.

Currently working on the next version with further refinements and extensions.

#Hydrogeochemistry #PHREEQC #GeochemicalModeling #ReproducibleResearch #EnvironmentalGeochemistry #rstats #qgis #SvystunovaGully #mining

🧪 Davis vs Pitzer: stress-testing a geochemical model

In hydrogeochemical modeling, there is a common concern:
are “simpler” activity models reliable for saline waters, or do they break down beyond their formal limits?

In my study, I primarily use minteq.v4.dat (Davies equation) in PHREEQC. Formally, Davies is recommended up to ionic strength ~0.5–0.7, while Pitzer theory is considered the gold standard for brines.

Instead of assuming, I tested it.

What I did
- Recalculated >1000 real water samples
- Compared saturation indices (SI) computed with:
- Davies-based MINTEQ.V4
- Pitzer formulation (pitzer.dat)

What I found
For calcite, results are nearly identical:
- points lie close to the 1:1 line
- median SI differs by <0.1
- no statistically significant difference (Mann–Whitney p = 0.22)

This is not an argument against Pitzer.
It is evidence that Davies-based models can remain robust for certain systems when their

#Hydrogeochemistry #PHREEQC #GeochemicalModeling #ReproducibleResearch #SvystunovaGully

The “Hydrochemical Gap”: why manual saturation indices don’t work
Many Kryvbas studies still estimate calcite saturation as if high-salinity waters behaved like distilled water.
This plot shows the opposite.

🔵 Blue line: textbook Ksp for calcite.
🔴 Red points: actual ion activity product (IAP) for real samples (PHREEQC).

Key facts:

1️⃣ At 30–40 g/L salinity, real IAP is 1.5–1.6 log units lower — meaning free Ca²⁺ activity is ~38× lower than concentration-based formulas predict.

2️⃣ The dip at 2–5 g/L marks the mixing zone, where ionic strength and ion-pairing change non-linearly and break simple calculations.

3️⃣ Main reason:
• high ionic strength suppresses activity,
• Ca²⁺ forms complexes,
• multicomponent interactions are non-linear.

Result:
Manual SI calculations suggest “precipitation,” but real activity shows undersaturation and aggressive water.
Only thermodynamic modelling reveals the true behaviour.

#Hydrogeochemistry #Geochemistry #Mining #PHREEQC #RStats #SvystunovaGully #Kryvbas #SvystunovaGully

🧪 The “Oversaturation Illusion” in Kryvbas Mine Waters

While modeling Kryvbas water chemistry (R + PHREEQC), I found a fundamental issue in how saturation is often evaluated.

We usually calculate calcite equilibrium from ion concentrations — fine for fresh water.
But Kryvbas mine waters are brines, where ionic strength and complexation dominate.

📉 Results from ~1000 samples (minteq.v4):
- Once salinity exceeds ~3 g/L, Ca²⁺ activity drops sharply.
- At 15–20 g/L, calcium activity coefficient is ≈ 0.35.

Meaning: more than half of the “calcium concentration” is inert — a dead load that cannot form precipitates.

This explains why traditional methods predicted oversaturation where the water was actually aggressive and dissolving rocks.

Modeling: PHREEQC + minteq.v4 (US EPA), Davis equation.

#Hydrogeochemistry #WaterChemistry #PHREEQC #Geochemistry #Groundwater
#Mining #Tailings #IonActivity #Thermodynamics #Kryvbas #OpenScience #RStats #SvystunovaGully

🎉 A small milestone for my hydrogeochemistry monograph about #SvystunovaGully

Today, my Zenodo record passed 1,500 downloads across both released versions — something I never expected when I began this project as a personal, curiosity-driven exploration.

It’s an independent study on groundwater–rock interaction, technogenic metasomatism, and carbonate–water equilibria, based largely on PHREEQC thermodynamic modelling and geospatial data.

I’m grateful to everyone who took an interest, skimmed, downloaded, or shared it.
Your attention gives this work a life I didn’t imagine it would ever have.

❗ And one more thing: in the monograph’s preamble I explicitly thank the developers of free and open-source software. Without their tools — R, PHREEQC, QGIS, LaTeX, Linux, JabRef and many others — this research would simply not have been possible.

#Hydrogeochemistry #Geochemistry #Groundwater #PHREEQC #FOSS #OpenScience #RStats #QGIS #Zenodo #EnvironmentalScience #Carbonates #Metasomatism #GeoData #WaterPollution

💧 Aggressive groundwater and the long memory of contamination

Even after more than a decade of monitoring, the groundwater in this carbonate aquifer remains undersaturated with respect to calcite — consistently showing negative SI values across nearly all observation wells.

This means that the system is still chemically aggressive toward limestone, slowly dissolving the host rock.
The impoundment, operating since 1976, continues to influence the aquifer — with the geometry of possible subsurface voids still uncertain to this day.

📘 All supporting data and modeling results are included in the draft monograph:
🔗 https://zenodo.org/records/16741148

🧪 Data & visualization: PHREEQC + R + QGIS

#Hydrogeochemistry #Geochemistry #MineWater #Groundwater #PHREEQC #EnvironmentalGeoscience #IndependentResearch #OpenScience #QGIS #RStats #Aquifer #Zenodo #SvystunovaGully #GroundwaterContamination #FOSS