“…More importantly, under the perspectives of health protection and pollution control, toxic metabolites contribute to dietary risk assessment [25] as well as surface water quality risk assessment [61]. Recent developments in coupling reaction networks with ecohydrological processes have improved the capability to predict herbicide dynamics in spite of the increased complexity in model structure (e.g., PRZM [17], MACRO [39], SWAT [4], HYDRUS [47], MODFLOW-RT3D [40], TOUGHREACT [87], and BRTSim [54]). These types of simulations need to be endowed with uncertainty and sensitivity analyses to account for errors in data collection, parameter value estimation, and model structure, which usually result in nonlinear model responses and unforeseen outcomes [20,70,76,88].…”