Long‐term assessment of nutrient flow pathway dynamics and in‐stream fate in a temperate karst agroecosystem watershed

Nutrient dynamics in karst agroecosystems remain poorly understood, in part due to limited long‐term nested datasets that can discriminate upland and in‐stream processes. We present a 10‐year dataset from a karst watershed in the Inner‐Bluegrass Region of central Kentucky, consisting of nitrate (nitrate‐N [NO₃^?]), dissolved reactive phosphorus (DRP), total organic carbon (TOC), and total ammoniacal‐N (TAN) measurements at nested spring and stream sites as well as flowrate at the watershed outlet. Hydrograph separation techniques were coupled with multiple linear regression and Empirical Mode Decomposition time‐series analysis to determine significance of seasonal processes and to generate continuous estimates of nutrient pathway loadings. Further, we used model results of benthic algae growth and decomposition dynamics from a nearby watershed to assess if transient storage in algal biomass could explain differences in spring and downstream watershed nutrient loading. Results highlight statistically significant seasonality for all nutrients at stream sites, but only for NO₃^? at springs with longitudinal variability showing significant decreases occurring from spring to stream sites for NO₃^? and DRP, and significant increases for TOC and TAN. Pathway loading analysis highlighted the importance of slow flow pathways to source approximately 70% of DRP and 80% of NO₃^?. Results for in‐stream dynamics suggest that benthic autotroph dynamics can explain summer deviations for TOC, TAN, and DRP but not NO₃^?. Regarding upland dynamics, our findings agree well with existing perceptions in karst for N pathways and upland source seasonality but deviate from perceptions that karst conduits are retentive of P, reflecting the limited buffering capacity of the soil profile and conduit sediments in the Inner‐Bluegrass. Regarding in‐stream fate, our findings highlighted the significance of seasonally driven nutrient processing in the bedrock‐controlled streambed to influence nutrient fluxes at the watershed outlet. Contrary to existing perceptions, we found high N attenuation and an unexplained NO₃^? sink in the bedrock stream, leading us to postulate that floating macrophytes facilitate high rates of denitrification.     

Assessing the potential of drone‐based thermal infrared imagery for quantifying river temperature heterogeneity

Climate change is altering river temperature regimes, modifying the dynamics of temperature‐sensitive fishes. The ability to map river temperature is therefore important for understanding the impacts of future warming. Thermal infrared (TIR) remote sensing has proven effective for river temperature mapping, but TIR surveys of rivers remain expensive. Recent drone‐based TIR systems present a potential solution to this problem. However, information regarding the utility of these miniaturised systems for surveying rivers is limited. Here, we present the results of several drone‐based TIR surveys conducted with a view to understanding their suitability for characterising river temperature heterogeneity. We find that drone‐based TIR data are able to clearly reveal the location and extent of discrete thermal inputs to rivers, but thermal imagery suffers from temperature drift‐induced bias, which prevents the extraction of accurate temperature data. Statistical analysis of the causes of this drift reveals that drone flight characteristics and environmental conditions at the time of acquisition explain ~66% of the variance in TIR sensor drift. These results shed important light on the factors influencing drone‐based TIR data quality and suggest that further technological development is required to enable the extraction of robust river temperature data. Nonetheless, this technology represents a promising approach for augmenting in situ sensor capabilities and improved quantification of advective inputs to rivers at intermediate spatial scales between point measurements and “conventional” airborne or satellite remote sensing.     

Water use strategies of two co‐occurring tree species in a semi‐arid karst environment

