Plants versus streams: Their groundwater-mediated competition at “El Morro,” a developing catchment in the dry plains of Argentina

Our understanding of how groundwater mediates evapotranspiration/streamflow partitioning is still fragmented and catchment studies under changing vegetation conditions can provide a useful frame for integration. We explored this partition in a flat sedimentary dry catchment in central Argentina in which the replacement of native vegetation with rainfed crops was accompanied by the abrupt formation of groundwater-fed streams by subsurface erosion (i.e., sapping) episodes. Historical records indicated widespread water table rises (~0.3 m y⁻¹ on average). Groundwater level and stream baseflow fluctuated seasonally with minima in the warm rainy season, indicating that evaporative discharge rather than rainfall shapes saturated flows. Diurnal groundwater level fluctuations showed that plant uptake was widespread where water tables are shallow (<3 m) but restricted to deep-rooted Prosopis forests where they are deep (7–10 m). MODIS and LANDSAT NDVI revealed a long-term greening for native vegetation, new wetlands included, but not for croplands, suggesting more limited evapotranspiration-groundwater level regulation under agriculture. Close to the deepest (20 m) and most active incisions, groundwater level and greenness declined and stream baseflow showed no seasonal fluctuations, hinting decoupling from evapotranspiration. Intense ecological and geomorphological transformations in this catchment exposed the interplay of five mechanisms governing evapotranspiration/streamflow partition including (a) unsaturated uptake and both (b) riparian and (c) distributed uptake from the saturated zone by plants, as well as (d) deepening incisions and (e) sediment deposits over riparian zones by streams. Acknowledging the complex interplay of these mechanisms with groundwater is crucial to predict and manage future hydrological changes in the dry plains of South America.     

Ecohydrological evolution model on riparian vegetation in hyperarid regions and its validation in the lower reach of Tarim River

The evolution of the ecohydrological system driven by external climatic forcing and internal feedbacks between vegetation and hydrology, which is more remarkable in arid and semiarid regions, has attracted substantial research attention in recent years. To examine critically the state‐of‐the‐art assumptions and dynamic equations used in the evolution study of an ecohydrological system, the rule of proceeding from simplicity to complexity should be followed. The riparian vegetation ecohydrological system in hyperarid regions (e.g. the lower Tarim River) can serve as a starting point given its simplicity, which has been seldom examined before in terms of system evolution. Further, the water transfer practice from 2000 to 2006 in the lower Tarim River serves as a valuable prototype experiment for model validation. This is because the remarkable changes in groundwater and vegetation in the area have taken place within a shorter period and thus can be easily observed. In the present study, the ecohydrological evolution model on riparian vegetation (ERV model) in hyperarid regions was proposed by coupling groundwater movement and vegetation dynamics. In the ERV model, the groundwater table serves as a critical feedback variable that determines the vegetation dynamics (colonization and mortality) and is determined by vegetation transpiration other than groundwater movement. The monitored groundwater table by wells and satellite‐observed vegetation coverage from the Moderate Resolution Imaging Spectroradiometer are used for model validation. The simulation results show the good performance of the ERV model with uncalibrated parameters. It was also calibrated manually using a multiobjective method, and the fine‐tuned parameters are close to the uncalibrated ones, indicating the robustness of the model. The analysis shows further that the increased evapotranspiration is substantially due to the water transfer and thus the vegetation growth, which implies the importance of ecohydrological coupling for long‐term hydrological modelling. Copyright © 2011 John Wiley & Sons, Ltd.     

Effect of elevation resolution on evapotranspiration simulations using MODFLOW

