Managed aquifer recharge (MAR) refers to a suite of methods that is increasingly being applied worldwide for sustainable groundwater management to tackle drinking or irrigation water shortage or to restore and maintain groundwater ecosystems. The potential for MAR is far from being exhausted, not only due to geological and hydrogeological conditions or technical and economic feasibility but also due to its lack of acceptance by the public and policymakers. One approach to enable the safe and accepted use of MAR could be to provide comprehensive risk management, including the identification, analysis, and evaluation of potential risks related to MAR. This article reviews current MAR risk assessment methodologies and guidelines and summarizes possible hazards and related processes. It may help planners and operators select the appropriate MAR risk assessment approaches and support the risk identification process. In addition to risk assessment (and subsequent risk treatment) related to the MAR implementation phase, this review also addresses risk assessment for MAR operation. We also highlight the limitations and lessons learned from the application and development of risk assessment methodologies. Moreover, developments are recommended in the area of MAR-related risk assessment methodologies and regulation. Depending on data availability, collected methodologies may be applicable for MAR sites worldwide. Integr Environ Assess Manag 2022;18:1513–1529. © 2022 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC). KEY POINTS The past two decades have seen a marked increase in the literature on MAR-related risk assessment, and a further increase is expected, since more attractive strategies are being used worldwide for integrated groundwater management. The development of holistic risk management plans can improve social acceptability and contribute to safer implementation and operation of MAR schemes. In connection with holistic risk management, MAR can be seen as a viable option for reaching the United Nations Sustainable Development Goal #6. Widely lacking regulations worldwide for risk assessment of MAR schemes is identified to hinder MAR implementation and safe operation worldwide.

Understanding transport and fate processes in the subsurface is of fundamental importance to identify the leaching potentials of herbicides or other compounds to groundwater resources. HYDRUS‐1D was used to simulate water flow and solute transport in arable land lysimeters. Simulations were compared to observed drainage rates and stable water isotopes (δ ¹⁸ O) in the drainage. Four different model setups were investigated and statistically evaluated for their model performance to identify dominant processes for water flow characterization in the vadose zone under similar cultivation management and climatic conditions. The studied lysimeters contain soil cores dominated by sandy gravel (Ly1) and clayey sandy silt (Ly2), both cropped with maize located in Wielenbach, Germany. First, a single‐porosity setup was chosen. For Ly1, modeling results were satisfactory, but for Ly2, the damping observed in the isotope signature of the drainage could not be fully covered. By considering immobile water with a dual‐porosity setup for Ly2, model performance improved. This could be due to a higher fraction of fine pores in Ly2 available for water storage, leading to mixing processes of isotopically enriched summer precipitation and lighter winter water. Accounting for isotopic evaporation fractionation processes in both model setups did not lead to improved model performance. Consequentially, the difference in soil hydraulic properties between the two lysimeters seems to impact water flow processes. Knowledge of such differences is crucial to prevent contamination and mitigate potential risks to soil and groundwater.

The catastrophic impact of Hurricane Dorian in September 2019 was unprecedented for the island of Grand Bahama. Flooding in the western portion of the island damaged pine ecosystems, inundated the soil and groundwater with salt water, and disrupted potable water service throughout the island. More than two years post-Hurricane Dorian, the freshwater lenses that the island relies on for potable water are still inundated with salt water. This collaborative paper summarizes all efforts of researchers and practitioners to evaluate the freshwater lenses, as well as their associated ecosystems, that serve as the main source of drinking water for the island of Grand Bahama. Hydrogeologic and vegetation assessments were conducted on the two primary wellfields that provide 95% of the drinking water to the island, over the span of two and a half years from the immediate aftermath of Hurricane Dorian through present day. While salinity and total dissolved solid concentrations in groundwater have declined, present levels indicate that the full recovery of the freshwater lenses may take decades. Forest assessments indicate that in Wellfield 6, which was the primary source of potable water pre-Hurricane Dorian, the pine forests suffered significant damage with complete pine mortality and little regeneration of pine trees occurring, which could impact the underlying freshwater lens. Lessons learned from this event underscore the vulnerability of water resources in The Bahamas and the critical need for adaptation strategies to improve resilience to future extreme events and climate change.

Characterizing unsaturated water flow in the subsurface is a requirement for understanding effects of droughts on agricultural production or impacts of climate change on groundwater recharge. By employing an improved lumped‐parameter model (LPM) approach that mimics variable flow we have interpreted stable water isotope data (δ ¹⁸ O and δ ² H), taken over 3 years at a lysimeter site located in Germany. Lysimeter soil cores were characterized by sandy gravel (Ly1) and clayey sandy silt (Ly2), and both lysimeters were vegetated with maize. Results were compared with numerical simulation of unsaturated flow and stable water isotope transport using HYDRUS‐1D. In addition, both approaches were extended by the consideration of preferential flow paths. Application of the extended LPM, and thus varying flow and transport parameters, substantially improved the description of stable water isotope observations in lysimeter seepage water. In general, findings obtained from the extended LPM were in good agreement to numerical modeling results. However, observations were more difficult to describe mathematically for Ly2, where the periodicity of seasonal stable water isotope fluctuation in seepage water was not fully met by numerical modeling. Furthermore, an extra isotopic upshift improved simulations for Ly2, probably controlled by stable water isotope exchange processes between mobile soil water and quasi‐immobile water within stagnant zones. Finally, although LPM requires less input data compared with numerical models, both approaches achieve comparable decision‐support integrity. The extended LPM approach can thus be a powerful tool for soil and groundwater management approaches.

This study investigates the contamination potential of herbicides to groundwater with the help of numerical modeling (HYDRUS‐1D) and stable carbon isotopes for characterizing biodegradation. Four herbicides, metolachlor, terbuthylazine, prosulfuron, and nicosulfuron, were applied over a period of 4.5 years on two lysimeters located in Wielenbach, Germany, and monitored by lysimeter drainage. These lysimeters contained soil cores dominated by sandy gravel (Ly1) and clayey sandy silt (Ly2) and were both cropped with maize ( Zea mays ). In the preceding study, we characterized flow within the lysimeters by using stable water isotopes and unsaturated flow models. Building up on these findings, models were extended for describing reactive transport of the herbicides and investigating process contributions. At the end of the experiment, 0.9%–15.9% of the applied herbicides (up to 20.9% if including metabolites) were recovered by lysimeter drainage. Metabolite formation and accumulation was observed, and biodegradation was also indicated by small changes in carbon isotope signals (δ ¹³ C) between applied and leached herbicides. Model setups could describe the dynamics of herbicide concentrations in lysimeter drainage well. Concentration peaks in drainage were partly also linked with strong precipitation events, indicating preferential flow influence. The soil core with the coarser texture (Ly1) showed less herbicide leaching than the finer texture (Ly2), which can be explained by a larger mobile phase in Ly1. Overall, our approaches and findings contribute to the understanding of multi‐process herbicide transport in the vadose zone and leaching potentials to groundwater, where δ ¹³ C can provide valuable hints for microbial degradation.