Numerical modelling of nitrate transport in fractured porous media under non-isothermal conditions

Subsurface contamination is a frequent occurrence in fractured porous systems, posing a potential threat for the groundwater contamination. Tracking the movement of these contaminants is an inherent aspect of effective remediation strategy. The non-isothermal conditions prevailing in the subsurface environment further add to the complexity of the existing scenario. The current study focuses on simulating the concentration profiles of nitrogen species in a fracture-matrix system under non-isothermal conditions. The kinetics and biochemical thermodynamics of nitrogen transformation reactions were explicitly modelled in this study by adopting a finite differential numerical scheme. The numerical results clearly depicted the spatial-temporal profiles of the concentration of all the species in response to the observed peak values. Considering the sensitivity of the model parameters, an increase in flow velocity triggered the migration of all nitrogen species in the fracture, while an increase in matrix porosity reduced the concentration by enhancing the chemical reactions. An increase in fracture aperture also could trigger the denitrification process in the fracture to reduce the nitrate-nitrogen contamination in the fracture. The temperature variation between 25 °C and 45 °C in the fracture and the matrix essentially reduced the availability of nitrate-nitrogen and nitrogen gas in the fracture under non-isothermal conditions. Hence, an increase in the temperature coefficient can reduce the spike of nitrate-nitrogen and nitrogen gas in fracture by minimizing such transformation rates.