Surface elevations represented in MODFLOW head-dependent packages are usually derived from digital elevation models (DEMs) that are available at much high resolution. Conventional grid refinement techniques to simulate the model at DEM resolution increases computational time, input file size, and in many cases are not feasible for regional applications. This research aims at utilizing the increasingly available high resolution DEMs for effective simulation of evapotranspiration (ET) in MODFLOW as an alternative to grid refinement techniques. The source code of the evapotranspiration package is modified by considering for a fixed MODFLOW grid resolution and for different DEM resolutions, the effect of variability in elevation data on ET estimates. Piezometric head at each DEM cell location is corrected by considering the gradient along row and column directions. Applicability of the research is tested for the lower Rio Grande (LRG) Basin in southern New Mexico. The DEM at 10 m resolution is aggregated to resampled DEM grid resolutions which are integer multiples of MODFLOW grid resolution. Cumulative outflows and ET rates are compared at different coarse resolution grids. Results of the analysis conclude that variability in depth-to-groundwater within the MODFLOW cell is a major contributing parameter to ET outflows in shallow groundwater regions. DEM aggregation methods for the LRG Basin have resulted in decreased volumetric outflow due to the formation of a smoothing error, which lowered the position of water table to a level below the extinction depth.     

Use of Nested Flow Models and Interpolation Techniques for Science‐Based Management of the Sheyenne National Grassland,North Dakota,USA

Noxious weeds threaten the Sheyenne National Grassland (SNG) ecosystem and therefore herbicides have been used for control. To protect groundwater quality, the herbicide application is restricted to areas where the water table is less than 10 feet (3.05 m) below the ground surface in highly permeable soils, or less than 6 feet (1.83 m) below the ground surface in low permeable soils. A local MODFLOW model was extracted from a regional GFLOW analytic element model and used to develop depth‐to‐groundwater maps in the SNG that are representative for the particular time frame of herbicide applications. These maps are based on a modeled groundwater table and a digital elevation model (DEM). The accuracy of these depth‐to‐groundwater maps is enhanced by an artificial neural networks (ANNs) interpolation scheme that reduces residuals at 48 monitoring wells. The combination of groundwater modeling and ANN improved depth‐to‐groundwater maps, which in turn provided more informed decisions about where herbicides can or cannot be safely applied.     

Extinction depth and evapotranspiration from ground water under selected land covers

In many landscapes, vegetation extracts water from both the unsaturated and the saturated zones. The partitioning of evapotranspiration (ET) into vadose zone evapotranspiration and ground water evapotranspiration (GWET) is complex because it depends on land cover and subsurface characteristics. Traditionally, the GWET fraction is assumed to decay with increasing depth to the water table (DTWT), attaining a value of 0 at what is termed the extinction depth. A simple assumption of linear decay with depth is often used but has never been rigorously examined using unsaturated-saturated flow simulations. Furthermore, it is not well understood how to relate extinction depths to characteristics of land cover and soil texture. In this work, variable saturation flow theory is used to simulate GWET for three land covers and a range of soil properties under drying soil conditions. For a water table within half a meter of the land surface, nearly all ET is extracted from ground water due to the close hydraulic connection between the unsaturated and the saturated zones. For deep-rooted vegetation, the decoupling of ground water and vadose zone was found to begin at water table depths between 30 and 100 cm, depending on the soil texture. The decline of ET with DTWT is better simulated by an exponential decay function than the commonly used linear decay. A comparison with field data is consistent with the findings of this study. Tables are provided to vary the extinction depth for heterogeneous landscapes with different vegetation cover and soil properties.     

Analytical Solutions of Travel Time to a Pumping Well with Variable Evapotranspiration

Analytical solutions of groundwater travel time to a pumping well in an unconfined aquifer have been developed in previous studies, however, the change in evapotranspiration was not considered. Here, we develop a mathematical model of unconfined flow toward a discharge well with redistribution of groundwater evapotranspiration for travel time analysis. Dependency of groundwater evapotranspiration on the depth to water table is described using a linear formula with an extinction depth. Analytical solutions of groundwater level and travel time are obtained. For a typical hypothetical example, these solutions perfectly agree with the numerical simulation results based on MODFLOW and MODPATH. As indicated in a dimensionless framework, a lumped parameter which is proportional to the pumping rate controls the distributions of groundwater evapotranspiration rate and the travel time along the radial direction.     

Evapotranspiration and its main controlling mechanism over the desert riparian forests in the lower Tarim River Basin