The flow of precipitation from the surface through to groundwater in karst systems is a complex process involving storage in the unsaturated zone and diffuse and preferential recharge pathways. The processes associated with this behaviour are not well understood, despite the prevalence of karst aquifers being used as freshwater supplies. As a result, uncertainty regarding the ecohydrological processes in this geological setting remains large. In response to the need to better understand the impact of woody vegetation on groundwater recharge, annual evapotranspiration (ET) rates and tree water sources were measured for two years above a shallow, fresh karst aquifer. Water use strategies of the co‐occurring Eucalyptus diversifolia subsp. diversifolia Bonpl. and Allocasuarina verticillata (Lam.) L. Johnson were investigated using a monthly water balance approach, in conjunction with measurement of the stable isotopes of water, leaf water potentials and soil matric potentials. The results suggest that it is unlikely groundwater resources are required to sustain tree transpiration, despite its shallow proximity to the soil surface, and that similarities exist between ET losses and the estimated long‐term average rainfall for this area. Irrespective of stand and morphological differences, E. diversifolia and A. verticillata ET rates showed remarkable convergence, demonstrating the ability of these co‐occurring species to maximise their use of the available precipitation, which avoids the requirement to differentiate between these species when estimating ET at a landscape scale. We conclude that the water holding capacity of porous geological substrates, such as those associated with karst systems, will play an important role in equilibrating annual rainfall variability and should be considered when assessing ecohydrological links associated with karst systems. Copyright © 2013 John Wiley & Sons, Ltd.     

Diffraction imaging of sub‐vertical fractures and karst with full‐resolution 3D Ground‐Penetrating Radar

Vertical fractures with openings of less than one centimetre and irregular karst cause abundant diffractions in Ground‐Penetrating Radar (GPR) records. GPR data acquired with half‐wavelength trace spacing are uninterpretable as they are dominated by spatially undersampled scattered energy. To evaluate the potential of high‐density 3D GPR diffraction imaging a 200 MHz survey with less than a quarter wavelength grid spacing (0.05 m × 0.1 m) was acquired at a fractured and karstified limestone quarry near the village of Cassis in Southern France. After 3D migration processing, diffraction apices line up in sub‐vertical fracture planes and cluster in locations of karstic dissolution features. The majority of karst is developed at intersections of two or more fractures and is limited in depth by a stratigraphic boundary. Such high‐resolution 3D GPR imaging offers an unprecedented internal view of a complex fractured carbonate reservoir model analogue. As seismic and GPR wave kinematics are similar, improvements in the imaging of steep fractures and irregular voids at the resolution limit can also be expected from high‐density seismic diffraction imaging.     

Aeolian fetch distance and secondary airflow effects: the influence of micro‐scale variables on meso‐scale foredune development

Unsuccessful attempts to use process‐scale models to predict long‐term aeolian sediment transport patterns have long been a feature of aeolian research. It has been proposed that one approach to overcome these problems is to identify micro‐scale variables that are important at longer timescales. This paper assesses the contribution of two system variables (secondary airflow patterns and fetch distance) to medium‐term (months to years) dune development. The micro‐scale importance of these variables had been established during previous work at the site (Magilligan Strand, Northern Ireland). Three methods were employed. First, sand drift potentials were calculated using 2 years of regional wind data and a sediment transport model. Second, wind data and large trench traps (2 m length × 1 m width × 1·5 m depth) were used to assess the actual sediment transport patterns over a 2‐month period. Third, a remote‐sensing technique for the identification of fetch distance, a saltation impact sensor (Safire) and wind data were utilized to gauge, qualitatively, sediment transport patterns over a 1‐month period. Secondary airflow effects were found to play a major role in the sediment flux patterns at these timescales, with measured and predicted rates matching closely during the trench trap study. The results suggest that fetch distance is an unimportant variable at this site. Copyright © 2007 John Wiley & Sons, Ltd.     

A review of low‐cost space‐borne data for flood modelling: topography,flood extent and water level

During the last two decades, remote sensing data have led to tremendous progress in advancing flood inundation modelling. In particular, low‐cost space‐borne data can be invaluable for large‐scale flood studies in data‐scarce areas. Various satellite products yield valuable information such as land surface elevation, flood extent and water level, which could potentially contribute to various flood studies. An increasing number of research studies have been dedicated to exploring those low‐cost data towards building, calibration and evaluation, and remote‐sensed information assimilation into hydraulic models. This paper aims at reviewing these recent scientific efforts on the integration of low‐cost space‐borne remote sensing data with flood modelling. Potentials and limitations of those data in flood modelling are discussed. This paper also introduces the future satellite missions and anticipates their likely impacts in flood modelling. Copyright © 2015 John Wiley & Sons, Ltd.     

Robotic photosieving from low‐cost multirotor sUAS: a proof‐of‐concept

