Modeling heavy metal water-sediment partitioning using a nonlinear bifurcation approach

Livija Cveticanin; Vesna Gvoić; Dejan Krčmar; Aleksandra Tubić; Miodrag Zuković; Mihaela Timofti; Miljana Prica; Livija Cveticanin; Vesna Gvoić; Dejan Krčmar; Aleksandra Tubić; Miodrag Zuković; Mihaela Timofti; Miljana Prica

doi:10.3934/environsci.2026034

AIMS Environmental Science

2026, Volume 13, Issue 5: 865-895. doi: 10.3934/environsci.2026034

Previous Article Next Article

Research article

Modeling heavy metal water-sediment partitioning using a nonlinear bifurcation approach

1.
Department of Mechanics, Faculty of Technical Sciences, University of Novi Sad, Trg Dositeja Obradovića 6, 21000 Novi Sad, Serbia
2.
GAMF Faculty of Engineering and Computer Science, John von Neumann University, H-6000 Kecskemet, Izsak ut 10, Hungary
3.
Department of Graphic Engineering and Design, Faculty of Technical Sciences, University of Novi Sad, Trg Dositeja Obradovića 6, 21000 Novi Sad, Serbia
4.
Department of Chemistry, Biochemistry and Environmental Protection, Faculty of Sciences, University of Novi Sad, Trg Dositeja Obradovića 3, 21000 Novi Sad, Serbia
5.
Department of Applied Sciences, Cross-Border Faculty, Dunarea de Jos University of Galati, Str. Domnească, nr. 111, Galati, Romania
6.
REXDAN Research Infrastructure, Dunarea de Jos University of Galati, Str. Domnească, nr. 111, Galati, Romania

Received: 27 July 2026 Revised: 24 August 2026 Accepted: 15 September 2026 Published: 21 September 2026

Heavy metal transport between the water column and sediment is commonly described using linear or empirical relationships, which may not capture nonlinear feedback mechanisms and abrupt changes in contamination states. In this study, we proposed a nonlinear cubic equilibrium model for describing metal partitioning between water and sediment and for identifying critical transition conditions associated with bifurcation behavior.
The proposed model established a nonlinear relationship between the input metal concentration in water I_i and the accumulated sediment concentration C_i. Using experimental data for copper (Cu) and nickel (Ni) transport from the literature data, model parameters were estimated by the least-squares method. The obtained response curves exhibited an S-shaped structure, indicating the existence of multiple equilibrium states, bistability, and hysteresis. Fold points obtained from the condition dI_i/dC_i = 0 defined critical thresholds at which saddle-node bifurcations occur and abrupt transitions between contamination states become possible.
The bifurcation analysis with respect to the nonlinear parameter p_i further identified critical values separating bistable and monostable regimes. The results demonstrated that Cu and Ni systems may exhibit history-dependent behavior, where the sediment contamination level is not determined solely by the current pollution input but also by the previous evolution of the system.
The proposed framework extends conventional pollutant transport models by incorporating nonlinear dynamics, multistability, and regime-shift mechanisms. It provides a mathematical basis for identifying tipping points in water-sediment contamination systems and offers a potential tool for assessing environmental risks associated with sudden changes in metal accumulation.
The main contribution of this work is the demonstration that nonlinear bifurcation analysis provides insight into critical states, multistability, hysteresis, and regime transitions in water-sediment contamination systems, extending conventional concentration-based approaches toward a deeper understanding of nonlinear environmental dynamics.
- heavy metal transport,
- water-sediment interaction,
- nonlinear cubic model,
- bifurcation analysis,
- saddle-node bifurcation,
- hysteresis,
- multistability,
- tipping point,
- regime shifts,
- sediment contamination
Citation: Livija Cveticanin, Vesna Gvoić, Dejan Krčmar, Aleksandra Tubić, Miodrag Zuković, Mihaela Timofti, Miljana Prica. Modeling heavy metal water-sediment partitioning using a nonlinear bifurcation approach[J]. AIMS Environmental Science, 2026, 13(5): 865-895. doi: 10.3934/environsci.2026034

Related Papers:

Abstract

Heavy metal transport between the water column and sediment is commonly described using linear or empirical relationships, which may not capture nonlinear feedback mechanisms and abrupt changes in contamination states. In this study, we proposed a nonlinear cubic equilibrium model for describing metal partitioning between water and sediment and for identifying critical transition conditions associated with bifurcation behavior.

The proposed model established a nonlinear relationship between the input metal concentration in water I_i and the accumulated sediment concentration C_i. Using experimental data for copper (Cu) and nickel (Ni) transport from the literature data, model parameters were estimated by the least-squares method. The obtained response curves exhibited an S-shaped structure, indicating the existence of multiple equilibrium states, bistability, and hysteresis. Fold points obtained from the condition dI_i/dC_i = 0 defined critical thresholds at which saddle-node bifurcations occur and abrupt transitions between contamination states become possible.

The bifurcation analysis with respect to the nonlinear parameter p_i further identified critical values separating bistable and monostable regimes. The results demonstrated that Cu and Ni systems may exhibit history-dependent behavior, where the sediment contamination level is not determined solely by the current pollution input but also by the previous evolution of the system.

The proposed framework extends conventional pollutant transport models by incorporating nonlinear dynamics, multistability, and regime-shift mechanisms. It provides a mathematical basis for identifying tipping points in water-sediment contamination systems and offers a potential tool for assessing environmental risks associated with sudden changes in metal accumulation.

The main contribution of this work is the demonstration that nonlinear bifurcation analysis provides insight into critical states, multistability, hysteresis, and regime transitions in water-sediment contamination systems, extending conventional concentration-based approaches toward a deeper understanding of nonlinear environmental dynamics.

