Search

Adolf and Günther Thiem, father and son, left behind a methodological legacy that many current hydrogeologists are probably unaware of. It goes far beyond the Dupuit–Thiem analytical model for pump test analysis, which is connected to their name. Methods, which we use on a day-to-day basis today, such as isopotential maps, tracer tests, and vertical wells, were amongst the many contributions which the Thiems either developed or improved. Remarkably, this was not done in a university context but rather as a by-product of their practical work of designing and building water supply schemes in countries all over Europe. Some of these waterworks are still active. Both Thiems were also great science communicators. Their contributions were read and applied in many countries, especially in the USA, through a personal connection between Günther Thiem and Oscar Edward Meinzer, the leading United States Geological Survey (USGS) hydrogeologist of the time.

Identifying and quantifying recharge processes linked to ephemeral surface water features is challenging due to their episodic nature. We use a combination of well-established near-surface geophysical methods to provide evidence of a surface and groundwater connection under a small ephemeral recharge feature in a flat, semi-arid region near Adelaide, Australia. We use a seismic survey to obtain P-wave velocity through travel-time tomography and S-wave velocity through the multichannel analysis of surface waves. The ratios between P-wave and S-wave velocities are used to calculate Poisson's ratio, which allow us to infer the position of the water table. Separate geophysical surveys were used to obtain electrical conductivity measurements from time-domain electromagnetics and water contents from downhole nuclear magnetic resonance. The geophysical observations provide evidence to support a groundwater mound underneath a subtle ephemeral surface water feature. Our results suggest that recharge is localized and that small-scale ephemeral features may play an important role in groundwater recharge. Furthermore, we show that a combined geophysical approach can provide a perspective that helps shape the hydrogeological conceptualization of a semi-arid region.

Study region: The ungauged, agriculturally dominated La Vi River Basin, Vietnam. Study focus: Groundwater abstraction for food and industrial production is increasing globally, putting pressure on groundwater resources and associated ecosystems. In many countries, monitoring of abstraction is poorly organised, resulting in a paucity of data, particularly in developing regions. Therefore, alternative approaches to estimate groundwater withdrawals are necessary. In this study, two soft-data approaches for indirect catchment-scale groundwater abstraction estimation are developed using: (1) local knowledge through a qualitative field survey of groundwater level fluctuations and groundwater withdrawals, and (2) land use data combined with local knowledge on cropping and irrigation practices. New hydrological insights for the region: The approaches are tested and applied for the La Vi River Basin, for the 2016 dry season. Total dry season estimated abstractions of 31.07 × 106 m3 and 36.19 × 106 m3 resulted from the two approaches. The advantage of the second approach is the spatial distribution of the estimated groundwater abstraction, aligning highly intensive abstractions with intensive agricultural areas. Despite high uncertainty in both estimates, this quantitative estimate gives valuable information for water managers, and the relatively good agreement between the methods provides trust in the estimates. The approaches are cost-effective and computationally simple solutions for estimating groundwater abstraction in data-poor regions.

In most hydrological systems, evapotranspiration (ET) and precipitation are the largest components of the water balance, which are difficult to estimate, particularly over complex terrain. In recent decades, the advent of remotely sensed data based ET algorithms and distributed hydrological models has provided improved spatially upscaled ET estimates. However, information on the performance of these methods at various spatial scales is limited. This study compares the ET from the MODIS remotely sensed ET dataset (MOD16) with the ET estimates from a SWAT hydrological model on graduated spatial scales for the complex terrain of the Sixth Creek Catchment of the Western Mount Lofty Ranges, South Australia. ET from both models was further compared with the coarser-resolution AWRA-L model at catchment scale. The SWAT model analyses are performed on daily timescales with a 6-year calibration period (2000–2005) and 7-year validation period (2007–2013). Differences in ET estimation between the SWAT and MOD16 methods of up to 31, 19, 15, 11 and 9 % were observed at respectively 1, 4, 9, 16 and 25 km2 spatial resolutions. Based on the results of the study, a spatial scale of confidence of 4 km2 for catchment-scale evapotranspiration is suggested in complex terrain. Land cover differences, HRU parameterisation in AWRA-L and catchment-scale averaging of input climate data in the SWAT semi-distributed model were identified as the principal sources of weaker correlations at higher spatial resolution.

Profiles of temperature time series are commonly used to determine hyporheic flow patterns and hydraulic dynamics in the streambed sediments. Although hyporheic flows are 3-D, past research has focused on determining the magnitude of the vertical flow component and how this varies spatially. This study used a portable 56-sensor, 3-D temperature array with three heat pulse sources to measure the flow direction and magnitude up to 200 mm below the water–sediment interface. Short, 1 min heat pulses were injected at one of the three heat sources and the temperature response was monitored over a period of 30 min. Breakthrough curves from each of the sensors were analysed using a heat transport equation. Parameter estimation and uncertainty analysis was undertaken using the differential evolution adaptive metropolis (DREAM) algorithm, an adaption of the Markov chain Monte Carlo method, to estimate the flux and its orientation. Measurements were conducted in the field and in a sand tank under an extensive range of controlled hydraulic conditions to validate the method. The use of short-duration heat pulses provided a rapid, accurate assessment technique for determining dynamic and multi-directional flow patterns in the hyporheic zone and is a basis for improved understanding of biogeochemical processes at the water–streambed interface.

Numerous basin aquifers in arid and semi-arid regions of the world derive a significant portion of their recharge from adjacent mountains. Such recharge can effectively occur through either stream infiltration in the mountain-front zone (mountain-front recharge, MFR) or subsurface flow from the mountain (mountain-block recharge, MBR). While a thorough understanding of recharge mechanisms is critical for conceptualizing and managing groundwater systems, distinguishing between MFR and MBR is difficult. We present an approach that uses hydraulic head, chloride and electrical conductivity (EC) data to distinguish between MFR and MBR. These variables are inexpensive to measure, and may be readily available from hydrogeological databases in many cases. Hydraulic heads can provide information on groundwater flow directions and stream–aquifer interactions, while chloride concentrations and EC values can be used to distinguish between different water sources if these have a distinct signature. Such information can provide evidence for the occurrence or absence of MFR and MBR. This approach is tested through application to the Adelaide Plains basin, South Australia. The recharge mechanisms of this basin have long been debated, in part due to difficulties in understanding the hydraulic role of faults. Both hydraulic head and chloride (equivalently, EC) data consistently suggest that streams are gaining in the adjacent Mount Lofty Ranges and losing when entering the basin. Moreover, the data indicate that not only the Quaternary aquifers but also the deeper Tertiary aquifers are recharged through MFR and not MBR. It is expected that this finding will have a significant impact on the management of water resources in the region. This study demonstrates the relevance of using hydraulic head, chloride and EC data to distinguish between MFR and MBR.

