Assessing regional‐scale spatio‐temporal patterns of groundwater–surface water interactions using a coupled SWAT‐MODFLOW model

Interaction between groundwater and surface water in watersheds has significant impacts on water management and water rights, nutrient loading from aquifers to streams, and in‐stream flow requirements for aquatic species. Of particular importance are the spatial patterns of these interactions. This study explores the spatio‐temporal patterns of groundwater discharge to a river system in a semi‐arid region, with methods applied to the Sprague River Watershed (4100 km²) within the Upper Klamath Basin in Oregon, USA. Patterns of groundwater–surface water interaction are explored throughout the watershed during the 1970–2003 time period using a coupled SWAT‐MODFLOW model tested against streamflow, groundwater level and field‐estimated reach‐specific groundwater discharge rates. Daily time steps and coupling are used, with groundwater discharge rates calculated for each model computational point along the stream. Model results also are averaged by month and by year to determine seasonal and decadal trends in groundwater discharge rates. Results show high spatial variability in groundwater discharge, with several locations showing no groundwater/surface water interaction. Average annual groundwater discharge is 20.5 m³/s, with maximum and minimum rates occurring in September–October and March–April, respectively. Annual average rates increase by approximately 0.02 m³/s per year over the 34‐year period, negligible compared with the average annual rate, although 70% of the stream network experiences an increase in groundwater discharge rate between 1970 and 2003. Results can assist with water management, identifying potential locations of heavy nutrient mass loading from the aquifer to streams and ecological assessment and planning focused on locations of high groundwater discharge. Copyright © 2016 John Wiley & Sons, Ltd.     

Effects of forest roads on the hydrological response of a small‐scale mountain watershed in Greece

The effects of land use changes on the ecology and hydrology of natural watersheds have long been debated. However, less attention has been given to the hydrological effects of forest roads. Although less studied, several researchers have claimed that streamflow changes related to forest roads can cause a persistent and pervasive effect on hillslope hydrology and the functioning of the channel system. The main potential direct effects of forest roads on natural watersheds hydrologic response are runoff production on roads surfaces due to reduced infiltration rates, interruption of subsurface flow by road cutslopes and rapid transfer of the produced runoff to the stream network through roadside ditches. The aforementioned effects may significantly modify the total volume and timing of the hillslope flow to the stream network. This study uses detailed field data, spatial data, hydro‐meteorological records, as well as numerical simulation to investigate the effects of forest roads on the hydrological response of a small‐scale mountain experimental watershed, which is situated in the east side of Penteli Mountain, Attica, Greece. The results of this study highlight the possible effects of forest roads on the watersheds hydrological response that may significantly influence direct runoff depths and peak flow rates. It is demonstrated that these effects can be very important in permeable watersheds and that more emphasis should be given on the impact of roads on the watersheds hydrological response. Copyright © 2014 John Wiley & Sons, Ltd.     

Piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events on the Loess Plateau,China

Establishing a universal watershed‐scale erosion and sediment yield prediction model represents a frontier field in erosion and soil/water conservation. The research presented here was conducted on the Chabagou watershed, which is located in the first sub‐region of the hill‐gully area of the Loess Plateau, China. A back‐propagation artificial neural model for watershed‐scale erosion and sediment yield was established, with the accuracy of the model, then compared with that of multiple linear regression. The sensitivity degree of various factors to erosion and sediment yield was quantitatively analysed using the default factor test. On the basis of the sensitive factors and the fractal information dimension, the piecewise prediction model for erosion and sediment yield of individual rainfall events was established and further verified. The results revealed the back‐propagation artificial neural network model to perform better than the multiple linear regression model in terms of predicting the erosion modulus, with the former able to effectively characterize dynamic changes in sediment yield under comprehensive factor conditions. The sensitivity of runoff erosion power and runoff depth to the erosion and sediment yield associated with individual rainfall events was found to be related to the complexity of surface topography. The characteristics of such a hydrological response are thus closely related to topography. When the fractal information dimension is greater than the topographic threshold, the accuracy of prediction using runoff erosion power is higher than that of using runoff depth. In contrast, when the fractal information dimension is smaller than the topographic threshold, the accuracy of prediction using runoff depth is higher than that of using runoff erosion power. The developed piecewise prediction model for watershed‐scale erosion and sediment yield of individual rainfall events, which introduces runoff erosion power and runoff depth using the fractal information dimension as a boundary, can be considered feasible and reliable and has a high prediction accuracy. Copyright © 2013 John Wiley & Sons, Ltd.     

The potential for agricultural land use change to reduce flood risk in a large watershed