Measurement of riverbed material grain sizes is now a routine part of fieldwork in fluvial geomorphology and lotic ecology. In the last decade, several authors have proposed remote sensing approaches of grain size measurements based on terrestrial and aerial imagery. Given the current rise of small unmanned aerial system (sUAS) applications in geomorphology, there is now increasing interest in the application of these remotely sensed grain size mapping methods to sUAS imagery. However, success in this area has been limited owing to two fundamental problems: lack of constraint of image scale for sUAS imagery and blurring effects in sUAS images and resulting orthomosaics. In this work, we solve the former by showing that SfM‐photogrammetry can be used in a direct georeferencing (DG) workflow (i.e. with no ground validation) in order to predict image scale within margins of 3%. We then propose a novel approach of robotic photosieving of dry exposed riverbed grains that relies on near‐ground images acquired from a low‐cost sUAS and which does not require the presence of ground control points or visible scale objects. We demonstrate that this absence of scale objects does not affect photosieving outputs thus resulting in a low‐cost and efficient sampling method for surficial grains. Copyright © 2017 John Wiley & Sons, Ltd.     

Improved Recharge Estimation from Portable,Low‐Cost Weather Stations

Groundwater recharge estimation is a critical quantity for sustainable groundwater management. The feasibility and robustness of recharge estimation was evaluated using physical‐based modeling procedures, and data from a low‐cost weather station with remote sensor techniques in Southern Abbotsford, British Columbia, Canada. Recharge was determined using the Richards‐based vadose zone hydrological model, HYDRUS‐1D. The required meteorological data were recorded with a HOBO^TM weather station for a short observation period (about 1 year) and an existing weather station (Abbotsford A) for long‐term study purpose (27 years). Undisturbed soil cores were taken at two locations in the vicinity of the HOBO^TM weather station. The derived soil hydraulic parameters were used to characterize the soil in the numerical model. Model performance was evaluated using observed soil moisture and soil temperature data obtained from subsurface remote sensors. A rigorous sensitivity analysis was used to test the robustness of the model. Recharge during the short observation period was estimated at 863 and 816 mm. The mean annual recharge was estimated at 848 and 859 mm/year based on a time series of 27 years. The relative ratio of annual recharge‐precipitation varied from 43% to 69%. From a monthly recharge perspective, the majority (80%) of recharge due to precipitation occurred during the hydrologic winter period. The comparison of the recharge estimates with other studies indicates a good agreement. Furthermore, this method is able to predict transient recharge estimates, and can provide a reasonable tool for estimates on nutrient leaching that is often controlled by strong precipitation events and rapid infiltration of water and nitrate into the soil.     

Suppression of near‐surface scattered body‐to‐surface waves: a steerable and nonlinear filtering approach

We present an approach based on local‐slope estimation for the separation of scattered surface waves from reflected body waves. The direct and scattered surface waves contain a significant amount of seismic energy. They present great challenges in land seismic data acquisition and processing, particularly in arid regions with complex near‐surface heterogeneities (e.g., dry river beds, wadis/large escarpments, and karst features). The near‐surface scattered body‐to‐surface waves, which have comparable amplitudes to reflections, can mask the seismic reflections. These difficulties, added to large amplitude direct and back‐scattered surface (Rayleigh) waves, create a major reduction in signal‐to‐noise ratio and degrade the final sub‐surface image quality. Removal of these waves can be difficult using conventional filtering methods, such as an filter, without distorting the reflected signal. The filtering algorithm we present is based on predicting the spatially varying slope of the noise, using steerable filters, and separating the signal and noise components by applying a directional nonlinear filter oriented toward the noise direction to predict the noise and then subtract it from the data. The slope estimation step using steerable filters is very efficient. It requires only a linear combination of a set of basis filters at fixed orientation to synthesize an image filtered at an arbitrary orientation. We apply our filtering approach to simulated data as well as to seismic data recorded in the field to suppress the scattered surface waves from reflected body waves, and we demonstrate its superiority over conventional techniques in signal preservation and noise suppression.     

Quantifying River‐Groundwater Interactions of New Zealand's Gravel‐Bed Rivers: The Wairau Plain