The normalised difference vegetation index (NDVI) is a useful tool for studying vegetation activity and ecosystem performance at a large spatial scale. In this study we use the Gravity Recovery and Climate Experiment (GRACE) total water storage (TWS) estimates to examine temporal variability of the NDVI across Australia. We aim to demonstrate a new method that reveals the moisture dependence of vegetation cover at different temporal resolutions. Time series of monthly GRACE TWS anomalies are decomposed into different temporal frequencies using a discrete wavelet transform and analysed against time series of the NDVI anomalies in a stepwise regression. The results show that combinations of different frequencies of decomposed GRACE TWS data explain NDVI temporal variations better than raw GRACE TWS alone. Generally, the NDVI appears to be more sensitive to interannual changes in water storage than shorter changes, though grassland-dominated areas are sensitive to higher-frequencies of water-storage changes. Different types of vegetation, defined by areas of land use type, show distinct differences in how they respond to the changes in water storage, which is generally consistent with our physical understanding. This unique method provides useful insight into how the NDVI is affected by changes in water storage at different temporal scales across land use types.

Study focus: Groundwater is of strategic importance. The accurate estimation of groundwater recharge and assessing the fundamental controlling factors are therefore of utmost importance to protect groundwater systems. We used the spatially-distributed water-balance model WetSpass to estimate long-term average recharge in Flanders. We validated recharge rates with base flow estimates of 67 daily stream flow records using the hydrograph analyses. To this end we performed principal component analysis, multiple linear regression analysis and relative importance analysis to assess the controlling factors of the spatial variation of recharge and base flow with the influencing watershed characteristics. New hydrological insights for the region: The average resulting recharge is 235 mm/year and occurs mainly in winter. The overall moderate correlation between base flow estimates and modeled recharge rates indicates that base flow is a reasonable proxy of recharge. Groundwater recharge variation was explained in order of importance by precipitation, soil texture and vegetation cover; while base flow variation was strongly controlled by vegetation cover and groundwater depth. The results of this study highlight the important role of spatial variables in estimation of recharge and base flow. In addition, the prominent role of vegetation makes clear the potential importance of land-use changes on recharge and hence the need to include a proper strategy for land-use change in sustainable management of groundwater resources.

As the availability of spatially distributed data sets for distributed rainfall-runoff modelling is strongly increasing, more attention should be paid to the influence of the quality of the data on the calibration. While a lot of progress has been made on using distributed data in simulations of hydrological models, sensitivity of spatial data with respect to model results is not well understood. In this paper we develop a spatial sensitivity analysis method for spatial input data (snow cover fraction – SCF) for a distributed rainfall-runoff model to investigate when the model is differently subjected to SCF uncertainty in different zones of the model. The analysis was focussed on the relation between the SCF sensitivity and the physical and spatial parameters and processes of a distributed rainfall-runoff model. The methodology is tested for the Biebrza River catchment, Poland, for which a distributed WetSpa model is set up to simulate 2 years of daily runoff. The sensitivity analysis uses the Latin-Hypercube One-factor-At-a-Time (LH-OAT) algorithm, which employs different response functions for each spatial parameter representing a 4 × 4 km snow zone. The results show that the spatial patterns of sensitivity can be easily interpreted by co-occurrence of different environmental factors such as geomorphology, soil texture, land use, precipitation and temperature. Moreover, the spatial pattern of sensitivity under different response functions is related to different spatial parameters and physical processes. The results clearly show that the LH-OAT algorithm is suitable for our spatial sensitivity analysis approach and that the SCF is spatially sensitive in the WetSpa model. The developed method can be easily applied to other models and other spatial data.

In Europe, public water withdrawals make up on average 30% and in some cases up to 60% of total water withdrawals. These withdrawals are becoming increasingly important with growing population density; hence there is a need to understand the spatial and temporal trends involved. Pan-European public/municipal water withdrawals and consumption were mapped for 2006 and forecasted for 2030. Population and tourism density were assumed to be the main driving factors for withdrawals. Country-level statistics on public water withdrawals were disaggregated to a combined population and tourism density map (the "user" density map) computed for 2006. The methodology was validated using actual regional withdrawal statistics from France for 2006. The total absolute error (TAE) calculated was proven to be reduced by taking into account the tourism density in addition to the population density. In order to forecast the map to 2030 we considered a reference scenario where per capita withdrawals were kept constant in time. Although there are large variations from region to region, this resulted in a European average increase of water withdrawals of 16%. If we extrapolate the average reduction in per capita withdrawals seen between 2000 and 2008, we forecast a reduction in average total water withdrawals of 4%. Considering a scenario where all countries converge to an optimal water use efficiency, we see an average decrease of 28%.

Outcropping sediments can be used as easily accessible analogues for studying subsurface sediments, especially to determine the small-scale spatial variability of hydrogeological parameters. The use of cost-effective in situ measurement techniques potentially makes the study of outcrop sediments even more attractive. We investigate to what degree air-permeameter measurements on outcrops of unconsolidated sediments can be a proxy for aquifer saturated hydraulic conductivity (K) heterogeneity. The Neogene aquifer in northern Belgium, known as a major groundwater resource, is used as the case study. K and grain-size data obtained from different outcropping sediments are compared with K and grain-size data from aquifer sediments obtained either via laboratory analyses on undisturbed borehole cores (K and grain size) or via large-scale pumping tests (K only). This comparison shows a pronounced and systematic difference between outcrop and aquifer sediments. Part of this difference is attributed to grain-size variations and earth surface processes specific to outcrop environments, including root growth, bioturbation, and weathering. Moreover, palaeoenvironmental conditions such as freezing–drying cycles and differential compaction histories will further alter the initial hydrogeological properties of the outcrop sediments. A linear correction is developed for rescaling the outcrop data to the subsurface data. The spatial structure pertaining to outcrops complements that obtained from the borehole cores in several cases. The higher spatial resolution of the outcrop measurements identifies small-scale spatial structures that remain undetected in the lower resolution borehole data. Insights in stratigraphic and K heterogeneity obtained from outcrop sediments improve developing conceptual models of groundwater flow and transport.

We present a regional framework for an integrated and spatiotemporally distributed assessment of human-induced trends in the hydrology and the associated ecological health of a semi-arid basin where both human activities (i.e. agriculture) and natural ecosystems are highly groundwater dependent. To achieve this, we analysed the recent trends (from year 2000 to 2010) in precipitation, evapotranspiration (actual and potential) and vegetation greenness (i.e. NDVI) using a combination of satellite and ground-based observations. The trend assessment was applied for the semi-arid Konya Basin (Turkey), one of the largest endorheic basins in the world. The results revealed a consistent increasing trend of both yearly evapotranspiration (totally 63 MCM yr−1 from croplands) and mean NDVI (about 0.004 NDVI yr−1 in irrigated croplands), especially concentrating in the plain part of the basin, while no significant trends were observed for the precipitation and potential evapotranspiration variables. On the contrary, a consistent decreasing trend of both yearly evapotranspiration (totally −2.1 MCM yr−1) and mean NDVI (−0.001 NDVI yr−1) was observed in the wetlands, which also cannot be explained by trends in precipitation and potential evapotranspiration. The emerging picture suggest that the greening trend of the vegetation and increasing of evapotranspiration in the plain are related to land cover changes (i.e. conversion into irrigated croplands) and to the intensification of the supplementary irrigation for agriculture, which in turn caused drying out of some wetlands and the natural vegetation which mostly depend on the groundwater, the main source of irrigation water as well. Our study presented an example of the utility of spatially and temporally continuous RS data in assessing the regional trends in hydrological and ecological variables and their interactions in a spatially distributed manner in a semi-arid region, which can also be adapted to other regions. Such spatiotemporally distributed analysis at the basin level is particularly important considering that most of the water management interventions also take place at this scale.