Effects of agricultural land management practices on surface runoff are evident at local scales, but evidence for watershed‐scale impacts is limited. In this study, we used the Soil and Water Assessment Tool model to assess changes in downstream flood risks under different land uses for the large, intensely agricultural, Raccoon River watershed in Iowa. We first developed a baseline model for flood risk based on current land use and typical weather patterns and then simulated the effects of varying levels of increased perennials on the landscape under the same weather patterns. Results suggest that land use changes in the Raccoon River could reduce the likelihood of flood events, decreasing both the number of flood events and the frequency of severe floods. The duration of flood events were not substantially affected by land use change in our assessment. The greatest flood risk reduction was associated with converting all cropland to perennial vegetation, but we found that converting half of the land to perennial vegetation or extended rotations (and leaving the remaining area in cropland) could also have major effects on reducing downstream flooding potential. We discuss the potential costs of adopting the land use change in the watershed to illustrate the scale of subsidies required to induce large‐scale conversion to perennially based systems needed for flood risk reduction. Copyright © 2013 John Wiley & Sons, Ltd.     

Assessing the effects of land cover and future climate conditions on the provision of hydrological services in a medium‐sized watershed of Portugal

The separated and combined effects of land‐cover scenarios and future climate on the provision of hydrological services were evaluated in Vez watershed, northern Portugal. Soil and Water Assessment Tool was calibrated against daily discharge, sediments and nitrates, with good agreements between model predictions and field observations. Four hypothetical land‐cover scenarios were applied under current climate conditions (eucalyptus/pine, oak, agriculture/vine and low vegetation). A statistical downscaling of four General Circulation Models, bias‐corrected with ground observations, was carried out for 2021–2040 and 2041–2060, using representative concentration pathway 4.5 scenario. Also, the combined effects of future climate conditions were evaluated under eucalyptus/pine and agriculture/vine scenario. Results for land cover revealed that eucalyptus/pine scenario reduced by 7% the annual water quantity and up to 17% in the summer period. Although climate change has only a modest effect on the reduction of the total annual discharge (?7%), the effect on the water levels during summer was more pronounced, between ?15% and ?38%. This study shows that climate change can affect the provision of hydrological services by reducing dry season flows and by increasing flood risks during the wet months. Regarding the combined effects, future climate may reduce the low flows, which can be aggravated with eucalyptus/pine scenario. In turn, peak flows and soil erosion can be offset. Future climate may increase soil erosion and nitrate concentration, which can be aggravated with agriculture scenario. Results moreover emphasize the need to consider both climate and land‐cover impacts in adaptation and land management options at the watershed scale. Copyright © 2015 John Wiley & Sons, Ltd.     

Watershed delineation using hydrographic features and a DEM in plain river network region

Watershed delineation is a required step when conducting any spatially distributed hydrological modelling. Automated approaches are often proposed to delineate a watershed based on a river network extracted from the digital elevation model (DEM) using the deterministic eight‐neighbour (D8) method. However, a realistic river network cannot be derived from conventional DEM processing methods for a large flat area with a complex network of rivers, lakes, reservoirs, and polders, referred to as a plain river network region (PRNR). In this study, a new approach, which uses both hydrographic features and DEM, has been developed to address the problems of watershed delineation in PRNR. It extracts the river nodes and determines the flow directions of the river network based on a vector‐based hydrographic feature data model. The river network, lakes, reservoirs, and polders are then used to modify the flow directions of grid cells determined by D8 approach. The watershed is eventually delineated into four types of catchments including lakes, reservoirs, polders, and overland catchments based on the flow direction matrix and the location of river nodes. Multiple flow directions of grid cells are represented using a multi‐direction encoding method, and multiple outflows of catchments are also reflected in the topology of catchments. The proposed approach is applied to the western Taihu watershed in China. Comparisons between the results obtained from the D8 approach, the ‘stream burning’ approach, and those from the proposed approach clearly demonstrate an improvement of the new approach over the conventional approaches. This approach will benefit the development of distributed hydrological models in PRNR for the consideration of different types and multiple inlets and outlets of catchments. Copyright © 2015 John Wiley & Sons, Ltd.     

Oxygen‐18 dynamics in precipitation and streamflow in a semi‐arid agricultural watershed,Eastern Washington,USA