Evapotranspiration(ET) and its controlling mechanism over the desert riparian forests in arid regions are the important scientific basis for the water resources managements of the lower reaches of the inland rivers of China. Nearly three years of continuous measurements of surface ET, soil water content at different depths and groundwater table over a typical Tamarix spp. stand and a typical Populus euphratica stand were conducted in the lower reach of the Tarim River. The ET seasonal trends in the growing season were controlled by plant phenology, and ET in non-growing season was weak. The diurnal variations of ET resulting from the comprehensive effects of all atmospheric factors were significantly related with reference ET. The spatial pattern of ET was determined by vegetation LAI, more vegetation coverage, more ET amount. Groundwater is the water source of surface ET, and the soil water in shallow layers hardly took part in the water exchange in the groundwatersoil-plant-air system. The temporal processes of ET over the Tamarix stand and the Populus stand were similar, but the water consumption of the well-grown Populus euphratica was higher than that of the well-grown Tamarix spp. Further analysis indicates that plant transpiration accounts for most of the surface ET, with soil evaporation weak and negligible; groundwater table is a crucial factor influencing ET over the desert riparian forests, groundwater influences the processes and amounts of ET by controlling the growth and spatial distribution of desert riparian forests; quantifying the water stress of desert riparian forests using groundwater table is more appropriate, rather than soil water content. Based on the understanding of ET and water movements in the groundwater-soil-plant-air system, a generalized framework expressing the water cycling and its key controlling mechanism in the lower reaches of the inland rivers of China is described, and a simple model to estimate water requirements of the desert riparian forests is presented.     

Shallow groundwater dynamics and its driving forces in extremely arid areas: a case study of the lower Heihe River in northwestern China

Shallow groundwater is an important source of water for the maintenance and restoration of ecosystems in arid environments, which necessitates a deeper understanding of its complex spatial and temporal dynamics driven by hydrological processes. This study explores the dominant hydrological processes that control the shallow groundwater dynamics in the Gobi Desert‐riparian‐oasis system of the lower Heihe River, a typical arid inland river basin located in northwestern China. The groundwater level and temperature were monitored in 14 shallow wells at 30‐min intervals during the 2010–2012 period. After combining this information with meteorological and hydrological data, a comprehensive analysis was conducted to understand the dynamic behaviour of the shallow groundwater system and to determine the dominant factors that control the groundwater flow processes. The results of the study indicate notably large temporal and spatial variations in both the groundwater level and temperature. Noticeable fluctuations in the groundwater level (0.5–1 m) and temperature (4–8 °C) were observed in the riparian zone, evidencing a clear river influence. In comparison, the groundwater fluctuations in the Gobi Desert were more stable (the annual variations of the water table were less than 0.5 m, and the water temperature varied by no more than 2 °C). Strong variations in the groundwater table (1.5–5.0 m/year) and temperature (1.5–6.5 °C), mainly caused by surface flood irrigation and groundwater pumping, were observed in the oasis area. The investigated sites were categorized into three types that reflect the dominant hydrological processes: (1) the riparian zone, dominated by riverbank filtration and groundwater evapotranspiration; (2) the Gobi Desert area, controlled by groundwater evaporation and lateral recharge; and (3) the oasis area, dominated by groundwater evapotranspiration as well as surface–groundwater interactions caused by human activities. Copyright © 2012 John Wiley & Sons, Ltd.     

Relationships between riparian evapotranspiration and groundwater depth along a semiarid irrigated river valley