For policy making and spatial planning, information is needed about the impacts of climate change on natural ecosystems. To provide this information, commonly hydrological and ecological models are used. We give arguments for our assessment that modelling only is insufficient for determining the impacts of climate changes on natural ecosystems at regional scales. Instead, we proposed a combination of hydrological simulations, a literature review and process-knowledge on climate-hydrology-vegetation interactions, to compile a sketch map that indicates climate change effects on a number of ecosystems in the Netherlands. Soon after a first version of our sketch map was published by a Dutch professional journal, copies appeared in policy documents, and also in a commercial and popular atlas of the Netherlands. Moreover, the map led to a question in the Dutch parliament about the sustainability of bog reserves under the future climate. Apparently, there was an urgent need for the information provided by the map. The map shows that climate change will presumably have the largest influence on ecosystems in the Netherlands that depend on precipitation as the major water source, like heathlands, dry grasslands, rain-fed moorland pools and raised bogs. Also highly susceptible are fens in reserves surrounded by deeply drained polders, because such fens depend on the inlet of surface water, of which quality is likely to deteriorate upon climate change. While the map is indicative for directions of change, in view of the uncertainties of our study, no conclusions should be drawn that may have far-reaching consequences, such as giving up certain nature targets that might no longer be feasible in the future climate. Instead, we advise to anticipate the potential threats from climate change by taking a number of adaptation measures that enhance the robustness of nature reserves. To improve climate change projections on hydrology and ecosystems, future research should especially focus on feedbacks of vegetation on the water balance, on processes that directly influence plant performance and on the ecological effects of weather extremes.

Monitoring evapotranspiration over land is highly dependent on the surface state and vegetation dynamics. Data from spaceborn platforms are desirable to complement estimations from land surface models. The success of daily evapotranspiration monitoring at continental scale relies on the availability, quality and continuity of such data. The biophysical variables derived from SEVIRI on board the geostationary satellite Meteosat Second Generation (MSG) and distributed by the Satellite Application Facility on Land surface Analysis (LSA-SAF) are particularly interesting for such applications, as they aimed at providing continuous and consistent daily time series in near-real time over Africa, Europe and South America. In this paper, we compare them to monthly vegetation parameters from a database commonly used in numerical weather predictions (ECOCLIMAP-I), showing the benefits of the new daily products in detecting the spatial and temporal (seasonal and inter-annual) variability of the vegetation, especially relevant over Africa. We propose a method to handle Leaf Area Index (LAI) and Fractional Vegetation Cover (FVC) products for evapotranspiration monitoring with a land surface model at 3–5 km spatial resolution. The method is conceived to be applicable for near-real time processes at continental scale and relies on the use of a land cover map. We assess the impact of using LSA-SAF biophysical variables compared to ECOCLIMAP-I on evapotranspiration estimated by the land surface model H-TESSEL. Comparison with in-situ observations in Europe and Africa shows an improved estimation of the evapotranspiration, especially in semi-arid climates. Finally, the impact on the land surface modelled evapotranspiration is compared over a north–south transect with a large gradient of vegetation and climate in Western Africa using LSA-SAF radiation forcing derived from remote sensing. Differences are highlighted. An evaluation against remote sensing derived land surface temperature shows an improvement of the evapotranspiration simulations.

As recognized in the European Water Framework Directive, groundwater-dependent wetlands and their diverse ecosystems have important functions which need to be protected. The vegetation in such habitats is often dependent on quality, quantity and patterns of river discharge and groundwater-surface water interaction on a local or reach scale. Since groundwater-surface water exchange studies on natural rivers and wetlands with organic soils are scarce, more functional analysis is needed. To this end we combined different field methods including piezometer nests, temperature as tracer and seepage meter measurements. Some of these measurements were used as inputs and/or as validation for the numerical 1-D heat transport model STRIVE. In transient mode the model was used to calculate spatially distributed vertical exchange fluxes from temperature profiles measured at the upper Biebrza River in Poland over a period of nine months. Time series of estimated fluxes and hydraulic head gradients in the hyporheic zone were used to estimate the temporal variability of groundwater-surface water exchange. This paper presents a hierarchical approach for quantifying and interpreting groundwater-surface water interaction in space and time. The results for the upper Biebrza show predominantly upward water fluxes, sections of recharge, however, exist along the reach. The fluxes depend more on hydraulic gradients than on riverbed conductivity. This indicates that the fluvio-plain scale is required for interpreting the exchange fluxes, which are estimated on a local scale. The paper shows that a conceptual framework is necessary for understanding the groundwater-surface water interaction processes, where the exchange fluxes are influenced by local factors like the composition of the riverbed and the position of the measurement on a local scale, and by regional factors like the hydrogeology and topography on a fluvio-plain scale. The hierarchical methodology increases the confidence in the estimated exchange fluxes and improves the process understanding. The accuracy of the measurements and related uncertainties, however, remain challenges for wetland environments. Gaining quantitative information on groundwater-surface water interaction can improve modeling confidence and as a consequence helps to develop effective procedures for management and conservation of valuable groundwater dependent wetlands.

Given the importance of groundwater for food production and drinking water supply, but also for the survival of groundwater dependent terrestrial ecosystems (GWDTEs) it is essential to assess the impact of climate change on this freshwater resource. In this paper we study with high temporal and spatial resolution the impact of 28 climate change scenarios on the groundwater system of a lowland catchment in Belgium. Our results show for the scenario period 2070–2101 compared with the reference period 1960–1991, a change in annual groundwater recharge between −20% and +7%. On average annual groundwater recharge decreases 7%. In most scenarios the recharge increases during winter but decreases during summer. The altered recharge patterns cause the groundwater level to decrease significantly from September to January. On average the groundwater level decreases about 7 cm with a standard deviation between the scenarios of 5 cm. Groundwater levels in interfluves and upstream areas are more sensitive to climate change than groundwater levels in the river valley. Groundwater discharge to GWDTEs is expected to decrease during late summer and autumn as much as 10%, though the discharge remains at reference-period level during winter and early spring. As GWDTEs are strongly influenced by temporal dynamics of the groundwater system, close monitoring of groundwater and implementation of adaptive management measures are required to prevent ecological loss.

