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.     

Lake District hillslopes

The Lake District is a region of great scenic beauty in north‐west England that has inspired artists and poets alike, and which comes high on the list of classic geological localities in Great Britain in terms of both bedrock and geomorphological features. With its inspiring views, the Lake District is often portrayed as the product of repeated glaciation, mainly because of the clarity of the erosional and depositional features that can be seen there. But since the last glaciers disappeared other processes have been modifying the landscape, processes that have superimposed their own signatures on to the glacial features. Hillslopes in particular have undergone significant changes, as a result of slope failures in both bedrock and superficial sediments. Although these landforms are not unknown, they have not received the same level of investigation as the glacial features, resulting in a limited appreciation of their spatial distribution and significance in reshaping the landscape. This article outlines the characteristics and origins of some slope failure types, and demonstrates that there is still much to learn about the Lake District landscape.     

Biostratigraphic and detrital zircon age constraints on the basement of the Himalayan Foreland Basin: Implications for a Proterozoic link to the Lesser Himalaya and cratonic India

The Himalayan Foreland Basin in the Ganga Valley is key to assessing the pre‐collision relationship between cratonic India and the Himalaya – the world's largest mountain chain. The subsurface Ganga Supergroup, representing the sedimentary basement of the Ganga Valley, has been interpreted as a northern extension of the Proterozoic Vindhyan Supergroup in cratonic India. This interpretation is contentious because the depositional age of the Ganga Supergroup is not resolved: whereas the lower Ganga Supergroup is widely regarded as Proterozoic, the upper Ganga Supergroup has been variously inferred to include Neoproterozoic, lower Palaeozoic, or Cretaceous strata. Here, we integrate biostratigraphic and detrital zircon data from drill cores to show that the entire Ganga Supergroup is likely Proterozoic and can be correlated with Proterozoic successions on the northern Indian craton and in the Lesser Himalaya. This helps redefine the first‐order stratigraphic architecture and indicates broad depositional continuity along the northern Indian margin during the Proterozoic.     

High-P (P = 1.5–1.8 GPa) blueschist from Elba: Implications for underthrusting and exhumation of continental units in the Northern Apennines

The Acquadolce Subunit on the Island of Elba, Italy, records blueschist facies metamorphism related to the Oligocene–early Miocene stages of continental collision in the Northern Apennines. The blueschist facies metamorphism is represented by glaucophane- and lawsonite-bearing metabasite associated with marble and calcschist. These rock types occur as lenses in a schistose complex representing foredeep deposits of early Oligocene age. Detailed petrological analyses on metabasic and metapelitic protoliths, involving mineral and bulk-rock chemistry coupled with P–T and P–T–X(Fe₂O₃) pseudosection modelling using PERPLE_X, show that the Acquadolce Subunit recorded nearly isothermal exhumation from peak pressure–temperature conditions of 1.5–1.8 GPa and 320–370°C. During exhumation, peak lawsonite- and possibly carpholite- or stilpnomelane-bearing assemblages were overprinted and partially obliterated by epidote-blueschist and, subsequently, albite-greenschist facies metamorphic assemblages. This study sheds new light on the tectonic evolution of Adria-derived metamorphic units in the Northern Apennines, by showing (a) the deep underthrusting of continental crust during continental collision and (b) rapid exhumation along ‘cold’ and nearly isothermal paths, compatible with syn-orogenic extrusion.     

An analytical approach to ascertain saturation‐excess versus infiltration‐excess overland flow in urban and reference landscapes

Uncontrolled overland flow drives flooding, erosion, and contaminant transport, with the severity of these outcomes often amplified in urban areas. In pervious media such as urban soils, overland flow is initiated via either infiltration‐excess (where precipitation rate exceeds infiltration capacity) or saturation‐excess (when precipitation volume exceeds soil profile storage) mechanisms. These processes call for different management strategies, making it important for municipalities to discern between them. In this study, we derived a generalized one‐dimensional model that distinguishes between infiltration‐excess overland flow (IEOF) and saturation‐excess overland flow (SEOF) using Green–Ampt infiltration concepts. Next, we applied this model to estimate overland flow generation from pervious areas in 11 U.S. cities. We used rainfall forcing that represented low‐ and high‐intensity events and compared responses among measured urban versus predevelopment reference soil hydraulic properties. The derivation showed that the propensity for IEOF versus SEOF is related to the equivalence between two nondimensional ratios: (a) precipitation rate to depth‐weighted hydraulic conductivity and (b) depth of soil profile restrictive layer to soil capillary potential. Across all cities, reference soil profiles were associated with greater IEOF for the high‐intensity set of storms, and urbanized soil profiles tended towards production of SEOF during the lower intensity set of storms. Urban soils produced more cumulative overland flow as a fraction of cumulative precipitation than did reference soils, particularly under conditions associated with SEOF. These results will assist cities in identifying the type and extent of interventions needed to manage storm water produced from pervious areas.     

