Three dimensional solute transport numerical modeling of salinity fluctuations in a coastal aquifer

E. D.P. Perera, K. Jinno, A. Tsutsumi, Y. Hiroshiro

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Abstract

Over-exploitation of groundwater in coastal aquifers is one of the major reasons for saltwater intrusion. Saltwater intrusion contaminates the freshwater resource in coastal aquifers by increasing the salinity levels of the groundwater. A small quantity of saltwater is enough to contaminate the large quantity of freshwater in a coastal aquifer. Therefore, saltwater intrusion should be considered seriously. Numerical models are needed to enhance the understanding of saltwater intrusion and its related phenomena extensively. The Motooka region in the Fukuoka prefecture in Japan is a coastal area, where groundwater is the main water resource for green house agriculture and domestic use. Over-exploitation of groundwater in this coastal region has resulted in saltwater intrusion and thus in the contamination of the freshwater aquifer. In addition to the more obvious effects of saltwater intrusion, fluctuations in salinity caused by such intrusion is a crucial problem to address, since even slight changes in salinity of the water use for agricultural purposes significantly affect the crop's growth and yield. So far, a research on the salinity fluctuations with groundwater pumping and their effects on seasonal recharge of groundwater in the Motooka region has not been conducted. Therefore, in this study a three-dimensional density-dependent solute transport flow model is developed to simulate the salinity fluctuations due to groundwater pumping. In the present numerical study, the emphasis is on the development of conceptual, mathematical and numerical model of variable density flow and solute transport and its application to simulate the salinity fluctuation due to groundwater pumping at different rates. The model is based on the "transition zone"approach, which considers the interface as a miscible zone where freshwater and saltwater is mixing while maintaining a density gradient across the freshwater/saltwater interface. The transition zone approach requires simultaneous solutions of the governing water flow and solute transport equations. To this end, the model incorporates three fundamental equations in flow and solute transport, namely Darcy's law, general groundwater flow equation and the advection dispersion solute transport equation. The groundwater flow equation and solute transport equation are coupled by the equation of state to produce the salt concentration at each time step for whole flow domain. The finite difference method is used as the numerical technique to solve the partial differential equations of flow and transport under an implicit scheme. The method of characteristics is applied to solve the advection term in the solute transport equation. A non-uniform discretized grid system is adopted in the flow domain allocating relatively small grid sizes to pumping well locations. To achieve reliable results, relevant and important hydro-geological parameters are assigned to the numerical model after considering the hydrological situation of the Motooka region. Different boundary conditions are assigned considering the dominant hydrological processes those are believed to be in effect in the selected area. The numerical results obtained from the model demonstrate the salinity variation due to groundwater pumping and seasonal recharge rates from year 2001 to 2007 under the influence of saltwater intrusion. The results also reveal that model is capable of correctly simulating the physical processes. A comparison of the measured and modeled electric conductivities shows reasonable agreement.

Original languageEnglish
Title of host publication18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation
Subtitle of host publicationInterfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings
Pages4128-4134
Number of pages7
Publication statusPublished - Dec 1 2009
Event18th World IMACS Congress and International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, MODSIM09 - Cairns, QLD, Australia
Duration: Jul 13 2009Jul 17 2009

Publication series

Name18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings

Other

Other18th World IMACS Congress and International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, MODSIM09
CountryAustralia
CityCairns, QLD
Period7/13/097/17/09

Fingerprint

Solute Transport
Solute transport
Salinity
Ground Water
Numerical Modeling
Aquifers
Groundwater
Fluctuations
Three-dimensional
Transport Equation
Numerical models
Groundwater Flow
Groundwater flow
Advection
Exploitation
Water
Non-uniform Grid
Greenhouse
Darcy's Law
Model

All Science Journal Classification (ASJC) codes

  • Computational Theory and Mathematics
  • Computational Mathematics
  • Modelling and Simulation

Cite this

Perera, E. D. P., Jinno, K., Tsutsumi, A., & Hiroshiro, Y. (2009). Three dimensional solute transport numerical modeling of salinity fluctuations in a coastal aquifer. In 18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings (pp. 4128-4134). (18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings).

Three dimensional solute transport numerical modeling of salinity fluctuations in a coastal aquifer. / Perera, E. D.P.; Jinno, K.; Tsutsumi, A.; Hiroshiro, Y.

18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings. 2009. p. 4128-4134 (18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings).

Research output: Chapter in Book/Report/Conference proceedingConference contribution