New Zealand's gravel‐bed rivers have deposited coarse, highly conductive gravel aquifers that are predominantly fed by river water. Managing their groundwater resources is challenging because the recharge mechanisms in these rivers are poorly understood and recharge rates are difficult to predict, particularly under a more variable future climate. To understand the river‐groundwater exchange processes in gravel‐bed rivers, we investigate the Wairau Plain Aquifer using a three‐dimensional groundwater flow model which was calibrated using targeted field observations, “soft” information from experts of the local water authority, parameter regularization techniques, and the model‐independent parameter estimation software PEST. The uncertainty of simulated river‐aquifer exchange flows, groundwater heads, spring flows, and mean transit times were evaluated using Null‐space Monte‐Carlo methods. Our analysis suggests that the river is hydraulically perched (losing) above the regional water table in its upper reaches and is gaining downstream where marine sediments overlay unconfined gravels. River recharge rates are on average 7.3 m³/s, but are highly dynamic in time and variable in space. Although the river discharge regularly hits 1000 m³/s, the net exchange flow rarely exceeds 12 m³/s and seems to be limited by the physical constraints of unit‐gradient flux under disconnected rivers. An important finding for the management of the aquifer is that changes in aquifer storage are mainly affected by the frequency and duration of low‐flow periods in the river. We hypothesize that the new insights into the river‐groundwater exchange mechanisms of the presented case study are transferable to other rivers with similar characteristics.     

Spectral‐temporal characterization of riverflow variability in England and Wales for the period 1865–2002

We have investigated riverflow variability in England and Wales by examining the reconstructed monthly discharge time series from fifteen catchments in these regions for the period 1865–2002. The riverflow fluctuations exhibit a strong annual cycle. The flow in the annual cycle is found to be intermittent, with the degree of intermittency varying from one catchment to another. An intermittent flow is characterized by bursts of high discharge separated by intervals with low or no discharge. By applying a continuous wavelet transform to the time series, we have identified the occurrence of intermittency in the annual cycle. The riverflow activity is also found to exhibit variations at interannual and quasi‐decadal time scales. These variations may be linked to large‐scale climatic processes such as the North Atlantic Oscillation (NAO). We have used the kurtosis of the probability density functions of the various time series as a measure of the degree of intermittency. An intermittent flow is characterized by a peaked (super‐Gaussian) probability density function with kurtosis in excess of 3. A higher value of kurtosis signifies a higher degree of intermittency. Intermittent fluctuations are more difficult to predict accurately than persistent oscillations, i.e., those lasting continuously over a long time interval. Copyright © 2009 John Wiley & Sons, Ltd.     

Assessing aeolian beach‐surface dynamics using a remote sensing approach

A remote sensing technique for assessing beach surface moisture was used to provide insight into beach‐surface evolution during an aeolian event. An experiment was carried out on 21 October 2007 at Greenwich Dunes, Prince Edward Island National Park, Canada, during which cameras were mounted on a mast on the foredune crest at a height of about 14 m above the beach. Maps of beach surface moisture were created based on a calibrated relationship between surface brightness from the photographs and surface moisture content measured in situ at points spaced every 2.5 m along a transect using a Delta‐T moisture probe. A time sequence of maps of surface moisture content captured beach surface evolution through the transport event at a spatial and temporal resolution that would be difficult to achieve with other sampling techniques such as impedance probes. Erosion of the foreshore and berm crest resulted in an increase in surface moisture content in these areas as the wetter underlying sediments were exposed. Flow expansion downwind of the berm crest led to sand deposition and a consequent decrease in surface moisture content. Remote sensing systems such as the one presented here allow observations of the combined evolution of beach surface moisture, shoreline position, and fetch distances during short‐term experiments and hence provide a comprehensive rendering of sediment erosion and transport processes. Copyright © 2012 John Wiley & Sons, Ltd.     

Soil creep and convex‐upward velocity profiles: theoretical and experimental investigation of disturbance‐driven sediment transport on hillslopes