Evapotranspiration (ET) from riparian vegetation can be difficult to estimate due to relatively abundant water supply, spatial vegetation heterogeneity, and interactions with anthropogenic influences such as shallower groundwater tables, increased salinity, and nonpoint source pollution induced by irrigation. In semiarid south-eastern Colorado, reliable ET estimates are scarce for the riparian corridor that borders the Arkansas River. This work investigates relationships between the riparian ecosystem along the Arkansas River and an underlying alluvial aquifer using ET estimates from remotely sensed data and modelled water table depths. Results from a calibrated, finite-difference groundwater model are used to estimate weekly water table fluctuations in the riparian ecosystem from 1999 to 2009, and estimates of ET are calculated using the Operational Simplified Surface Energy Balance (SSEBop) model with over 200 Landsat scenes covering over 30 km² of riparian ecosystem along a 70-km stretch of the river. Comparison of calculated monthly SSEBop ET to estimated alfalfa reference ET from local micrometeorological station data indicated statistically significant high linear correspondence (R² = .87). Daily calculated SSEBop ET showed statistically significant moderate linear correspondence with data from a local weighing lysimeter (R² = .59). Simulated monthly SSEBop ET values were larger in drier years compared with wetter years, and ET variability was also larger in drier years. Peak ET most commonly occurred during the month of June for all 11 years of analysis. Relationships between ET and water table depth showed that peak monthly ET was highest when groundwater depths were less than about 3 m, and ET values were significantly lower for groundwater depths greater than 3 m. Negative sample Spearman correlation highlighted riparian corridor locations where ET increased as a result of decreased groundwater depths across years with different hydroclimatic conditions. This study shows how a combination of remotely sensed riparian ET estimates and a regional groundwater model can improve our understanding of linkages between riparian consumptive use and near-river groundwater conditions influenced by irrigation return flow and different climatic drivers.     

IHMS—Integrated Hydrological Modelling System. Part 1. Hydrological processes and general structure

A newly Integrated Hydrological Modelling System (IHMS) has been developed to study the impact of changes in climate, land use and water management on groundwater and seawater intrusion (SWI) into coastal areas. The system represents the combination of three models, which can, if required, be run separately. It has been designed to assess the combined impact of climate, land use and groundwater abstraction changes on river, drainage and groundwater flows, groundwater levels and, where appropriate, SWI. The approach is interdisciplinary and reflects an integrated water management approach. The system comprises three packages: the Distributed Catchment Scale Model (DiCaSM), MODFLOW (96 and 2000) and SWI models. In addition to estimating all water balance components, DiCaSM, produces the recharge data that are used as input to the groundwater flow model of the US Geological Survey, MODFLOW. The latter subsequently generates the head distribution and groundwater flows that are used as input to the SWI model, SWI. Thus, any changes in land use, rainfall, water management, abstraction, etc. at the surface are first handled by DiCaSM, then by MODFLOW and finally by the SWI. The three models operate at different spatial and temporal scales and a facility (interface utilities between models) to aggregate/disaggregate input/output data to meet a desired spatial and temporal scale was developed allowing smooth and easy communication between the three models. As MODFLOW and SWI are published and in the public domain, this article focuses on DiCaSM, the newly developed unsaturated zone DiCaSM and equally important the interfacing utilities between the three models. DiCaSM simulates a number of hydrological processes: rainfall interception, evapotranspiration, surface runoff, infiltration, soil water movement in the root zone, plant water uptake, crop growth, stream flow and groundwater recharge. Input requirements include distributed data sets of rainfall, land use, soil types and digital terrain; climate data input can be either distributed or non‐distributed. The model produces distributed and time series output of all water balance components including potential evapotranspiration, actual evapotranspiration, rainfall interception, infiltration, plant water uptake, transpiration, soil water content, soil moisture (SM) deficit, groundwater recharge rate, stream flow and surface runoff. This article focuses on details of the hydrological processes and the various equations used in DiCaSM, as well as the nature of the interface to the MODFLOW and SWI models. Furthermore, the results of preliminary tests of DiCaSM are reported; these include tests related to the ability of the model to predict the SM content of surface and subsurface soil layers, as well as groundwater levels. The latter demonstrates how the groundwater recharge calculated from DiCaSM can be used as input into the groundwater model MODFLOW using aggregation and disaggregation algorithms (built into the interface utility). SWI has also been run successfully with hypothetical examples and was able to reproduce the results of some of the original examples of Bakker and Schaars ( 2005 ). In the subsequent articles, the results of applications to different catchments will be reported. Copyright © 2010 John Wiley & Sons, Ltd.     