Perera, EDP, Jinno, K, Tsutsumi, A & Hiroshiro, Y 2009, Three dimensional solute transport numerical modeling of salinity fluctuations in a coastal aquifer. in 18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings. 18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings, pp. 4128-4134, 18th World IMACS Congress and International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, MODSIM09, Cairns, QLD, Australia, 7/13/09.
Perera EDP, Jinno K, Tsutsumi A, Hiroshiro Y. Three dimensional solute transport numerical modeling of salinity fluctuations in a coastal aquifer. In 18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings. 2009. p. 4128-4134. (18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings).
Perera, E. D.P. ; Jinno, K. ; Tsutsumi, A. ; Hiroshiro, Y. / Three dimensional solute transport numerical modeling of salinity fluctuations in a coastal aquifer. 18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings. 2009. pp. 4128-4134 (18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings).
@inproceedings{41715da9c8a8498297039df041a344a8,
title = "Three dimensional solute transport numerical modeling of salinity fluctuations in a coastal aquifer",
abstract = "Over-exploitation of groundwater in coastal aquifers is one of the major reasons for saltwater intrusion. Saltwater intrusion contaminates the freshwater resource in coastal aquifers by increasing the salinity levels of the groundwater. A small quantity of saltwater is enough to contaminate the large quantity of freshwater in a coastal aquifer. Therefore, saltwater intrusion should be considered seriously. Numerical models are needed to enhance the understanding of saltwater intrusion and its related phenomena extensively. The Motooka region in the Fukuoka prefecture in Japan is a coastal area, where groundwater is the main water resource for green house agriculture and domestic use. Over-exploitation of groundwater in this coastal region has resulted in saltwater intrusion and thus in the contamination of the freshwater aquifer. In addition to the more obvious effects of saltwater intrusion, fluctuations in salinity caused by such intrusion is a crucial problem to address, since even slight changes in salinity of the water use for agricultural purposes significantly affect the crop's growth and yield. So far, a research on the salinity fluctuations with groundwater pumping and their effects on seasonal recharge of groundwater in the Motooka region has not been conducted. Therefore, in this study a three-dimensional density-dependent solute transport flow model is developed to simulate the salinity fluctuations due to groundwater pumping. In the present numerical study, the emphasis is on the development of conceptual, mathematical and numerical model of variable density flow and solute transport and its application to simulate the salinity fluctuation due to groundwater pumping at different rates. The model is based on the {"}transition zone{"}approach, which considers the interface as a miscible zone where freshwater and saltwater is mixing while maintaining a density gradient across the freshwater/saltwater interface. The transition zone approach requires simultaneous solutions of the governing water flow and solute transport equations. To this end, the model incorporates three fundamental equations in flow and solute transport, namely Darcy's law, general groundwater flow equation and the advection dispersion solute transport equation. The groundwater flow equation and solute transport equation are coupled by the equation of state to produce the salt concentration at each time step for whole flow domain. The finite difference method is used as the numerical technique to solve the partial differential equations of flow and transport under an implicit scheme. The method of characteristics is applied to solve the advection term in the solute transport equation. A non-uniform discretized grid system is adopted in the flow domain allocating relatively small grid sizes to pumping well locations. To achieve reliable results, relevant and important hydro-geological parameters are assigned to the numerical model after considering the hydrological situation of the Motooka region. Different boundary conditions are assigned considering the dominant hydrological processes those are believed to be in effect in the selected area. The numerical results obtained from the model demonstrate the salinity variation due to groundwater pumping and seasonal recharge rates from year 2001 to 2007 under the influence of saltwater intrusion. The results also reveal that model is capable of correctly simulating the physical processes. A comparison of the measured and modeled electric conductivities shows reasonable agreement.",
author = "Perera, {E. D.P.} and K. Jinno and A. Tsutsumi and Y. Hiroshiro",
year = "2009",
month = "12",
day = "1",
language = "English",
isbn = "9780975840078",
series = "18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings",
pages = "4128--4134",
booktitle = "18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation",

}

TY - GEN

T1 - Three dimensional solute transport numerical modeling of salinity fluctuations in a coastal aquifer

AU - Perera, E. D.P.

AU - Jinno, K.

AU - Tsutsumi, A.

AU - Hiroshiro, Y.

PY - 2009/12/1

Y1 - 2009/12/1