The movement of unconsolidated materials near the Earth's surface is often driven by disturbances that occur at a range of spatial and temporal scales. The nature of these disturbances ranges from highly variable, such as tree turnover, to periodic and predictable, such as frost heave or creep. To explore the effect of probabilistic disturbances on surface processes, we formulated a granular creep model with analogy to rate process theory (RPT) used for chemical reactions. According to the theory, individual particles must be energized to a height greater than adjacent particles in order for grain dilation and transport to occur. The height of neighbouring particles (which is akin to activation energy in chemical reactions) varies with slope angle such that energy barriers get smaller in the downslope direction as slopes steepen. When slopes approach the friction‐limited angle of repose, the height of energy barriers approaches zero and grains ?ow in the absence of disturbance. An exponential function is used to describe the probability distribution of particle excitation height although alternative distributions are possible. We tested model predictions of granular dynamics in an experimental sandpile. In the sandpile, acoustic energy serves as the disturbance agent such that grains dilate and shear in response. Particle velocities are controlled by the frequency of energy pulses that result in grain displacement. Using tracer particles, we observed a convex‐upward velocity pro?le near the surface of the sandpile, consistent with predictions of our RPT‐based velocity model. In addition, we depth‐integrated the velocity model to predict how ?ux rates vary with inclination of the sandpile and observed non‐linear ?ux–gradient curves consistent with model predictions. By varying the acoustic energy level in the experimental sandpile, we documented changes in the rate of grain movement; similar changes in modelled velocities were achieved by varying the exponent of the particle excitation probability distribution. The general agreement between observed and modelled granular behaviour in our simple laboratory sandpile supports the utility of RPT‐based methods for modelling transport processes (e.g. soil creep, frost heave, and till deformation), thus enabling us to account for the probabilistic nature of disturbances that liberate sediment in natural landscapes. Copyright © 2004 John Wiley & Sons, Ltd.     

A multimethod approach to inform epikarst drip discharge modelling: Implications for palaeo‐climate reconstruction

Resolving the hydrological processes that form speleothems and the palaeo‐climate archives that they contain is difficult. Typical approaches to hydrological investigation are not suited to karst landscapes, geophysics are seldom applied, drip monitoring and modelling have limitations, and ignoring potential hydrological impacts can result in a proxy record that does not reflect the external environment. We aim to understand the processes and controls that have created a palaeo‐climate proxy record preserved in a speleothem (JC001) in the “Grotto of Oddities,” part of the Jersey Cave at the Yarrangobilly Caves, Australia, to infer the likely nature and resolution of this record. Electrical resistivity tomography (ERT), traditional surveying, and drip discharge monitoring (April 2013 to February 2015) were used to investigate the structure and hydrology of the epikarst overlying the Grotto of Oddities. Data collected through these methods were then used to construct a physically informed and parsimonious drip hydrology model. Geophysics showed that changes in hillslope above the Grotto of Oddities are collocated with a region of low resistivity, which forms an epikarstic reservoir acting to supply enhanced discharge to the speleothem. Drip monitoring showed hysteretic behaviour with a distinct threshold response, and a simple drip classification indicated that the speleothem associated with the drip has the potential to record palaeo‐seasonality or an annual–decadal signal. Discharge modelling indicated discharge was comprised of quick and slow flow, and that discharge is probably perennial. These multimethod results together indicate that the speleothem likely represents a palaeo‐climate record of a length and resolution unprecedented for nonglacial areas of the Southern Hemisphere and for Australia in particular and will significantly enhance current knowledge of the climate of southeast Australia. Although ERT methods have previously been applied in the karst landscape, to our knowledge, this represents the first application of these multiple methods in combination as an approach to assess the fidelity of a speleothem, based on an understanding of the hydrological processes for palaeo‐climate reconstruction.     

Water and air dynamics within a deep vadose zone of a karst massif: Observations from the Lukina jama–Trojama cave system (−1,431 m) in Dinaric karst (Croatia)

Extreme heterogeneity of karst systems makes them very challenging to study. Various processes within the system affect its global response, usually measured at karst springs. Research conducted in caves provides a unique opportunity for in situ analysis of separate processes in karst underground. The aim of the present study was to research the water and air dynamics within a deep karst system. Air and water basic physical parameters across the Lukina jama–Trojama cave system (?1,431 m) were continuously monitored during a 1‐year period. Recorded hydrograph of the siphon lake at the bottom of the cave was used to interpret the characteristics of an unexplored phreatic/epiphreatic conduit network. Water origin in the siphon was determined based on temperature and electrical conductivity. Air temperature and humidity monitoring revealed a strong inflow of air of sub‐zero temperature into the upper portion of the cave during winter. Cave passage morphology was interpreted as the main determinant of air dynamics, which caused ice to accumulate extensively in the upper portions of the cave and caused the temperature on the top of the homothermic zone to be significantly below the mean outside temperature. Air dynamics also lowered the temperature of water flowing through the cave vadose zone and feeding the phreatic zone of the massif. The pronounced temperature difference between the phreatic zone and the top of the homothermic zone probably contributed to the thermal gradient observed in the cave, which is steeper than in ice‐free caves in the area. Our results enabled the development of a conceptual model that describes coupling between air and water dynamics in the cave system and its surroundings.     