In this work we assess the uncertainty in modelling the groundwater flow for the Pampa del Tamarugal Aquifer (PTA) – North Chile using a novel and fully integrated multi-model approach aimed at explicitly accounting for uncertainties arising from the definition of alternative conceptual models. The approach integrates the Generalized Likelihood Uncertainty Estimation (GLUE) and Bayesian Model Averaging (BMA) methods. For each member of an ensemble M of potential conceptualizations, model weights used in BMA for multi-model aggregation are obtained from GLUE-based likelihood values. These model weights are based on model performance, thus, reflecting how well a conceptualization reproduces an observed dataset D. GLUE-based cumulative predictive distributions for each member of M are then aggregated obtaining predictive distributions accounting for conceptual model uncertainties. For the PTA we propose an ensemble of eight alternative conceptualizations covering all major features of groundwater flow models independently developed in past studies and including two recharge mechanisms which have been source of debate for several years. Results showed that accounting for heterogeneities in the hydraulic conductivity field (a) reduced the uncertainty in the estimations of parameters and state variables, and (b) increased the corresponding model weights used for multi-model aggregation. This was more noticeable when the hydraulic conductivity field was conditioned on available hydraulic conductivity measurements. Contribution of conceptual model uncertainty to the predictive uncertainty varied between 6% and 64% for ground water head estimations and between 16% and 79% for ground water flow estimations. These results clearly illustrate the relevance of conceptual model uncertainty.

Land-use changes are frequently indicated to be one of the main human-induced factors influencing the groundwater system. For land-use change, groundwater research has mainly focused on the change in water quality thereby neglecting changes in quantity. The objective of this paper is to assess the impact of land-use changes, from 2000 until 2020, on the hydrological balance and in particular on groundwater quantity, as results from a case study in the Kleine Nete basin, Belgium. New is that this study tests a methodology, which couples a land-use change model with a water balance and a steady-state groundwater model. Four future land-use scenarios (A1, A2, B1 and B2) based on the Special Report on Emission Scenarios (SRES) are modelled with the CLUE-S model. Water balance components, groundwater level and baseflow are simulated using the WetSpass model in conjunction with a steady-state MODFLOW groundwater flow model. Results show that the average recharge decreases with 2.9, 1.6, 1.8 and 0.8% for scenario A1, A2, B1 and B2, respectively, over the 20 covered years. The predicted reduction in recharge results in a small decrease of the average groundwater level in the basin, ranging from 2.5 cm for scenario A1 to 0.9 cm for scenario B2, and a reduction of the baseflow with maximum 2.3% and minimum 0.7% for scenario A1 and B2, respectively. Although these averages appear to indicate small changes in the groundwater system, spatial analysis shows that much larger changes are located near the major cities in the study area. Hence, spatial planning should take better account of effects of land-use change on the groundwater system and define mitigating actions for reducing the negative impacts of land-use change.

Abstract  The effects of climate change on the groundwater systems in the Grote-Nete catchment, Belgium, covering an area of 525 km2, is modeled using wet (greenhouse), cold or NATCC (North Atlantic Thermohaline Circulation Change) and dry climate scenarios. Low, central and high estimates of temperature changes are adopted for wet scenarios. Seasonal and annual water balance components including groundwater recharge are simulated using the WetSpass model, while mean annual groundwater elevations and discharge are simulated with a steady-state MODFLOW groundwater model. WetSpass results for the wet scenarios show that wet winters and drier summers are expected relative to the present situation. MODFLOW results for wet high scenario show groundwater levels increase by as much as 79 cm, which could affect the distribution and species richness of meadows. Results obtained for cold scenarios depict drier winters and wetter summers relative to the present. The dry scenarios predict dry conditions for the whole year. There is no recharge during the summer, which is mainly attributed to high evapotranspiration rates by forests and low precipitation. Average annual groundwater levels drop by 0.5 m, with maximum of 3.1 m on the eastern part of the Campine Plateau. This could endanger aquatic ecosystem, shrubs, and crop production. [ABSTRACT FROM AUTHOR]

Este artículo presenta el resultado de una investigación de conductos subterráneos, los cuales conectan dolinas y la resurgencia de un río ciego en la cuenca tropical Nam La, compuesta por calizas altamente karstificadas, en el noroeste de Vietnam. El Río Nam La desaparece en el subsuelo por medio de varias dolinas cerca de la salida de la cuenca. Durante la estación lluviosa esto resulta en inundaciones aguas arriba de las dolinas. Una hipótesis es que el Río Nam La resurge en un manantial de caverna grande a 4.5 km al oriente de la salida de la cuenca. Se llevó a cabo un estudio multi-temático de las conexiones posibles entre las dolinas y la resurgencia para investigar la estructura geológica, tectónica, análisis estructural de cavernas, y series de tiempo de descarga. Se evalúa y se demuestra también la existencia de conductos subterráneos por medio de experimentos con trazadores. En base a análisis de lineamientos se predice la localización de los conductos subterráneos. Se utiliza un mapa de densidad de longitudes de lineamientos construido a partir de sensores remotos para seguir las rutas de las dolinas a la resurgencia, teniendo longitudes más cortas pero densidad de lineamientos más altas. Esto da por resultado una predicción en vista de planta de conductos subterráneos que se ajustan con el desarrollo de cavernas y fallas. El funcionamiento de los conductos se explica posteriormente mediante el análisis de registro de inundaciones de una dolina cercana la cual actúa como un reservorio de almacenamiento temporal que mitiga inundaciones en el área de salida de la cuenca. [ABSTRACT FROM AUTHOR]

For a karstified limestone area in NW Vietnam, the relationship between the distribution of lineaments and borehole specific capacity is determined, resulting in the conclusion that not only the borehole geomorphological-hydrogeological position but also the lineament distribution influences the specific capacity.No significant spatial well yield patterns are evident in this highly fractured-karstified region. The supposition is that lineaments caused by geotectonic activities affect the local variability in borehole specific capacity. Sixteen pumping tests in conjunction with a comprehensive lineament analysis are used to prove this relationship. The boreholes and lineaments are classified into two groups according to their similarity in geomorphological-hydrogeological features. Lineaments tend to be less detectable in discharge areas (lowland, wide and flat valleys) in contrast to the high density in recharge areas (highland narrow-mountainous ravines). In addition, the presence of a stream network in the former can act as a recharge source to the underlain karstic groundwater system. Consequently, boreholes that are in the discharge areas with a lower density of lineaments often produce high yield. For recharge areas with a high density of lineaments, a good correlation is found between lineament density and borehole specific capacity. [ABSTRACT FROM AUTHOR]