References

[1]	Bai J, Zhao J, Zhang Z, et al. (2022). Assessment and a review of research on surface water quality modeling. Ecol Model 466: 109888.https://doi.org/10.1016/j.ecolmodel.2022.109888 doi: 10.1016/j.ecolmodel.2022.109888
[2]	Wang Q, Li S, Jia P, et al. (2013) A review of surface water quality models. Sci World J 231768.https://doi.org/10.1155/2013/231768 doi: 10.1155/2013/231768
[3]	Fu B, Merritt WS, Croke BFW, et al. (2019) Catchment scale water quality and erosion models: Review and synthesis. Environ Model Softw 114: 75-97.https://doi.org/10.1016/j.envsoft.2018.12.008 doi: 10.1016/j.envsoft.2018.12.008
[4]	Geng M, Qian Z, Jiang H, et al. (2024) Assessing the impact of water-sediment factors on water quality to guide river-connected lake water environment improvement. Sci Total Environ 912: 168866.https://doi.org/10.1016/j.scitotenv.2023.168866 doi: 10.1016/j.scitotenv.2023.168866
[5]	Zhou Y, Wang X, Li W, et al. (2023) Water quality evaluation and pollution source apportionment of surface water in a major city in Southeast China using multi-statistical analyses and machine learning models. Int J Environ Res Public Health 20: 881.https://doi.org/10.3390/ijerph20010881 doi: 10.3390/ijerph20010881
[6]	Xu X, Lai T, Jahan S, et al. (2022) Water and sediment analyse using predictive models. arXiv arXiv: 2203.03422.
[7]	Pimiento MA, Anta J, Torres A (2025) Machine learning predictive modelling for sediment risk indices within an urbanized river channel. Hazard Adv 18: 100708.https://doi.org/10.1016/j.hazadv.2025.100708 doi: 10.1016/j.hazadv.2025.100708
[8]	Lai T, Farid F, Kuan Y, et al. (2025) Transfer learning for assessing heavy metal pollution in seaports sediments. arXiv arXiv: 2506.22096.
[9]	Luviano Soto I, Concha Sánchez Y, Raya A (2025) Water quality polluted by total suspended solids classified within an artificial neural network approach. Water Qual Res J 60: 214-228.https://doi.org/10.2166/wqrj.2024.061 doi: 10.2166/wqrj.2024.061
[10]	Yan T, Shen S-L, Zhou A (2022) Indices and models of surface water quality assessment: Review and perspectives. Environ Pollut 308: 119611.https://doi.org/10.1016/j.envpol.2022.119611 doi: 10.1016/j.envpol.2022.119611
[11]	Shanbehzadeh S, Dastjerdi MV, Hassanzadeh A, et al. (2014) Heavy metals in water and sediment: A case study of Tembi River. J Environ Public Health 858720.https://doi.org/10.1155/2014/858720 doi: 10.1155/2014/858720
[12]	Prajapati S, Sabokruhie P, Brinkmann M, et al. (2023) Modelling transport and fate of copper and nickel across the South Saskatchewan River using WASP-TOXI. Water 15: 265.https://doi.org/10.3390/w15020265 doi: 10.3390/w15020265
[13]	Nayfeh AH, Mook DT (1975) Nonlinear Oscilations, Eds., New York: John Wiley & Sons.
[14]	Majd M. M, Kordzadeh-Kermnai V, Ghalandari V et al. (2022) Adsorption isotherm models: A comprehensive and systematic review (2010-2020). Sci Total Environ 812: 151334.https://doi.org/10.1016/j.scitotenv.2021.151334 doi: 10.1016/j.scitotenv.2021.151334
[15]	Wang J, Guo X (2020) Adsorption isotherm models: Classification, physical meaning, application and solving method. Chemosphere 258: 127279.https://doi.org/10.1016/j.chemosphere.2020.127279 doi: 10.1016/j.chemosphere.2020.127279
[16]	Kornijów R (2026) Asymmetric responses to resource overload: A bow-tie perspective on lake eutrophication and human obesity. Theory Bioscihttps://doi.org/10.1007/s12064-026-00473-y
[17]	Zhang H, Li H, Huo S, et al (2025) Navigating tipping points: A complex systems framework for anticipating lake ecosystem collapse. Fundam Reshttps://doi.org/10.1016/j.fmre.2025.07.004
[18]	USEPA. Method 3051A: Microwave Assisted Acid Digestion of Sediments, Sludges, Soils, and Oils. SW-846. U. S. Environmental Protection Agency, 2007. Available from: https://www.epa.gov/hw-sw846/sw-846-test-method-3051a-microwave-assisted-acid-digestion-sediments-sludges-soils-and
[19]	USEPA. Method 7000B: Flame Atomic Absorption Spectrophotometry. SW-846. U. S. Environmental Protection Agency, 2007. Available from: https://www.epa.gov/hw-sw846/sw-846-test-method-7000b-flame-atomic-absorption-spectrophotometry
[20]	USEPA. Method 7010: Graphite Furnace Atomic Absorption Spectrophotometry, Revision 0. SW-846: Test Methods for Evaluating Solid Waste, Physical/Chemical Methods. U. S. Environmental Protection Agency, 2007. Washington, DC, USA. Available from: USEPA. Method 7010: Graphite Furnace Atomic Absorption Spectrophotometry. SW-846. U. S. Environmental Protection Agency, 2007. Available from: https://www.epa.gov/hw-sw846/sw-846-test-method-7010-graphite-furnace-atomic-absorption-spectrophotometry

Reader Comments

Your name:*

Email:*
© 2026 the Author(s), licensee AIMS Press. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0)