Effects of helical flow in one‐dimensional modelling of sediment distribution at river bifurcations

Complex flow processes at river bifurcations and the influence of the layout of a bifurcation make it difficult to predict sediment distribution over the downstream branches in case bedload transport dominates. In one‐dimensional models we need a nodal point relationship that prescribes the distribution of sediment over the downstream branches. We have identified which factors need to be included in such a relationship for the division of bedload transport at bifurcations. Next, irrotational flow theory for idealized geometries has been used to derive a simple physics‐based nodal point relationship that accounts for the effects of helical flow in the situation that a channel takes off under an angle from a straight main channel. This first step towards a complete nodal point relationship is applicable to bedload transport situations if the flow is clearly curved and if there is no pronounced bed topography. The relationship has been tested against data from a unique set of laboratory measurements, numerical data and data from a scale model of the Rhine bifurcation at Pannerden in the Netherlands. We find that the derived model yields a reasonable prediction of the sediment division over the downstream branches, and yields better predictions than the Wang et al. model for the situation considered. Considering the relative complexity and limited accuracy of the nodal point relationship for the effect of helical flow alone, however, we conclude thatderiving a practical physics‐based 1‐D relationship including all relevant processes is not feasible. We therefore recommend 2‐D or 3‐D modelling for all cases in general where morphological evolution depends on the division of bedload transport at bifurcations. Copyright © 2012 John Wiley & Sons, Ltd.     

Short‐term versus medium‐term monitoring for detecting gully‐erosion variability in a Mediterranean environment

This study investigates how medium‐term gully‐development data differ from short‐term data, and which factors influence their spatial and temporal variability at nine selected actively retreating bank gullies situated in four Spanish basin landscapes. Small‐format aerial photographs using unmanned, remote‐controlled platforms were taken at the gully sites in short‐term intervals of one to two years over medium‐term periods of seven to 13 years and gully change during each period was determined using stereophotogrammetry and a geographic information system. Results show a high variability of annual gully retreat rates both between gullies and between observation periods. The mean linear headcut retreat rates range between 0·02 and 0·26 m a^–1. Gully area loss was between 0·8 and 22 m² a^–1 and gully volume loss between 0·5 to 100 m³ a^–1, of which sidewall erosion may play a considerable part. A non‐linear relationship between catchment area and medium‐term gully headcut volume change was found for these gullies. The short‐term changes observed at the individual gullies show very high variability: on average, the maximum headcut volume change observed in 7–13 years was 14·3 times larger than the minimum change. Dependency on precipitation varies but is clearly higher for headcuts than sidewalls, especially in smaller and less disturbed catchments. The varying influences of land use and human activities with their positive or negative effects on runoff production and connectivity play a dominant role in these study areas, both for short‐term variability and medium‐term difference in gully development. The study proves the value of capturing spatially continuous, high‐resolution three‐dimensional data using small‐format aerial photography for detailed gully monitoring. Results confirm that short‐term data are not representative of longer‐term gully development and demonstrate the necessity for medium‐ to long‐term monitoring. However, short‐term data are still required to understand the processes – particularly human activity at varying time scales – causing fluctuations in gully erosion rates. Copyright © 2011 John Wiley & Sons, Ltd.     

Optical remote mapping of rivers at sub‐meter resolutions and watershed extents