The dune ecosystem plays a significant role in the global carbon cycle. The Horqin Sandy Land is a typical semi-arid fragile ecosystem in northern China. Understanding the magnitudes and dynamics of carbon dioxide fluxes within this region is essential for understanding the carbon balance. Used 6 years (2013-2018) measurements from an eddy-covariance system, we analyzed the dynamic patterns of net ecosystem carbon exchange (NEE), gross primary production (GPP), and ecosystem respiration (Reco) of the dune ecosystem in Horqin Sandy Land and examined their responses to climate factors with a focus on the precipitation. The results showed that the NEE of the dune ecosystem fluctuated from -166 to 100 gCO 2 ·m -2 ·year -1 across the 6 growing seasons, with an average of -56 gCO 2 ·m -2 ·year -1 . The precipitation was not a key factor influencing the carbon flux variability. During the mid-growth stage, GPP was primarily affected by the effective precipitation frequency (R 2 ranging from 0.65 to 0.85, P < 0.05), followed by fractional vegetation cover (R 2 ranging from 0.65 to 0.68, P < 0.05). However, in the early and late growth stages, temperature predominantly drove the carbon flux (R 2  = 0.75, P < 0.01). The interannual variability of carbon flux can be predominantly elucidated by phenological indicators such as CO 2 uptake (CU start ), end of CO 2 uptake (CU end ), CO 2 uptake period (CUP), and Spring lag. The results demonstrated the dune ecosystem is a weak carbon sink in semi-arid ecosystems. Furthermore, we emphasized the significance of effective precipitation frequency in regulating carbon fluxes. Our results provide a foundational understanding of the carbon balance in semi-arid ecosystems., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024. Published by Elsevier Ltd.)

Water-table maps are fundamental to hydrogeological studies and a manual, hand-drawn method is still commonly used to produce them. Despite this, the accuracy and variability of such maps have received little attention in international literature. In a unique experiment, 63 groundwater professionals drew water-table equipotential contours based on the same dataset of point measurements and were asked to infer flow directions and predict groundwater elevations at predefined locations. The root mean squared error (RMSE) for the average map calibration data was 10.5 m, which is accuracy comparable to numerical groundwater models. This study confirmed that to produce hand-drawn water-table maps, practitioners seek to not only fit the spatial data, but also to conform to their own cognitive model of hydrogeological concepts and processes. The calibration accuracy increased with experience; from a RMSE of 13.3 m for practitioners with 0-3 years of experience to a RMSE of 9.2 m for those with four or more years. Despite considerable variability in the style of the hand-drawn water-table maps, the maps were consistent in their representation of the dominant regional groundwater flow directions. There was less consensus, however, in predicting the direction of surface water-groundwater interaction for a stream reach. Hand-drawn water-table mapping remains useful and valid, especially as a starting point for hydrogeological conceptualization, yet further work is required to resolve issues around transparency, repeatability, and reproducibility., (© 2024 Commonwealth of Australia and The Authors. Groundwater published by Wiley Periodicals LLC on behalf of National Ground Water Association.)

Rain use efficiency (RUE) quantifies the ecosystem's capacity to use precipitation water to assimilate atmospheric CO 2 . The spatial distribution of RUE and its drivers across the Australian continent is largely unknown. This knowledge gap limits our understanding of the possible contribution of Australian ecosystems to global carbon assimilation. This study investigates the spatial distribution of RUE across diverse terrestrial ecosystems in Australia. The results show that RUE ranges from 0.43 (1st percentile) to 3.10 (99th percentile) g C m -2  mm -1 with a continental mean of 1.19 g C m -2  mm -1 . About 68 % of the spatiotemporal variability of RUE can be explained by a multiple linear regression model primarily contributed by climatic predictors. Benchmarked by the model estimation, drainage-diverging/converging landscapes tend to have reduced/increased RUE. The model also revealed the impact of increasing atmospheric CO 2 concentration on RUE. The continental mean RUE would increase by between 29.3 and 64.8 % by the end of this century under the SSP5-8.5 scenario in which the CO 2 concentration is projected to double from the present level. This increase in projected RUE is attributed to the assumed greening effect of increasing CO 2 concentration, which does not consider the saturation of CO 2 fertilisation effect and the warming effect on increasing wildfire occurrence. Under the SSP1-2.6 scenario, RUE would decrease by about 7 %. This study provides baseline RUEs of various ecosystems in Australia for investigating the impacts of human interferences and climate change on the capacity of Australian vegetation to assimilate atmospheric CO 2 under given precipitation., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2024 The Authors. Published by Elsevier B.V. All rights reserved.)

Under the influence of climate change and human activities, the problem of grassland degradation is becoming increasingly severe. Detection of changes in vegetation cover is crucial for a better understanding of the interaction between humans and ecosystems. This study maps changes in vegetation cover using the Google Earth Engine (GEE). We used 36 years of Landsat satellite imagery (1985-2020) in the Xilin River Basin, China, to classify grassland conditions and validated the results with field observation data. The overall classification of the model accuracy assessment was 83.3%. The Dynamic Reference Vegetation Cover Method (DRCM) was adopted to remove the effect of interannual variation of rainfall, allowing to focus on the impact of human activities on vegetation cover changes. The results identify five categories of vegetation cover changes: significantly increased, potentially increased, stable, potentially decreased, and significantly decreased. The reference level is derived from the most persistent land surface coverage across different grassland types and all years. Overall, 9.3% of the study area had a significant increase in vegetation cover, 14.2% a potential increase, 48.6% of the area showed a stable vegetation condition, 9.8% showed a potential decrease, and 18.1% a significant decrease in vegetation cover. The largest proportion of combined potential and significant reduction was 35.2% for desert grassland, where the vegetation faced the most severe reduction. This study will provide a basis for identifying grassland degradation and developing scientific management policies., Competing Interests: Declaration of competing interest The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper., (Copyright © 2023 Elsevier Ltd. All rights reserved.)

Despite the potential of remote sensing for monitoring reservoir operation, few studies have investigated the extent to which reservoir releases can be inferred across different spatial and temporal scales. Through evaluating 21 reservoirs in the highly regulated Greater Mekong region, remote sensing imagery was found to be useful in estimating daily storage volumes for within-year and over-year reservoirs (correlation coefficients [CC] ≥ 0.9, normalized root mean squared error [NRMSE] ≤ 31%), but not for run-of-river reservoirs (CC < 0.4, 40% ≤ NRMSE ≤ 270%). Given a large gap in the number of reservoirs between global and local databases, the proposed framework can improve representation of existing reservoirs in the global reservoir database and thus human impacts in hydrological models. Adopting an Integrated Reservoir Operation Scheme within a multi-basin model was found to overcome the limitations of remote sensing and improve streamflow prediction at ungauged cascade reservoir systems where previous modeling approaches were unsuccessful. As a result, daily regulated streamflow was predicted competently across all types of reservoirs (median values of CC = 0.65, NRMSE = 8%, and Kling-Gupta efficiency [KGE] = 0.55) and downstream hydrological stations (median values of CC = 0.94, NRMSE = 8%, and KGE = 0.81). The findings are valuable for helping to understand the impacts of reservoirs and dams on streamflow and for developing more useful adaptation measures to extreme events in data sparse river basins., (© 2022. The Authors.)