N2 - Over-exploitation of groundwater in coastal aquifers is one of the major reasons for saltwater intrusion. Saltwater intrusion contaminates the freshwater resource in coastal aquifers by increasing the salinity levels of the groundwater. A small quantity of saltwater is enough to contaminate the large quantity of freshwater in a coastal aquifer. Therefore, saltwater intrusion should be considered seriously. Numerical models are needed to enhance the understanding of saltwater intrusion and its related phenomena extensively. The Motooka region in the Fukuoka prefecture in Japan is a coastal area, where groundwater is the main water resource for green house agriculture and domestic use. Over-exploitation of groundwater in this coastal region has resulted in saltwater intrusion and thus in the contamination of the freshwater aquifer. In addition to the more obvious effects of saltwater intrusion, fluctuations in salinity caused by such intrusion is a crucial problem to address, since even slight changes in salinity of the water use for agricultural purposes significantly affect the crop's growth and yield. So far, a research on the salinity fluctuations with groundwater pumping and their effects on seasonal recharge of groundwater in the Motooka region has not been conducted. Therefore, in this study a three-dimensional density-dependent solute transport flow model is developed to simulate the salinity fluctuations due to groundwater pumping. In the present numerical study, the emphasis is on the development of conceptual, mathematical and numerical model of variable density flow and solute transport and its application to simulate the salinity fluctuation due to groundwater pumping at different rates. The model is based on the "transition zone"approach, which considers the interface as a miscible zone where freshwater and saltwater is mixing while maintaining a density gradient across the freshwater/saltwater interface. The transition zone approach requires simultaneous solutions of the governing water flow and solute transport equations. To this end, the model incorporates three fundamental equations in flow and solute transport, namely Darcy's law, general groundwater flow equation and the advection dispersion solute transport equation. The groundwater flow equation and solute transport equation are coupled by the equation of state to produce the salt concentration at each time step for whole flow domain. The finite difference method is used as the numerical technique to solve the partial differential equations of flow and transport under an implicit scheme. The method of characteristics is applied to solve the advection term in the solute transport equation. A non-uniform discretized grid system is adopted in the flow domain allocating relatively small grid sizes to pumping well locations. To achieve reliable results, relevant and important hydro-geological parameters are assigned to the numerical model after considering the hydrological situation of the Motooka region. Different boundary conditions are assigned considering the dominant hydrological processes those are believed to be in effect in the selected area. The numerical results obtained from the model demonstrate the salinity variation due to groundwater pumping and seasonal recharge rates from year 2001 to 2007 under the influence of saltwater intrusion. The results also reveal that model is capable of correctly simulating the physical processes. A comparison of the measured and modeled electric conductivities shows reasonable agreement.

AB - Over-exploitation of groundwater in coastal aquifers is one of the major reasons for saltwater intrusion. Saltwater intrusion contaminates the freshwater resource in coastal aquifers by increasing the salinity levels of the groundwater. A small quantity of saltwater is enough to contaminate the large quantity of freshwater in a coastal aquifer. Therefore, saltwater intrusion should be considered seriously. Numerical models are needed to enhance the understanding of saltwater intrusion and its related phenomena extensively. The Motooka region in the Fukuoka prefecture in Japan is a coastal area, where groundwater is the main water resource for green house agriculture and domestic use. Over-exploitation of groundwater in this coastal region has resulted in saltwater intrusion and thus in the contamination of the freshwater aquifer. In addition to the more obvious effects of saltwater intrusion, fluctuations in salinity caused by such intrusion is a crucial problem to address, since even slight changes in salinity of the water use for agricultural purposes significantly affect the crop's growth and yield. So far, a research on the salinity fluctuations with groundwater pumping and their effects on seasonal recharge of groundwater in the Motooka region has not been conducted. Therefore, in this study a three-dimensional density-dependent solute transport flow model is developed to simulate the salinity fluctuations due to groundwater pumping. In the present numerical study, the emphasis is on the development of conceptual, mathematical and numerical model of variable density flow and solute transport and its application to simulate the salinity fluctuation due to groundwater pumping at different rates. The model is based on the "transition zone"approach, which considers the interface as a miscible zone where freshwater and saltwater is mixing while maintaining a density gradient across the freshwater/saltwater interface. The transition zone approach requires simultaneous solutions of the governing water flow and solute transport equations. To this end, the model incorporates three fundamental equations in flow and solute transport, namely Darcy's law, general groundwater flow equation and the advection dispersion solute transport equation. The groundwater flow equation and solute transport equation are coupled by the equation of state to produce the salt concentration at each time step for whole flow domain. The finite difference method is used as the numerical technique to solve the partial differential equations of flow and transport under an implicit scheme. The method of characteristics is applied to solve the advection term in the solute transport equation. A non-uniform discretized grid system is adopted in the flow domain allocating relatively small grid sizes to pumping well locations. To achieve reliable results, relevant and important hydro-geological parameters are assigned to the numerical model after considering the hydrological situation of the Motooka region. Different boundary conditions are assigned considering the dominant hydrological processes those are believed to be in effect in the selected area. The numerical results obtained from the model demonstrate the salinity variation due to groundwater pumping and seasonal recharge rates from year 2001 to 2007 under the influence of saltwater intrusion. The results also reveal that model is capable of correctly simulating the physical processes. A comparison of the measured and modeled electric conductivities shows reasonable agreement.

UR - http://www.scopus.com/inward/record.url?scp=80053030406&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=80053030406&partnerID=8YFLogxK

M3 - Conference contribution

AN - SCOPUS:80053030406

SN - 9780975840078

T3 - 18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation: Interfacing Modelling and Simulation with Mathematical and Computational Sciences, Proceedings

SP - 4128

EP - 4134

BT - 18th World IMACS Congress and MODSIM09 International Congress on Modelling and Simulation

ER -