Description: Subsurface flows, particularly hyporheic exchange fluxes, driven by streambed topography, permeability, channel gradient and dynamic flow conditions provide prominent ecological services such as nitrate removal from streams and aquifers. Stream flow dynamics cause strongly nonlinear and often episodic contributions of nutrient concentrations in river-aquifer systems. Using a fully coupled transient flow and reactive transport model, we investigated the denitrification potential of hyporheic zones during peak-flow events. The effects of streambed permeability, channel gradient and bedform amplitude on the spatio-temporal distribution of nitrate and dissolved organic carbon in streambeds and the associated denitrification potential were explored. Distinct peak-flow events with different intensity, duration and hydrograph shape were selected to represent a wide range of peak-flow scenarios. Our results indicated that the specific hydrodynamic characteristics of individual flow events largely determine the average positive or negative nitrate removal capacity of hyporheic zones, however the magnitude of this capacity is controlled by geomorphological settings (i.e. channel slope, streambed permeability and bedform amplitude). Specifically, events with longer duration and higher intensity were shown to promote higher nitrate removal efficiency with higher magnitude of removal efficiency in the scenarios with higher slope and permeability values. These results are essential for better assessment of the subsurface nitrate removal capacity under the influence of flow dynamics and particularly peak-flow events in order to provide tailored solutions for effective restoration of interconnected river-aquifer systems.

Description: Highlights • A new numerical model for permafrost in alpine regions. • Importance of lateral fluxes in mountain permafrost modeling. • Influence of unsaturated conditions on freezing processes. • Development of mountain permafrost during warming scenario. Abstract Alpine permafrost environments are highly vulnerable and sensitive to changes in regional and global climate trends. Thawing and degradation of permafrost has numerous adverse environmental, economic, and societal impacts. Mathematical modeling and numerical simulations provide powerful tools for predicting the degree of degradation and evolution of subsurface permafrost as a result of global warming. A particularly significant characteristic of alpine environments is the high variability in their surface geometry which drives large lateral thermal and fluid fluxes along topographic gradients. The combination of these topography-driven fluxes and unsaturated ground makes alpine systems markedly different from Arctic permafrost environments and general geotechnical ground freezing applications, and therefore, alpine permafrost demands its own specialized modeling approaches. In this work, we present a multi-physics permafrost model tailored to subsurface processes of alpine regions. In particular, we resolve the ice–water phase transitions, unsaturated conditions, and capillary actions, and account for the impact of the evolving pore space through freezing and thawing processes. Moreover, the approach is multi-dimensional, and therefore, inherently resolves the topography-driven horizontal fluxes. Through numerical case studies based on the elevation profiles of the Zugspitze (DE) and the Matterhorn (CH), we show the strong influence of lateral fluxes in 2D on active layer dynamics and the distribution of permafrost.