This study evaluates the potential of kriging-based (kriging and kriging-logistic) and machine learning models (MARS, GBRT, and ANN) in predicting the effluent arsenic concentration of a wastewater treatment plant. Two distinct input combination scenarios were established, using seven quantitative and qualitative independent influent variables. In the first scenario, all of the seven independent variables were taken into account for constructing the data-driven models. For the second input scenario, the forward selection k-fold cross-validation method was employed to select effective explanatory influent parameters. The results obtained from both input scenarios show that the kriging-logistic and machine learning models are effective and robust. However, using the feature selection procedure in the second scenario not only made the architecture of the model simpler and more effective, but also enhanced the performance of the developed models (e.g., around 7.8% performance enhancement of the RMSE). Although the standard kriging method provided the least good predictive results (RMSE = 0.18 ug/l and NSE=0.75), it was revealed that the kriging-logistic method gave the best performance among the applied models (RMSE = 0.11 ug/l and NSE=0.90)., (© 2021. The Author(s), under exclusive licence to Springer-Verlag GmbH Germany, part of Springer Nature.)

There is a paucity of information regarding the interaction between GONPs and natural aquifer sediments. Therefore, batch and column experiments were carried out to determine the transport, retention and attachment behavior of GONPs with the surfaces of native aquifer sediments. The experiments were performed with sediments comprising contrasting mineralogical features (sand grains, quartz and limestone sediments), at different temperatures, ionic strength and compositions. Uniquely, this research also investigated the effect of natural biofilm on the retention behavior of nanoparticles in porous media. The retention rate of GONPs at 22 °C was higher than at 4 °C. Moreover, there was greater retention of GONPs onto the surfaces of collectors at higher ionic strengths and cation valence. The retention profiles (RPs) of GONPs in pristine porous media at low ionic strength were linear, which contrasted with hyper-exponential shape of RPs at high ionic strength. The size-distribution analysis of retained GONPs showed decreasing particle diameter with increasing distance from the column inlet at high ionic strength and equal diameter at low ionic strengths. The GONP retention rate was higher for natural porous media than for sand, due to the presence of metal oxides heterogeneities. The presence of biofilm on porous media increased the retention rate of GONPs when compared to the porous media in the absence of biofilm., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

This research was conducted to evaluate the effect of co-transport of different-sized microorganisms on graphene oxide nanoparticles (GONPs) transport and retention in saturated pristine and biofilm-conditioned limestone columns. The transport and retention behavior of GONPs was studied in columns in the presence of MS2 -as a nano-sized- and Escherichia coli (E.coli) -as a micro-sized- microorganisms at low and high ionic strength conditions. Results showed no changes in GONPs transport and retention at high ionic strength in the presence of MS2 or E. coli, which was attributed to the effect of high concentration of divalent cation on aggregation of nanoparticles and microorganisms. Furthermore, simultaneous enhanced transport and decreased retention of GONPs in column was observed in the co-presence of microorganisms at low ionic strength. Results revealed that the main mechanism governing increasing GONPs transport in porous media was occupation of reactive surface sites of collectors by microorganisms, which prevented attachment of nanoparticles. The pre-saturation of columns with MS2 and E. coli caused increasing transport of GONPs in the columns, due to the occupation of surface reactive sites. Moreover, conditioning limestone collectors with natural biofilm resulted in the same rates of nanoparticle elution and retention (i.e., in the presence or absence of microorganisms) by straining of GONPs in the inlet end of columns which shows that the biofilm acts as a bio-filter against discharging nanoparticles into the effluents. Finally, from the obtained results, it can be postulated that the presence of microorganisms in a MAR site causes risk of groundwater pollution by toxic nanoparticles., Competing Interests: Declaration of competing interest No., (Copyright © 2020 Elsevier Ltd. All rights reserved.)

Background: Virus, as nano-sized microorganisms are prevalent in aquifers, which threaten groundwater quality and human health wellbeing. Virus inactivation by attachment onto the limestone surfaces is a determining factor in the transport and retention behavior of virus in carbonaceous aquifers., Methods: In the present study, the inactivation of MS2 -as a model virus- by attachment onto the surfaces of limestone grains was investigated in a series of batch experiments under different conditions such as limestone particle size distribution (0.25-0.50, 0.5-1 and 1-2 mm), treated wastewater and RO water, temperature (4 and 22 °C), initial MS2 concentrations (10 3 -10 7 PFU/mL) and static and dynamic conditions. The experimental data of MS2 inactivation was also fitted to a non-linear kinetic model with shoulder and tailing. The characteristics of biofilm on the surfaces of limestone aquifer materials were assessed using scanning electron microscopy (SEM) and confocal laser scanning microscopy (CLSM)., Results: The inactivation rate of virus decreased with increasing the adsorbent diameter. Furthermore, virus inactivation was greater at room temperature (22 °C) than 4 °C, in both static and dynamic conditions. The inactivation of virus via attachment onto the limestone aquifer materials in dynamic conditions was higher than under static conditions. In addition, fitting the experimental data with a kinetic model showed that virus inactivation was high at higher temperature, smaller limestone grains and dynamic conditions. Moreover, the experiments with treated wastewater showed that in authentic aqueous media, the virus inactivation was considerably higher than in RO water, due to the presence of either monovalent or divalent cations and surface roughness created by biofilms., Conclusion: Finally, in terms of managed aquifer recharge systems, the presence of biofilm increases bacteria and virus retention onto the aquifer surfaces. Graphical abstract., Competing Interests: Conflict of interestThe authors confirm that they have no real or perceived conflict of interest which might influence the results of the research. The final version of manuscript has been approved by all authors., (© Springer Nature Switzerland AG 2019.)

First-order half-lives for 26 trace organic compounds (TrOCs) were determined in the hyporheic zone (HZ) and along a 3 km reach of a first-order stream in South Australia during both dry and wet seasons. Two salt tracer experiments were conducted and evaluated using a transient storage model to characterize seasonal differences in stream residence time and transient storage. Lagrangian and time-integrated surface water sampling were conducted to calculate half-lives in the surface water. Half-lives in the HZ were calculated using porewater samples obtained from a modified mini-point sampler and hyporheic residence times measured via active heat-pulse sensing. Half of the investigated TrOCs (e.g., oxazepam, olmesartan, candesartan) were not significantly removed along both the investigated river stretch and the sampled hyporheic flow paths. The remaining TrOCs (e.g., metformin, guanylurea, valsartan) were found to be significantly removed in the HZ and along the river stretch with relative removals in the HZ correlating to reach-scale relative removals. Using the modeled transport parameters, it was estimated that wet season reach-scale removal of TrOCs was predominately caused by removal in the HZ when the intensity of hyporheic exchange was also higher. Factors that increase HZ exchange are thus likely to promote in-stream reactivity of TrOCs.