Groundwater system response at sites with different agricultural land uses: case of the Guarani Aquifer outcrop area,Brotas/SP-Brazil

A system identification approach can be incorporated in groundwater time series analysis, revealing information concerning the local hydrogeological situation. The aim of this work was to analyse water table fluctuations in an outcrop area of the Guarani Aquifer System (GAS) in Brotas/SP, Brazil, using data from a groundwater monitoring network. The water table dynamic was modelled using continuous time series models that reference the hydrogeological system upon which they operate. The model’s climatological inputs of precipitation and evapotranspiration generate impulse response (IR) functions with parameters that can be related to the physical conditions concerning the hydrological processes involved. The interpretation of the model parameters from two sets of monitoring wells selected at different land-use sites is presented, exemplifying the effect of different water table depths and the distance to the nearest drainage location. Systematic trends of water table depths were also identified from model parameters at specific periods and related to plant development, crop harvest and land-use changes.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR L. Ruiz     

Using MODFLOW 2000 to model ET and recharge for shallow ground water problems

In environments with shallow ground water elevation, small changes in the water table can cause significant variations in recharge and evapotranspiration fluxes. Particularly, where ground water is close to the soil surface, both recharge and evapotranspiration are regulated by a thin unsaturated zone and, for accuracy, must be represented using nonconstant and often nonlinear relationships. The most commonly used ground water flow model today, MODFLOW, was originally designed with a modular structure with independent packages representing recharge and evaporation processes. Systems with shallow ground water, however, may be better represented using either a recharge function that varies with ground water depth or a continuous recharge and evapotranspiration function that is dependent on depth to water table. In situations where the boundaries between recharging and nonrecharging cells change with time, such as near a seepage zone, a continuous ground water flux relationship allows recharge rates to change with depth rather than having to calculate them at each stress period. This research article describes the modification of the MODFLOW 2000 recharge and segmented evapotranspiration packages into a continuous recharge-discharge function that allows ground water flux to be represented as a continuous process, dependent on head. The modifications were then used to model long-term recharge and evapotranspiration processes on a saline, semiarid floodplain in order to understand spatial patterns of salinization, and an overview of this process is given.     

Pathlines of pollutant migration in groundwater of the Visakhapatnam urban area,India

The groundwater regime in the north‐western part of the Visakhapatnam urban area was polluted as early as 1981 by discharge of untreated industrial effluent from a Hindustan Polymers Limited (HPL) plant. A total dissolved solids (TDS) concentration of surficial effluent up to 6500 mg/l and of groundwater in the range 3000–4200 mg/l has been reported in the environs of the HPL plant during May 1992. Groundwater occurs under water table conditions. The groundwater flow model was simulated as a single layer aquifer, using MODFLOW and FLOWPATH computer models. Aquifer parameters were estimated and the value of effective porosity assumed. The water table configuration was mapped in May 1981 using 33 observation wells. The calibrated hydraulic head distribution was used to compute the velocity field using the effective porosity values. The flow paths of groundwater migration from the source locations up to 2002 AD were predicted. Pathlines of particles in the groundwater indicated a predominant north‐east and south‐west migration of groundwater pollution in the area. Copyright © 1999 John Wiley & Sons, Ltd.     

Groundwater use by vegetation in a tropical savanna riparian zone (Daly River, Australia)

Soil water matric potentials (Ψ_m) and the deuterium (δ²H) composition at natural abundance levels of xylem water, soil water, river water and groundwater were used to evaluate whether trees use groundwater during the dry season in the riparian zone of the Daly River (Northern Territory, Australia). Groundwater was a significant source of water for plant transpiration, probably accounting for more than 50% of the water transpired during the dry season. Groundwater use occurred either when trees used water from the capillary fringe or when low Ψ_m induced by soil water uptake lifted groundwater in the vadose zone. Several water use strategies were inferred within the riparian plant community. Melaleuca argentea W. Fitzg and Barringtonia acutangula (L.) Gaertn. appeared to be obligate phreatophytes as they used groundwater almost exclusively and were associated with riverbanks and lower terraces with shallow (<5 m) water tables. Several species appeared to be facultative phreatophytes (including Cathorium umbellatum (Vahl.) Kosterm. and Acacia auriculiformis A. Cunn. ex Benth.) and tended to rely more heavily on soil water with increased elevation in the riparian zone. The levee-bound Corymbia bella K.D. Hill and L.A.S. Johnson mostly used soil water and is either a facultative phreatophyte or a non-phreatophyte. The temporal variability in groundwater utilisation by the trees is unclear because the study focused on the end of the dry season only. A decline in the regional water table as a result of groundwater pumping may affect the health of riparian zone vegetation in the Daly River because groundwater use is significant during the dry season.     