Understanding flow pathways and mechanisms that generate streamflow is important to understanding agrochemical contamination in surface waters in agricultural watersheds. Two environmental tracers, δ¹⁸O and electrical conductivity (EC), were monitored in tile drainage (draining 12 ha) and stream water (draining nested catchments of 6‐5700 ha) from 2000 to 2008 in the semi‐arid agricultural Missouri Flat Creek (MFC) watershed, near Pullman Washington, USA. Tile drainage and streamflow generated in the watershed were found to have baseline δ¹⁸O value of ?14·7‰ (VSMOW) year round. Winter precipitation accounted for 67% of total annual precipitation and was found to dominate streamflow, tile drainage, and groundwater recharge. ‘Old’ and ‘new’ water partitioning in streamflow were not identifiable using δ¹⁸O, but seasonal shifts of nitrate‐corrected EC suggest that deep soil pathways primarily generated summer streamflow (mean EC 250 µS/cm) while shallow soil pathways dominated streamflow generation during winter (EC declining as low as 100 µS/cm). Using summer isotopic and EC excursions from tile drainage in larger catchment (4700‐5700 ha) stream waters, summer in‐stream evaporation fractions were estimated to be from 20% to 40%, with the greatest evaporation occurring from August to October. Seasonal watershed and environmental tracer dynamics in the MFC watershed appeared to be similar to those at larger watershed scales in the Palouse River basin. A 0·9‰ enrichment, in shallow groundwater drained to streams (tile drainage and soil seepage), of δ¹⁸O values from 2000 to 2008 may be evidence of altered precipitation conditions due to the Pacific Decadal Oscillation (PDO) in the Inland Northwest. Copyright © 2009 John Wiley & Sons, Ltd.     

Granitoid emplacement by multiple sheeting during Variscan dextral transpression: The Saint-Laurent – La Jonquera pluton (Eastern Pyrenees)

The structural study of the Saint-Laurent – La Jonquera pluton (Eastern Pyrenees), a Variscan composite laccolithic intrusion emplaced in metasedimentary and gneissic rocks of the Roc de Frausa dome, by means of the anisotropy of magnetic susceptibility (AMS) technique has allowed the determination of the nature and orientation of its magmatic fabrics. The magmatic foliation has a predominant NE–SW strike and the mean lineation is also NE–SW trending with a shallow plunge. A strain gradient is measured so that the tonalites to granodiorites that form the basal parts of the pluton, and are intruded into amphibolite-facies metamorphic rocks, recorded the highest anisotropies, whereas the monzogranites and leucogranites, emplaced into upper crustal, low-grade metamorphic rocks, are weakly deformed. These results point to the synkinematic sequential emplacement of multiple granitoid sheets, from less to more differentiated magmatic stages, during the Late Carboniferous D₂ event characterized by an E–W-trending dextral transpression. The magmatic foliation appears locally disturbed by the effects of two tectonic events. The first of them (D₃) produced mylonitization of granitoids along NW–SE retrograding shear zones and open folds in the host Ediacaran metasediments of the Roc de Frausa massif, likely during late Variscan times. Interference between D₂ and D₃ structures was responsible for the dome geometry of the whole Roc de Frausa massif. The second and last perturbation consisted of local southward tilting of the granitoids coupled to the Mesozoic–Cenozoic cover during the Alpine.     

The diagenetic history of Oligocene-Miocene sandstones of the Austrian north Alpine foreland basin

Diagenesis is an essential tool to reconstruct the development of reservoir rocks. Diagenetic processes - precipitation and dissolution - have an influence on pore space. The present paper aims to study the diagenetic history of deep-marine sandstones of the Austrian Alpine Foreland Basin. To reach that goal, sediment petrology and diagenetic features of more than 110 sandstone samples from water- and gas-bearing sections from gas fields within the Oligocene-Miocene Puchkirchen Group and Hall Formation has been investigated. Special emphasis was put on samples in the vicinity of the gas-water contact (GWC). The sediment petrography of sandstones of Puchkirchen Group and Hall Formation is similar; hence their diagenesis proceeded the same way. In fact, primary mineralogy was controlled by paleo-geography with increasing transport distance and diverse detrital input.Sediment petrographically, investigated sandstones from the water-bearing horizon seemed quite comparable to the gas-bearing sediments. In general, they can be classified as feldspatic litharenites to litharenites and display porosities of up to 30% and permeabilities of up to 1300 mD. The carbon and oxygen isotopic composition of bulk carbonate cements from these sandstones range from−3.8 to +2.2 and from −7.5 to +0.2‰ [VPDB]. However, near the Gas-Water Contact (GWC) a horizon with low porosities (<3%) and permeabilities (<0.1 mD) is present. This zone is completely cemented with calcite, which has a blocky/homogenous morphology. A slight, but significant negative shift in δ18O isotopy (−2.5‰) is evident.During early diagenesis the first carbonate generations formed. First a fibrous calcite and afterwards a micritic calcite precipitated. Further siliciclastic minerals, such as quartz and feldspar (K-feldspar and minor plagioclase), exhibit corroded grains. Occasionally, clay minerals (illite; smectite, chlorite) formed as rims around detrital grains. Late diagenesis is indicated by the formation of a low permeable zone at the GWC.