In recent years, a large body of the literature has been devoted to study reactive transport of solutes in porous media based on pure Lagrangian formulations. Such approaches have also been extended to accommodate second‐order bimolecular reactions, in which the reaction rate is proportional to the concentrations of the reactants. Rather, in some cases, chemical reactions involving two reactants follow more complicated rate laws. Some examples are (1) reaction rate laws written in terms of powers of concentrations, (2) redox reactions incorporating a limiting term (e.g., Michaelis‐Menten), or (3) any reaction where the activity coefficients vary with the concentration of the reactants, just to name a few. We provide a methodology to account for complex kinetic bimolecular reactions in a fully Lagrangian framework where each particle represents a fraction of the total mass of a specific solute. The method, built as an extension to the second‐order case, is based on the concept of optimal Kernel Density Estimator, which allows the concentrations to be written in terms of particle locations, hence transferring the concept of reaction rate to that of particle location distribution. By doing so, we can update the probability of particles reacting without the need to fully reconstruct the concentration maps. The performance and convergence of the method is tested for several illustrative examples that simulate the Advection‐Dispersion‐Reaction Equation in a 1‐D homogeneous column. Finally, a 2‐D application example is presented evaluating the need of fully describing non‐bilinear chemical kinetics in a randomly heterogeneous porous medium.

Improving degradation rates of emerging organic compounds (EOCs) in groundwater is still a challenge. Although their degradation is not fully understood, it has been observed that some substances are preferably degraded under specific redox conditions. The coupling of Managed Aquifer Recharge with soil aquifer remediation treatment, by placing a reactive layer containing organic matter at the bottom of the infiltration pond, is a promising technology to improve the rate of degradation of EOCs. Its success is based on assuming that recharged water and groundwater get well mixed, which is not always true. It has been demonstrated that mixing can be enhanced by inducing chaotic advection through extraction‐injection‐engineering. In this work, we analyze how chaotic advection might enhance the spreading of redox conditions with the final aim of improving degradation of a mix of benzotriazoles: benzotriazole, 5‐methyl‐benzotriazole, and 5‐chloro‐benzotriazole. The degradation of the first two compounds was fastest under aerobic conditions whereas the third compound was best degraded under denitrification conditions. We developed a reactive transport model that describes how a recharged water rich in organic matter mixes with groundwater, how this organic matter is oxidized by different electron acceptors, and how the benzotriazoles are degraded attending for the redox state. The model was tested in different scenarios of recharge, both in homogenous and in heterogenous media. It was found that chaotic flow increases the spreading of the plume of recharged water. Consequently, different redox conditions coexist at a given time, facilitating the degradation of EOCs.

Abstract. Managed aquifer recharge (MAR) can be affected by many risks. Those risks are related to different technical and non-technical aspects of recharge, like water availability, water quality, legislation, social issues, etc. Many other works have acknowledged risks of this nature theoretically; however, their quantification and definition has not been developed. In this study, the risk definition and quantification has been performed by means of "fault trees" and probabilistic risk assessment (PRA). We defined a fault tree with 65 basic events applicable to the operation phase. After that, we have applied this methodology to six different managed aquifer recharge sites located in the Mediterranean Basin (Portugal, Spain, Italy, Malta, and Israel). The probabilities of the basic events were defined by expert criteria, based on the knowledge of the different managers of the facilities. From that, we conclude that in all sites, the perception of the expert criteria of the non-technical aspects were as much or even more important than the technical aspects. Regarding the risk results, we observe that the total risk in three of the six sites was equal to or above 0.90. That would mean that the MAR facilities have a risk of failure equal to or higher than 90 % in the period of 2-6 years. The other three sites presented lower risks (75, 29, and 18 % for Malta, Menashe, and Serchio, respectively).

No abstract