ALBERT

Description: We investigate how the choice of injection mode impacts transport properties in kilometer-scale three-dimensional discrete fracture networks (DFN). The choice of injection mode, resident or flux-weighted, is designed to mimic different physical phenomena. It has been hypothesized that solute plumes injected under resident conditions evolve to behave similarly to solutes injected under flux-weighted conditions. Previously, computational limitations have prohibited the large scale simulations required to investigate this hypothesis. We investigate this hypothesis by using a high performance DFN suite, dfnWorks , to simulate flow in kilometer-scale three-dimensional DFNs based on fractured granite at the Forsmark site in Sweden, and adopt a Lagrangian approach to simulate transport therein. Results show that after traveling through a pre-equilibrium region both injection methods exhibit linear scaling of the first moment of travel time and power law scaling of the breakthrough curve with similar exponents, slightly larger than two. The physical mechanisms behind this evolution appear to be the combination of in-network channeling of mass into larger fractures, which offer reduced resistance to flow, and in-fracture channeling, which results from the topology of the DFN. This article is protected by copyright. All rights reserved.

Description: Hydrologic ensemble forecasts driven by atmospheric ensemble prediction systems need statistical post-processing in order to account for systematic errors in terms of both location and spread. Runoff is an inherently multivariate process with typical events lasting from hours in case of floods to weeks or even months in case of droughts. This calls for multivariate post-processing techniques that yield well calibrated forecasts in univariate terms and ensure a realistic temporal dependence structure at the same time. To this end, the univariate ensemble model output statistics (EMOS) post-processing method is combined with two different copula approaches that ensure multivariate calibration throughout the entire forecast horizon. The domain of this study covers three sub-catchments of the river Rhine that represent different sizes and hydrological regimes: the Upper Rhine up to the gauge Maxau, the river Moselle up to the gauge Trier, and the river Lahn up to the gauge Kalkofen. In this study the two approaches to model the temporal dependence structure are ensemble copula coupling (ECC), which preserves the dependence structure of the raw ensemble, and a Gaussian copula approach (GCA), which estimates the temporal correlations from training observations. The results indicate that both methods are suitable for modelling the temporal dependencies of probabilistic hydrologic forecasts. This article is protected by copyright. All rights reserved.

Description: A three-dimensional mathematical model that describes transport of contaminant in a horizontal aquifer with simultaneous diffusion into a fractured clay formation is proposed. A group of semi-analytical solutions is derived based on specific initial and boundary conditions as well as various source functions. The analytical model solutions are evaluated by numerical Laplace inverse transformation and analytical Fourier inverse transformation. The model solutions can be used to study the fate and transport in a three-dimensional spatial domain in which a non-aqueous phase liquid exists as a pool atop a fractured low permeability clay layer. The non-aqueous phase liquid gradually dissolves into the groundwater flowing past the pool, while simultaneously diffusing into the fractured clay formation below the aquifer. Mass transfer of the contaminant into the clay formation is demonstrated to be significantly enhanced by the existence of the fractures, even though the volume of fractures is relatively small compared to the volume of the clay matrix. The model solution is a useful tool in assessing contaminant attenuation processes in a confined aquifer underlain by a fractured clay formation. This article is protected by copyright. All rights reserved.

Description: The objective of this study is to incorporate a time-dependent SCS CN method (SMA_CN) in Soil and Water Assessment Tool (SWAT) and compare its performance with the existing CN method in SWAT by simulating the hydrology of two agricultural watersheds in Indiana, United States. Results show that fusion of the SMA_CN method causes decrease in runoff volume and increase in profile soil moisture content, associated with larger groundwater contribution to the streamflow. In addition, the higher amount of moisture in the soil profile slightly elevates the actual evapotranspiration. The SMA-based SWAT configuration consistently produces improved goodness of fit scores and less uncertain outputs with respect to streamflow during both calibration and validation. The SMA_CN method exhibits better match with the observed data for all flow regimes, thereby addressing issues related to peak and low flow predictions by SWAT in many past studies. Comparison of the calibrated model outputs with field-scale soil moisture observations reveal that the SMA overhauling enables SWAT to represent soil moisture condition more accurately, with better response to the incident rainfall dynamics. While the results from the modification of the SCS method in SWAT are promising, more studies including watersheds with various physical and climatic settings are needed to validate the proposed approach. This article is protected by copyright. All rights reserved.

Description: We adapted Newton's Law of Cooling to model downstream water temperature change in response to stream-adjacent forest harvest on small and medium streams (average 327 ha in size) throughout the Oregon Coast Range, USA. The model requires measured stream gradient, width, depth and upstream control reach temperatures as inputs and contains two free parameters which were determined by fitting the model to measured stream temperature data. This model reproduces the measured downstream temperature responses to within 0.4 C ° for 15 of the 16 streams studied and provides insight into the physical sources of site-to-site variation among those responses. We also use the model to examine how the pre-to-post harvest change in daily maximum stream temperature depends on distance from the harvest reach. The model suggests that the pre-to-post harvest temperature change approximately 300  m downstream of the harvest will range from roughly 82% to less than 1% of that temperature change which occurred within the harvest reach, depending primarily on the downstream width, depth, and gradient. Using study-averaged values for these channel characteristics the model suggests that for a stream representative of those in the study, the temperature change approximately 300  m downstream of the harvest will be 56% of the temperature change which occurred within the harvest reach. This adapted Newton's Law of Cooling procedure represents a highly practical means for predicting stream temperature behavior downstream of timber harvests relative to conventional heat budget approaches, and is informative of the dominant processes affecting stream temperature. This article is protected by copyright. All rights reserved.