Target 6.4 of the recently adopted Sustainable Development Goals (SDGs) deals with the reduction of water scarcity. To monitor progress towards this target, two indicators are used: Indicator 6.4.1 measuring water use efficiency and 6.4.2 measuring the level of water stress (WS). This paper aims to identify whether the currently proposed indicator 6.4.2 considers the different elements that need to be accounted for in a WS indicator. WS indicators compare water use with water availability. We identify seven essential elements: 1) both gross and net water abstraction (or withdrawal) provide important information to understand WS; 2) WS indicators need to incorporate environmental flow requirements (EFR); 3) temporal and 4) spatial disaggregation is required in a WS assessment; 5) both renewable surface water and groundwater resources, including their interaction, need to be accounted for as renewable water availability; 6) alternative available water resources need to be accounted for as well, like fossil groundwater and desalinated water; 7) WS indicators need to account for water storage in reservoirs, water recycling and managed aquifer recharge. Indicator 6.4.2 considers many of these elements, but there is need for improvement. It is recommended that WS is measured based on net abstraction as well, in addition to currently only measuring WS based on gross abstraction. It does incorporate EFR. Temporal and spatial disaggregation is indeed defined as a goal in more advanced monitoring levels, in which it is also called for a differentiation between surface and groundwater resources. However, regarding element 6 and 7 there are some shortcomings for which we provide recommendations. In addition, indicator 6.4.2 is only one indicator, which monitors blue WS, but does not give information on green or green-blue water scarcity or on water quality. Within the SDG indicator framework, some of these topics are covered with other indicators., (Copyright © 2017 The Authors. Published by Elsevier B.V. All rights reserved.)

Cone penetration testing (CPT) is one of the most efficient and versatile methods currently available for geotechnical, lithostratigraphic and hydrogeological site characterization. Currently available methods for soil behaviour type classification (SBT) of CPT data however have severe limitations, often restricting their application to a local scale. For parameterization of regional groundwater flow or geotechnical models, and delineation of regional hydro- or lithostratigraphy, regional SBT classification would be very useful. This paper investigates the use of model-based clustering for SBT classification, and the influence of different clustering approaches on the properties and spatial distribution of the obtained soil classes. We additionally propose a methodology for automated lithostratigraphic mapping of regionally occurring sedimentary units using SBT classification. The methodology is applied to a large CPT dataset, covering a groundwater basin of ~60 km2 with predominantly unconsolidated sandy sediments in northern Belgium. Results show that the model-based approach is superior in detecting the true lithological classes when compared to more frequently applied unsupervised classification approaches or literature classification diagrams. We demonstrate that automated mapping of lithostratigraphic units using advanced SBT classification techniques can provide a large gain in efficiency, compared to more time-consuming manual approaches and yields at least equally accurate results.

Drought is an intermittent disturbance of the water cycle that profoundly affects the terrestrial carbon cycle. However, the response of the coupled water and carbon cycles to drought and the underlying mechanisms remain unclear. Here we provide the first global synthesis of the drought effect on ecosystem water use efficiency (WUE = gross primary production (GPP)/evapotranspiration (ET)). Using two observational WUE datasets (i.e., eddy-covariance measurements at 95 sites (526 site-years) and global gridded diagnostic modelling based on existing observation and a data-adaptive machine learning approach), we find a contrasting response of WUE to drought between arid (WUE increases with drought) and semi-arid/sub-humid ecosystems (WUE decreases with drought), which is attributed to different sensitivities of ecosystem processes to changes in hydro-climatic conditions. WUE variability in arid ecosystems is primarily controlled by physical processes (i.e., evaporation), whereas WUE variability in semi-arid/sub-humid regions is mostly regulated by biological processes (i.e., assimilation). We also find that shifts in hydro-climatic conditions over years would intensify the drought effect on WUE. Our findings suggest that future drought events, when coupled with an increase in climate variability, will bring further threats to semi-arid/sub-humid ecosystems and potentially result in biome reorganization, starting with low-productivity and high water-sensitivity grassland.

Shale gas is currently being explored in Europe as an alternative energy source to conventional oil and gas. There is, however, increasing concern about the potential environmental impacts of shale gas extraction by hydraulic fracturing (fracking). In this study, we focussed on the potential impacts on regional water resources within the Baltic Basin in Poland, both in terms of quantity and quality. The future development of the shale play was modeled for the time period 2015-2030 using the LUISA modeling framework. We formulated two scenarios which took into account the large range in technology and resource requirements, as well as two additional scenarios based on the current legislation and the potential restrictions which could be put in place. According to these scenarios, between 0.03 and 0.86% of the total water withdrawals for all sectors could be attributed to shale gas exploitation within the study area. A screening-level assessment of the potential impact of the chemicals commonly used in fracking was carried out and showed that due to their wide range of physicochemical properties, these chemicals may pose additional pressure on freshwater ecosystems. The legislation put in place also influenced the resulting environmental impacts of shale gas extraction. Especially important are the protection of vulnerable ground and surface water resources and the promotion of more water-efficient technologies.

This paper presents the Kalman Filtered Double Constraint Method (DCM-KF) as a technique to estimate the hydraulic conductivities in the grid blocks of a groundwater flow model. The DCM is based on two forward runs with the same initial grid block conductivities, but with alternating flux-head conditions specified on parts of the boundary and the wells. These two runs are defined as: (1) the flux run, with specified fluxes (recharge and well abstractions), and (2) the head run, with specified heads (measured in piezometers). Conductivities are then estimated as the initial conductivities multiplied by the fluxes obtained from the flux run and divided by the fluxes obtained from the head run. The DCM is easy to implement in combination with existing models (e.g., MODFLOW). Sufficiently accurate conductivities are obtained after a few iterations. Because of errors in the specified head-flux couples, repeated estimation under varying hydrological conditions results in different conductivities. A time-independent estimate of the conductivities and their inaccuracy can be obtained by a simple linear KF with modest computational requirements. For the Kleine Nete catchment, Belgium, the DCM-KF yields sufficiently accurate calibrated conductivities. The method also results in distinguishing regions where the head-flux observations influence the calibration from areas where it is not able to influence the hydraulic conductivity., (© 2014, National Ground Water Association.)

Diffuse phosphorus (P) export from agricultural land to surface waters is a significant environmental problem. It is critical to determine the natural background P losses from diffuse sources, but their identification and quantification is difficult. In this study, three headwater catchments with differing land use (arable, pasture and forest) were monitored for 3 years to quantify exports of dissolved (<0.45 µm) reactive P and total dissolved P. Mean total P exports from the arable catchment ranged between 0.08 and 0.28 kg ha(-1) year(-1). Compared with the reference condition (forest), arable land and pasture exported up to 11-fold more dissolved P. The contribution of dissolved (<0.45 µm) unreactive P was low to negligible in every catchment. Agricultural practices can exert large pressures on surface waters that are controlled by hydrological factors. Adapting policy to cope with these factors is needed for lowering these pressures in the future.

The estimation of surface-subsurface water interactions is complex and highly variable in space and time. It is even more complex when it has to be estimated in urban areas, because of the complex patterns of the land-cover in these areas. In this research a modeling approach with integrated remote sensing analysis has been developed for estimating water fluxes in urban environments. The methodology was developed with the aim to simulate fluxes of contaminants from polluted sites. Groundwater pollution in urban environments is linked to patterns of land use and hence it is essential to characterize the land cover in a detail. An object-oriented classification approach applied on high-resolution satellite data has been adopted. To assign the image objects to one of the land-cover classes a multiple layer perceptron approach was adopted (Kappa of 0.86). Groundwater recharge has been simulated using the spatially distributed WetSpass model and the subsurface water flow using MODFLOW in order to identify and budget water fluxes. The developed methodology is applied to a brownfield case site in Vilvoorde, Brussels (Belgium). The obtained land use map has a strong impact on the groundwater recharge, resulting in a high spatial variability. Simulated groundwater fluxes from brownfield to the receiving River Zenne were independently verified by measurements and simulation of groundwater-surface water interaction based on thermal gradients in the river bed. It is concluded that in order to better quantify total fluxes of contaminants from brownfields in the groundwater, remote sensing imagery can be operationally integrated in a modeling procedure., (Copyright © 2010 Elsevier B.V. All rights reserved.)