Short‐ and long‐term evapotranspiration rates at ecological restoration sites along a large river receiving rare flow events

Many large rivers around the world no longer flow to their deltas, due to ever greater water withdrawals and diversions for human needs. However, the importance of riparian ecosystems is drawing increasing recognition, leading to the allocation of environmental flows to restore river processes. Accurate estimates of riparian plant evapotranspiration (ET) are needed to understand how the riverine system responds to these rare events and achieve the goals of environmental flows. In 2014, historic environmental flows were released into the Lower Colorado River at Morelos Dam (Mexico); this once perennial but now dry reach is the final stretch to the mighty Colorado River Delta. One of the primary goals was to supply native vegetation restoration sites along the reach with water to help seedlings establish and boost groundwater levels to foster the planted saplings. Patterns in ET before, during, and after the flows are useful for evaluating whether this goal was met and understanding the role that ET plays in this now ephemeral river system. Here, diurnal fluctuations in groundwater levels and Moderate Resolution Imaging Spectroradiometer (MODIS) data were used to compare estimates of ET specifically at 3 native vegetation restoration sites during 2014 planned flow events, and MODIS data were used to evaluate long‐term (2002–2016) ET responses to restoration efforts at these sites. Overall, ET was generally 0–10 mm d^?1 across sites, and although daily ET values from groundwater data were highly variable, weekly averaged estimates were highly correlated with MODIS‐derived estimates at most sites. The influence of the 2014 flow events was not immediately apparent in the results, although the process of clearing vegetation and planting native vegetation at the restoration sites was clearly visible in the results.     

Discontinuous Steady‐State Analytical Solutions of the Boussinesq Equation and Their Numerical Representation by Modflow

Known analytical solutions of groundwater flow equations are routinely used for verification of computer codes. However, these analytical solutions (e.g., the Dupuit solution for the steady‐state unconfined unidirectional flow in a uniform aquifer with a flat bottom) represent smooth and continuous water table configurations, simulating which does not pose any significant problems for the numerical groundwater flow models, like MODFLOW. One of the most challenging numerical cases for MODFLOW arises from drying‐rewetting problems often associated with abrupt changes in the elevations of impervious base of a thin unconfined aquifer. Numerical solutions of groundwater flow equations cannot be rigorously verified for such cases due to the lack of corresponding exact analytical solutions. Analytical solutions of the steady‐state Boussinesq equation, associated with the discontinuous water table configurations over a stairway impervious base, are presented in this article. Conditions resulting in such configurations are analyzed and discussed. These solutions appear to be well suited for testing and verification of computer codes. Numerical solutions, obtained by the latest versions of MODFLOW (MODFLOW‐2005 and MODFLOW‐NWT), are compared with the presented discontinuous analytical solutions. It is shown that standard MODFLOW‐2005 code (as well as MODFLOW‐2000 and older versions) has significant convergence problems simulating such cases. The problems manifest themselves either in a total convergence failure or erroneous results. Alternatively, MODFLOW‐NWT, providing a good match to the presented discontinuous analytical solutions, appears to be a more reliable and appropriate code for simulating abrupt changes in water table elevations.     

Modelling the non‐linear hydrological behaviour of a small Mediterranean forested catchment