Description: Digital elevation models (DEMs) that are used in hydrological applications must be processed to remove sinks, mainly topographic depressions. Flow enforcement techniques include filling methods, which raise elevations within depressions, breaching, which carves channels through blockages, and hybrid methods. Despite previous research demonstrating the large impact to DEMs and subsequent analyses of depression filling, it is common practice apply this technique to flow enforcement. This is partly due to the greater efficiency of depression filling tools compared to breaching counterparts, which often limits breaching to applications of small- to moderate-sized DEMs. A new hybrid flow enforcement algorithm is presented in this study. The method can be run in complete breaching, selective breaching (either breached or filled), or constrained breaching (partial breaching) modes, allowing for greater flexibility in how practitioners enforce continuous flow paths. Algorithm performance was tested with DEMs of varying topography, spatial extents, and resolution. The sites included three moderate sized DEMs (52,000,000 to 190,000,000 cells) and three massive DEMs of the Iberian Peninsula, and the Amazon and Nile River basins, the largest containing nearly one billion cells. In complete breaching mode, the new algorithm required 87% of the time needed by a filling method to process the test DEMs, while the selective breaching and constrained breaching modes, operating with maximum breach depth constraints, increased run times by 8% and 27% respectively. Therefore, the new algorithm offers comparable performance to filling and the ability to process massive topographic data sets, while giving practitioners greater flexibility and lowering DEM impact. This article is protected by copyright. All rights reserved.

Description: Mountain snowpacks provide most of the annual discharge of western U.S. rivers, but the future of water resources in the western U.S. is tenuous, as climatic changes have resulted in earlier spring melts that have exacerbated summer droughts. Compounding changes to the physical environment are biotic disturbances including the mountain pine beetle (MPB), which has decimated millions of acres of western North American forests. At the watershed scale, MPB disturbance increases the peak hydrograph, and at the stand scale the ‘gray’ phase of MPB canopy disturbance decreases canopy snow interception, increases snow albedo, increases net shortwave radiation and decreases net longwave radiation versus the ‘red’ phase. Fewer studies have been conducted on the red phase of MPB disturbance, and in the mixed coniferous stands that may follow MPB-damaged forests. We measured the energy balance of four snowpacks representing different stages of MPB damage, management, and recovery: a lodgepole pine stand, a MPB-infested stand in the red phase, a mixed coniferous stand (representing one successional trajectory), and a clearcut (representing reactive management) in the Tenderfoot Creek Experimental Forest in Montana, USA. Net longwave radiation was lower in the MPB-infested stand despite higher basal area and plant area index of the other forests, suggesting that the dessicated needles serve as a less effective thermal buffer against longwave radiative losses. Eddy covariance observations of sensible and latent heat flux indicate that they are of similar but opposite magnitude, on the order of 20 MJ m −2 during the melt period. Further analyses reveal that net turbulent energy fluxes were near zero due to the temperature and atmospheric vapor pressure encountered during the melt period. Future research should place snow science in the context of forest succession and management, and address important uncertainties regarding the timing and magnitude of needlefall events. This article is protected by copyright. All rights reserved.

Description: This work focuses on the implementation of a Shallow Water-Exner model for compound natural channels with complex geometry and movable bed within the finite volume framework. The model is devised for compound channels modeling: cross-section overbanks are treated with fixed bed conditions, while the main channel is left free to modify its morphology. A capacitive approach is used for bedload transport modeling, in which the solid flow rates are estimated with bedload transport formulas. The model equations pose some numerical issues in the case of natural channels, where bedload transport may occur for both subcritical and supercritical flows and geometry varies in space. An explicit path-conservative scheme, designed to overcome all these issues, is presented in the paper. The scheme solves liquid and solid phases dynamics in a coupled manner, in order to correctly model near critical currents/channel interactions and is well-balanced, that is able to properly reproduce steady states. The Roe and Osher Riemann solvers are implemented, so as to take into account the spatial geometry variations of natural channels. The scheme reaches up to 2 nd order accuracy. Validation is performed with fixed and movable bed test cases whose analytical solution is known, and with flume experimental data. An application of the model to a real case study is also shown. This article is protected by copyright. All rights reserved.

Description: In this paper we use a physical modelling approach to explore the effect of lateral confinement on gravel bed river planform style, bed morphology, and sediment transport processes. A set of 27 runs was performed in a large flume (25 m long, 2.9 m wide), with constant longitudinal slope (0.01) and uniform grain size (1 mm), changing the water discharge (1.5 to 2.5 l/s) and the channel width (0.15 m to 1.5 m) to model a wide range of channel configurations, from narrow, straight, embanked channels to wide braided networks. The outcomes of each run were characterized by a detailed digital elevation model describing channel morphology, a map of dry areas and areas actively transporting sediment within the channel, and continuous monitoring of the amount of sediment transported through the flume outlet. Analysis reveals strong relationships between unit stream power and parameters describing the channel morphology. In particular, a smooth transition is observed between narrow channels with an almost rectangular cross section profile (with sediment transport occurring across the entire channel width) and complex braided networks where only a limited proportion (30%) of the bed is active. This transition is captured by descriptors of the bed elevation frequency distribution, e.g. standard deviation, skewness and kurtosis. These summary statistics represent potentially useful indicators of bed morphology that are compared with other commonly used summary indicators such as the braiding index and the type and number of bars. This article is protected by copyright. All rights reserved.