Continental export of Si to the coastal zone is closely linked to the ocean carbon sink and to the dynamics of phytoplankton blooms in coastal ecosystems. Presently, however, the impact of human cultivation of the landscape on terrestrial Si fluxes remains unquantified and is not incorporated in models for terrestrial Si mobilization. In this paper, we show that land use is the most important controlling factor of Si mobilization in temperate European watersheds, with sustained cultivation (>250 years) of formerly forested areas leading to a twofold to threefold decrease in baseflow delivery of Si. This is a breakthrough in our understanding of the biogeochemical Si cycle: it shows that human cultivation of the landscape should be recognized as an important controlling factor of terrestrial Si fluxes.

Remote sensing offers a cost efficient means for identifying and monitoring wetlands over a large area and at different moments in time. In this study, we aim at providing ecologically relevant information on characteristics of temporary and permanent isolated open water wetlands, obtained by standard techniques and relatively cheap imagery. The number, surface area, nearest distance, and dynamics of isolated temporary and permanent wetlands were determined for the Western Cape, South Africa. Open water bodies (wetlands) were mapped from seven Landsat images (acquired during 1987 - 2002) using supervised maximum likelihood classification. The number of wetlands fluctuated over time. Most wetlands were detected in the winter of 2000 and 2002, probably related to road constructions. Imagery acquired in summer contained fewer wetlands than in winter. Most wetlands identified from Landsat images were smaller than one hectare. The average distance to the nearest wetland was larger in summer. In comparison to temporary wetlands, fewer, but larger permanent wetlands were detected. In addition, classification of non-vegetated wetlands on an Envisat ASAR radar image (acquired in June 2005) was evaluated. The number of detected small wetlands was lower for radar imagery than optical imagery (acquired in June 2002), probably because of deterioration of the spatial information content due the extensive pre-processing requirements of the radar image. Both optical and radar classifications allow to assess wetland characteristics that potentially influence plant and animal metacommunity structure. Envisat imagery, however, was less suitable than Landsat imagery for the extraction of detailed ecological information, as only large wetlands can be detected. This study has indicated that ecologically relevant data can be generated for the larger wetlands through relatively cheap imagery and standard techniques, despite the relatively low resolution of Landsat and Envisat imagery. For the characterisation of very small wetlands, high spatial resolution optical or radar images are needed. This study exemplifies the benefits of integrating remote sensing and ecology and hence stimulates interdisciplinary research of isolated wetlands.

The amount and intensity of runoff on catchment scale are strongly determinedby the presence of impervious land-cover types, which are the predominant cover types inurbanized areas. This paper examines the impact of different methods for estimatingimpervious surface cover on the prediction of peak discharges, as determined by a fullydistributed rainfall-runoff model (WetSpa), for the upper part of the Woluwe Rivercatchment in the southeastern part of Brussels. The study shows that detailed informationon the spatial distribution of impervious surfaces, as obtained from remotely sensed data,produces substantially different estimates of peak discharges than traditional approachesbased on expert judgment of average imperviousness for different types of urban land use.The study also demonstrates that sub-pixel estimation of imperviousness may be a usefulalternative for more expensive high-resolution mapping for rainfall-runoff modelling atcatchment scale.

The increasing anthropogenic pressure on natural environments results in impacts that affect tropical forest areas and their biodiversity. Adverse impacts on terrestrial and oceanic environments often compound in the intertidal area, where mangrove forest ecosystems thrive. In tropical coastal areas of many developing countries where people depend on wood and other mangrove forest products and services, forest degradation leads to socioeconomic problems. At the same time, increasing freshwater needs in these areas are expected to cause additional problems. On the basis of remote sensing and ground truthing complemented by colonial archival material from the Dutch East India Company (1602-1800), we report that changes to the historic system of inland freshwater management have increased dramatically in recent times. Hydrological changes, such as interbasin transfers, have resulted in a qualitative ecological and socioeconomic degradation in three coastal lagoons in southern Sri Lanka. Variations in river hydrology have caused changes in the areas suitable as mangrove habitat and, thus, have resulted in an altered distribution. However, increases in mangrove area can mask the degradation of the site in terms of floristic composition, significance of the species, and biodiversity (this effect is termed "cryptic ecological degradation"). It is important that such changes be carefully monitored to ensure biological and socioeconomic sustainability.

The prediction of the location of ground water discharge areas is a key aspect for the protection and (re)development of ground water-dependent wetlands. Ground water discharge areas can be simulated with MODFLOW using the DRAIN package by setting the drain level equal to the topography, while the conductance is mostly set to an arbitrary high value. However, conceptual and practical problems arise in the calculation of the ground water discharge by the DRAIN package as calculated water tables above the land surface, difficult parameterization of the conductance, and large water balance errors. To overcome these problems, a new SEEPAGE package for MODFLOW is proposed. The basic idea of this package is an adaptable constant head cell. It has a variable head, unless the ground water rises above the seepage level, in which case it has a constant head cell. The estimation of the ground water discharge location along a homogeneous, isotropic, linear sloping profile is used to verify the model and to compare it to the DRAIN solution. In an application to three basins in Belgium, it is shown that the SEEPAGE package can be used in combination with the DRAIN package in situations where an upper boundary for a free water table and additional resistance for drainage is required. It is clearly demonstrated that the identification and delineation of regional ground water discharge areas is more accurate using the SEEPAGE package.

This paper presents a Bayesian Monte Carlo method for evaluating the uncertainty in the delineation of well capture zones and its application to a wellfield in a heterogeneous, multiaquifer system. In the method presented, Bayes' rule is used to update prior distributions for the unknown parameters of the stochastic model for the hydraulic conductivity, and to calculate probability-based weights for parameter realizations using head residuals. These weights are then assigned to the corresponding capture zones obtained using forward particle tracking. Statistical analysis of the set of weighted protection zones results in a probability distribution for the capture zones. The suitability of the Bayesian stochastic method for a multilayered system is investigated, using the wellfield Het Rot at Nieuwrode, Belgium, located in a three-layered aquifer system, as an example. The hydraulic conductivity of the production aquifer is modeled as a spatially correlated random function with uncertain parameters. The aquitard and overlying unconfined aquifer are assigned random, homogeneous conductivities. The stochastic results are compared with deterministic capture zones obtained with a calibrated model for the area. The predictions of the stochastic approach are more conservative and indicate that parameter uncertainty should be taken into account in the delineation of well capture zones.