Soil pipe flow tracer experiments: 1. Connectivity and transport characteristics

Much debate has occurred in catchment hydrology regarding the connectivity of flow paths from upslope areas to catchment outlets. This study was conducted in two catchments, one with three upper branches, in a loess soil with a fragipan that fosters lateral flow and exhibits an extensive distribution of soil pipe collapse features. The study aimed to determine the connectivity of multiple soil pipe networks as well as determine pipe flow velocities during storm events. Fluorescein dye was injected directly into soil pipes at the upper most pipe collapse feature of four different hillslopes. Breakthrough curves (BTC) were determined by sampling multiple pipe collapse features downslope. The BTCs were used to determine the ‘average’ (centre of mass) and ‘maximum’ (first arrival) flow velocities. This study confirmed that these catchments contain individual continuous soil pipe networks that extend over 190 m and connect the upper most hillslopes areas with the catchment outlet. While the flow paths are continuous, the individual pipe networks consist of alternating reaches of subsurface flow through soil pipes and reaches of surface flow through gullies formed by pipe collapses. In addition, flow can be occurring both through the subsurface soil pipes simultaneous with surface flow generated by artesian flow from the soil pipes. The pipe flow velocities were as high as 0.3 m/s, which was in the range of streamflow velocities. These pipe flow velocities were also in the range of velocities observed in pinhole erosion tests suggesting that these large, mature soil pipes are still actively eroding. Copyright © 2015 John Wiley & Sons, Ltd.     

Response of low‐angle dunes to variable flow

Current understanding of bedform dynamics is largely based on field and laboratory observations of bedforms in steady flow environments. There are relatively few investigations of bedforms in flows dominated by unsteadiness associated with rapidly changing flows or tides. As a consequence, the ability to predict bedform response to variable flow is rudimentary. Using high‐resolution multibeam bathymetric data, this study explores the dynamics of a dune field developed by tidally modulated, fluvially dominated flow in the Fraser River Estuary, British Columbia, Canada. The dunes were dominantly low lee angle features characteristic of large, deep river channels. Data were collected over a field ca 1·0 km long and 0·5 km wide through a complete diurnal tidal cycle during the rising limb of the hydrograph immediately prior to peak freshet, yielding the most comprehensive characterization of low‐angle dunes ever reported. The data show that bedform height and lee angle slope respond to variable flow by declining as the tide ebbs, then increasing as the tide rises and the flow velocities decrease. Bedform lengths do not appear to respond to the changes in velocity caused by the tides. Changes in the bedform height and lee angle have a counterclockwise hysteresis with mean flow velocity, indicating that changes in the bedform geometry lag changes in the flow. The data reveal that lee angle slope responds directly to suspended sediment concentration, supporting previous speculation that low‐angle dune morphology is maintained by erosion of the dune stoss and crest at high flow, and deposition of that material in the dune trough.     

Qualitative GIS: An Open Framework Using SpatiaLite and Open Source GIS

Qualitative GIS is a relatively new methodological approach for analyzing and visualizing qualitative data within a geographic context. Qualitative data can take many forms, including interviews, documents, photographs, and audio and video clips. Content analysis for example, is an effective qualitative method for analyzing text‐based data. We argue that basic concepts, (i.e. how to store data, data requirements, visualization techniques, and modes of analysis) within qualitative GIS have not been adequately defined, rendering difficult the replication of work performed and hindering the development of incremental knowledge in the field. Database management systems provide a means for storing, managing, and analyzing qualitative GIS data. A standardized and well‐designed open source database system provides a mechanism for qualitative GIS projects, ensuring consistency and project replication. Qualitative GIS data stored in a database allows for additional visualization options, such as geographic word clouds. To demonstrate the concepts we performed content analysis on Master Transportation Plans from Calgary and Montreal using SpatiaLite, an open source database system. We developed Structured Query Language (SQL) queries to generate and populate groups and theme tables within the SpatiaLite database. We present our database design and queries in the hopes that they will help others conducting qualitative GIS research.