A progressive perceptual understanding approach was used to identify a model structure able to represent the non‐linear behaviour of the hydrological cycle in a small intermittent Mediterranean stream. The initial lumped model structure consisting of a series of four connected water tanks (LU3) progressed to a model with five tanks (LU4), and finally to a semidistributed model structure (SD4) in which spatial variability of the evapotranspiration according to the vegetation cover and to the local aspect was considered. In the final model structure, which gave the best fit (Nash–Sutcliffe efficiency index = 0·78), an additional tank representing the riparian zone was included (SD4‐R). Results showed that the abrupt changes of the riparian water table during summer and the formation of a perched water table during the transition from dry to wet conditions were the main mechanisms leading to the non‐linear hydrological behaviour. The transpiration process from the saturated zone and the spatial variability of evapotranspiration resulted in key factors successfully representing the annual water balance. The spatial and temporal validations carried out for each of the four model structures considered in this study supported the hypothesis adopted during the calibration process. Copyright © 2008 John Wiley & Sons, Ltd.     

How important is groundwater availability and stream perenniality to riparian and floodplain tree growth?

Riparian vegetation is important for stream functioning and as a major landscape feature. For many riparian plants, shallow groundwater is an important source of water, particularly in areas where rainfall is low, either annually or seasonally, and when extended dry conditions prevail for all or part of the year. The nature of tree water relationships is highly complex. Therefore, we used multiple lines of evidence to determine the water sources used by the dominant tree species Eucalyptus camaldulensis (river red gum), growing in riparian and floodplain areas with varying depth to groundwater and stream perenniality. Dendrometer bands were used to measure diel, seasonal, and annual patterns of tree water use and growth. Water stable isotopes (δ²H and δ¹⁸O) in plant xylem, soil water, and groundwater were measured to determine spatial and temporal patterns in plant water source use. Our results indicated riparian trees located on relatively shallow groundwater had greater growth rates, larger diel responses in stem diameter, and were less reactive to extended dry periods, than trees in areas of deep groundwater. These results were supported by isotope analysis that suggested all trees used groundwater when soil water stores were depleted at the end of the dry season, and this was most pronounced for trees with shallow groundwater. Trees may experience more frequent periods of water deficit stress and undergo reduced productivity in scenarios where water table accessibility is reduced, such as drawdown from groundwater pumping activities or periods of reduced rainfall recharge. The ability of trees to adapt to changing groundwater conditions may depend on the speed of change, the local hydrologic and soil conditions as well as the species involved. Our results suggest that E. camaldulesis growing at our study site is capable of utilizing groundwater even to depths >10 m, and stream perenniality is likely to be a useful indicator of riparian tree use of groundwater.     

The impact of wildfire on baseflow recession rates in California

The effect of wildfire on peak streamflow and annual water yield has been investigated empirically in numerous studies. The effect of wildfire on baseflow recession rates, in contrast, is not well documented. The objective of this paper was to quantify the effect of wildfire on baseflow recession rates in California for both individual watersheds and for all the study watersheds collectively. Two additional variables, antecedent groundwater storage and potential evapotranspiration, were also investigated for their effect on baseflow recession rates and postfire baseflow recession rate response. Differences between prefire and postfire baseflow recession rates were modeled statistically in 8 watersheds using a mixed statistical model that accounted for fixed and random effects. For the all‐watershed model, antecedent groundwater storage, potential evapotranspiration, and wildfire were each found to be significant controls on baseflow recession rates. Wildfire decreased baseflow recession rates 52.5% (37.6% to 66.0%), implying that postfire reductions in above‐ground vegetation (e.g., decreased interception, decreased evapotranspiration) were a stronger control on baseflow recession rate change than hydrophobicity. At an individual watershed scale, baseflow recession rate response to wildfire was found to be sensitive to intraannual differences in antecedent groundwater storage in 2 watersheds, with the effect of wildfire on baseflow recession rates being greater with lower levels of antecedent groundwater storage. Examination of burn severity for a subset of the study watersheds pointed to riparian zone burn severity as a potential primary control on postfire recession rate change. This study demonstrates that wildfire may have a substantial impact on fluxes to and from groundwater storages, altering the rate at which baseflow recedes.