At watershed extents, our understanding of river form, process and function is largely based on locally intensive mapping of river reaches, or on spatially extensive but low density data scattered throughout a watershed (e.g. cross sections). The net effect has been to characterize streams as discontinuous systems. Recent advances in optical remote sensing of rivers indicate that it should now be possible to generate accurate and continuous maps of in‐stream habitats, depths, algae, wood, stream power and other features at sub‐meter resolutions across entire watersheds so long as the water is clear and the aerial view is unobstructed. Such maps would transform river science and management by providing improved data, better models and explanation, and enhanced applications. Obstacles to achieving this vision include variations in optics associated with shadows, water clarity, variable substrates and target–sun angle geometry. Logistical obstacles are primarily due to the reliance of existing ground validation procedures on time‐of‐flight field measurements, which are impossible to accomplish at watershed extents, particularly in large and difficult to access river basins. Philosophical issues must also be addressed that relate to the expectations around accuracy assessment, the need for and utility of physically based models to evaluate remote sensing results and the ethics of revealing information about river resources at fine spatial resolutions. Despite these obstacles and issues, catchment extent remote river mapping is now feasible, as is demonstrated by a proof‐of‐concept example for the Nueces River, Texas, and examples of how different image types (radar, lidar, thermal) could be merged with optical imagery. The greatest obstacle to development and implementation of more remote sensing, catchment scale ‘river observatories’ is the absence of broadly based funding initiatives to support collaborative research by multiple investigators in different river settings. Copyright © 2007 John Wiley & Sons, Ltd.     

Prediction of groundwater‐induced flooding in a chalk aquifer for future climate change scenarios

Current climate change models for the southeast UK predict changing rainfall patterns, with increased incidence of extreme events. The chalk aquifer in the UK and northern France is susceptible to groundwater‐induced flooding under such conditions. In this methodological study we apply a frequency domain analysis approach to the chalk aquifer to derive a transfer function between effective rainfall and groundwater level from 7 years of monitoring data from the North Heath Barn site, near Brighton. The derived transfer function was calibrated and validated against monitoring data and then used to predict groundwater level for rainfall models for high, medium and low emission scenarios from the UKCP09 database. The derived transfer function is most closely comparable to the linear aquifer model, despite evidence for both matrix and fracture or karst water flow in the chalk, with transmissivity and unconfined storativity at the catchment scale of 1548 m² day^?1 and 1.6 × 10^?2. The application of the transfer function to UKCP09 rainfall data suggests that groundwater‐induced flooding may be about four times more frequent by 2040–2069 compared with 1961–1990 and seven times more frequent by 2070–2099. The model data also suggest an increase in the duration of groundwater minima relative to the reference period. Compared to deterministic modelling which requires detailed knowledge of aquifer heterogeneity and processes, the transfer function approach, although with limitations, is simpler, incorporating these factors into the analysis through frequency and phase coefficients, and thus may have the potential for groundwater risk assessment in other areas. Copyright © 2015 John Wiley & Sons, Ltd.     

Validating the KAGIS black‐box GIS‐based model in a Mediterranean karst aquifer: Case of study of Mela aquifer (SE Spain)

Karst Aquifer GIS‐based model (KAGIS model) is developed and applied to Mela aquifer, a small karst aquifer located in a Mediterranean region (SE Spain). This model considers different variables, such as precipitation, land use, surface slope and lithology, and their geographical heterogeneity to calculate both, the run‐off coefficients and the fraction of precipitation which contributes to fill the soil water reservoir existing above the aquifer. Evapotranspiration uptakes deplete water, exclusively, from this soil water reservoir and aquifer recharge occurs when water in the soil reservoir exceeds the soil field capacity. The proposed model also obtains variations of the effective porosity in a vertical profile, an intrinsic consequence of the karstification processes. A new proposal from the Nash–Sutcliffe efficiency index, adapted to arid environments, is presented and employed to evaluate the model's ability to predict the water table oscillations. The uncertainty in the model parameters is determined by the Generalized Likelihood Uncertainty Estimation method. Afterwards, when KAGIS is calibrated, wavelet analysis is applied to the resulting data in order to evaluate the variability in the aquifer behaviour. Wavelet analysis reveals that the rapid hydrogeological response, typical of a wide variety of karst systems, is the prevailing feature of Mela aquifer. This study proves that KAGIS is a useful tool to quantify recharge and discharge rates of karst aquifers and can be effectively applied to develop a proper management of water resources in Mediterranean areas.