Annual Water Budget Research Articles

Extent of gross underestimation of precipitation in India

Abstract. The underestimation of precipitation (UoP) in the hilly and mountainous parts of South Asia is estimated by some studies to be as large as the observed precipitation (P). However, UoP has been analyzed to only a limited extent across India. To help bridge this gap, watershed-scale UoP was analyzed using various P datasets within a water imbalance analysis. Among these P datasets, the often-used Indian Meteorological Department (IMD) dataset is of primary interest. The gross UoP was identified by analyzing the extent of the imbalance in the annual water budget of watersheds corresponding to 242 river gauging stations for which quality-controlled data on catchment boundaries and streamflow are available. The water year (WY)-based volume of observed annual P was compared against the observed annual streamflow (R) and the satellite-based actual evapotranspiration (ET). Across many watersheds of both Northern and Peninsular India, spurious water imbalance scenarios (P≤R or P≪R+ET) were realized. It is shown that the management of water, such as groundwater extraction, reservoir storage and water diversion, is generally minimal compared to the annual P in such watersheds. It is also shown that annual changes in terrestrial water storage are minimal compared to the annual P in such watersheds. Assuming that data on R (and, to a lesser extent, ET) are reliable, it is concluded that UoP is very likely the cause of this imbalance. Inter-watershed groundwater flow (IGF) is assumed to be negligible. While the effect of IGF on R is unknown, examples are provided which show that IGF is unlikely to be the cause of the observed imbalance in certain watersheds. All 12 of the P datasets analyzed here suffer from UoP, but the extent of the UoP varies by dataset and region. The reanalysis-based datasets ERA5-Land and IMDAA are less affected by UoP than the IMD dataset. Based on the 30-year period of WY 1985–2014, P for the whole of India could be as much as 19 % (ERA5-Land) to 37 % (IMDAA) higher than that from the IMD, with substantial variability within years and river basins. The actual magnitude of UoP is speculated to be even greater. Moreover, trends seen in the IMD's P are not always present in ERA5-Land and IMDAA. Studies using IMD should exercise caution since UoP could lead to the misrepresentation of water budgets and long-term trends. Limitations of this study are discussed.

Wetland water balance estimation in an arid region using remote sensing technology and hydrological modelling

ABSTRACT Understanding the hydrological processes associated with wetland dynamics is fundamental to studying climate change impacts and the global water cycle. This study simulates the variation of water balance in the Ghorra Playa, Eastern Tunisia, over the period extending from September 2017 to August 2020. The playa is a seasonal habitat for a number of migratory bird species and its watershed is of vital agricultural importance. Water balance estimations were performed using two different approaches: remote sensing-based analyses (Sentinel-2B, Global Precipitation Measurements (GPM) and Famine Early Warning System Network (FEWS-NET)) and modelling founded on field data. The playa presents an average annual water balance of 1.1 million cubic metres. Water inflows come from direct rainfall, mostly in the fall and spring seasons. Groundwater flow into the playa significantly influenced the pattern of water flux (81%). The annual, seasonal and monthly water budget simulations show reasonably good agreement between the remote sensing-based analysis and the hydrological modelling exercise.

Effect of vegetation treatment and water stress on evapotranspiration in bioretention systems

Evapotranspiration is a key hydrological process for reducing stormwater runoff in bioretention systems, regardless of their physical configuration. Understanding the volumes of stormwater that can be returned to the atmosphere via evapotranspiration is, therefore, a key consideration in the design of any bioretention system. This study establishes the evapotranspiration dynamics of three common, structurally different, bioretention vegetation treatments (an Amenity Grass mix, and mono-cultures of Deschampsia cespitosa and Iris sibirica) compared with an un-vegetated control using lab-scale column experiments. Via continuous mass and moisture loss data, observed evapotranspiration rates were compared with those predicted by the FAO-56 Penman–Monteith model for five 14-day dry periods during Spring 2021, Summer 2021, and Spring 2022. Soil moisture reductions over the 14-day trials led to reduced rates of evapotranspiration. This necessitated the use of a soil moisture extraction function alongside a crop coefficient to represent actual evapotranspiration from FAO-56 Penman–Monteith reference evapotranspiration estimates. Crop coefficients (Kc) varied between 0.65 and 2.91, with a value of 1.0 identified as a recommended default value in the absence of treatment-specific empirical data. A continuous hydrological model with Kc=1.0 and a loading ratio of 10:1 showed that evapotranspiration could account for between 1 and 12% of the annual water budget for a bioretention system located in the UK and Ireland, increasing to a maximum of 35% when using the highest Kc observed (2.91).

Dynamics of estrogenic activity in an urban river receiving wastewater effluents: effect-based measurements with CALUX

Estrogenic substances (ES) in an urban river Zenne (BE) dominated by wastewater effluents were assessed over the course of one year. To measure the bioequivalent (BEQ) 17 β-estradiol (E2) concentrations of ES, the biological effect-based methodology - the Chemical-Activated LUciferase gene eXpression (CALUX) bioassay was used. Daily water discharges were collected from January 2015 to February 2016 at or near the sampling stations in the Brussels Capital Region. An annual water budget shows that approximately 50% of the Zenne River flow downstream is from wastewater effluent. The estrogenic activity and yearly average ES load in influents and effluents of wastewater treatment plants (WWTPs) located in the North and South, combined sewer overflows (CSOs) and the Zenne River, were assessed for upstream and downstream of two WWTPs of Brussels. Both WWTPs with activated sludge treatment remove more than 90% of the ES. The influent concentrations of ES at the South and North WWTPs ranged from 30-359 and 18-55 ng E2 eq./L, respectively. The effluent concentrations of ES ranged from 1.0-2.1 and 1.1-6.6 ng E2 eq./L at WWTP-S and -N, respectively. The yearly average ES loads were 0.05-0.14 and 0.39-1.5 g E2 eq./d for WWTP-S and -N, respectively. The temporal variation of E2-eq concentrations at the river stations Z3 and Z5 (upstream) ranged from 1 to 2 ng E2 eq./L, while the ES activity at sites Z9 and Z11 (downstream) varied from 2-17 ng E2 eq./L and from 1-8 ng/L ng E2 eq./L, respectively. The relative ES loads to the Zenne River are as follows: WWTPs (31%), CSOs (27%), upstream Zenne (15%), a missing source (14%), and local tributaries (13%). ES in the Zenne River behave in a pseudo-persistent manner because of continuous input from the WWTPs and slow degradation in the 18 km river stretch. The BEQ concentration of E2 exceeds the EU environmental quality standards (EQS) of 0.4 ng E2/L throughout the Zenne River.

Water Budget, Biological Water Use, and the Soil Hydrological Cycle across Typical Ecosystems of the Heihe River Basin

Quantification of the water budget of an arid inland river ecosystem is essential but still a challenge for the sustainable development of water resources. In situ observed data were used to analyze the monthly and annual water budgets and the soil hydrological cycle for six typical ecosystems in the Heihe River Basin (HRB). The two-source model was used to partition evapotranspiration (ET) into transpiration (T) and evaporation, after which the validated model was applied to quantitatively analyze the biological water use fraction [T/Ecosystem Water Supply (WS)] for different ecosystems. There were differences in the water budgets of the different ecosystems due to differences in climate, vegetation, soil, and external inputs. Precipitation in the HRB decreased from upstream to downstream, whereas there was a gradual increase in ET. External sources of water (e.g., natural runoff from upstream, irrigation in the middle reaches, and groundwater recharge in the lower reaches) to soil layers played an important role in regulating the water budgets of HRB ecosystems. Cropland obtained the maximum biological water use fraction (0.50), followed by Populus euphratica (0.49), alpine meadow (0.49), alpine swamp meadow (0.44), Tamarix ramosissima (0.42), and Kalidium foliatum (0.4). The soil water residence time (at a depth of 40 cm) varied from 14 d to 97 d (average of 60 d). The order of plant species in terms of soil water residence time was: K. foliatum (88 d) > T. ramosissima (72 d) > alpine meadow (68 d) > alpine swamp meadow (63 d) > cropland (53 d) > P. euphratica forest (20 d). Differences in the biological water use fraction and soil water residence time could be attributed to the characteristics of the water budget for each ecosystem. This study quantified the water budget, biological water use, and soil hydrological cycle across typical ecosystems in HRB, and can act as a reference for ecosystem management of the arid inland river basin.

Estimating the Water Budget of the Upper Blue Nile River Basin With Water and Energy Processes (WEP) Model

Understanding the spatial and temporal distributions and variations of basin water budget components is essential for effective water resources management. Due to a lack of basic hydro-meteorological information, the Upper Blue Nile River Basin (UBNRB) remains poorly understood in quantifying its hydrologic fluxes and associated dynamics. This study used a physically based distributed hydrologic model, WEP. We used multi-year land use information to better estimate the water budget components (evapotranspiration, runoff and storage) of the UBNRB. WEP simulation was validated at two main sections of the Upper Blue Nile river monthly from 1992 to 2014 (23 years). Results show that the basin stores a significant amount of water during the long rainy season (June to September) due to higher precipitation and limited evapotranspiration. However, it loses this storage through evapotranspiration during the dry season (October to February). The overall basin precipitation is 1,051 mm per year. Evapotranspiration accounts for 58% of the annual water budget, runoff is 25% and storage is 18%. The findings reported in this study can shed some light on understanding the UBNRB water budget dynamics and inform water management practitioners.

Dormant Season Vegetation Phenology and Eddy Fluxes in Native Tallgrass Prairies of the U.S. Southern Plains

Carbon dioxide (CO2) fluxes and evapotranspiration (ET) during the non-growing season can contribute significantly to the annual carbon and water budgets of agroecosystems. Comparative studies of vegetation phenology and the dynamics of CO2 fluxes and ET during the dormant season of native tallgrass prairies from different landscape positions under the same climatic regime are scarce. Thus, this study compared the dynamics of satellite-derived vegetation phenology (as captured by the enhanced vegetation index (EVI) and the normalized difference vegetation index (NDVI)) and eddy covariance (EC)-measured CO2 fluxes and ET in six differently managed native tallgrass prairie pastures during dormant seasons (November through March). During December–February, vegetation phenology (EVI and NDVI) and the dynamics of eddy fluxes were comparable across all pastures in most years. Large discrepancies in fluxes were observed during March (the time of the initiation of growth of dominant warm-season grasses) across years and pastures due to the influence of weather conditions and management practices. The results illustrated the interactive effects between prescribed spring burns and rainfall on vegetation phenology (i.e., positive and negative impacts of prescribed spring burns under non-drought and drought conditions, respectively). The EVI better tracked the phenology of tallgrass prairie during the dormant season than did NDVI. Similar EVI and NDVI values for the periods when flux magnitudes were different among pastures and years, most likely due to the satellite sensors’ inability to fully observe the presence of some cool-season C3 species under residues, necessitated a multi-level validation approach of using ground-truth observations of species composition, EC measurements, PhenoCam (digital) images, and finer-resolution satellite data to further validate the vegetation phenology derived from the Moderate Resolution Imaging Spectroradiometer (MODIS) during dormant seasons. This study provides novel insights into the dynamics of vegetation phenology, CO2 fluxes, and ET of tallgrass prairie during the dormant season in the U.S. Southern Great Plains.

Spatial patterns of soil salinity in the central Argentinean Dry Chaco

The Dry Chaco is a semi-arid ecoregion in South America that hosts one of the largest dry forests in the world, but expansion of dryland agriculture and cattle ranching led to gradual conversion of native vegetation to anthropogenic land cover. The potential impact of these newly established agricultural lands on the surrounding environment is of great concern. Local studies have shown that deforestation leads to changes in the soil-water balance, and can expedite groundwater rise and mobilization of water-soluble salts to the surface affecting plant growth and crop productivity. This study is a regional assessment of soil salinity and salinization processes in the central Argentinean Dry Chaco. It is based on extensive dataset of 492 surface and 142 subsurface samples taken along east-west transects across the Dry Chaco. Soil electrical conductivity (EC) was used as an indicator of salinity, and supplementary information on salt crusts was derived from Google Earth imagery. Subsurface salinity (i.e., measured at 100 cm depth) showed clear regional patterns in natural soil salinity that are related to the annual water budget and topography. Besides primary salinization due to arrested drainage and landscape stagnation, anthropogenic activities increased secondary salinization, especially in the agricultural areas with shallow groundwater tables and irrigated croplands. In addition, the study demonstrated that remotely-sensed vegetation indices such as the seasonal variation in the density of green vegetation are particularly suitable to monitor regional variation in soil salinity. Our results show that the extent of future dryland salinization in the Dry Chaco will mainly depend on whether areas prone to natural soil salinity are further protected from deforestation, and the magnitude and rate of groundwater rise after deforestation as conditioned by local climate and geomorphology. This better understanding of soil salinity patterns and how they are affected as a result of anthropogenic activities is important for the implementation of appropriate and effective measures to prevent severe salinization.

Improving the representation of forests in hydrological models

Forests play a critical role in the hydrologic cycle, impacting the surface and groundwater dynamics of watersheds through transpiration, interception, shading, and modification of the atmospheric boundary layer. It is therefore critical that forest dynamics are adequately represented in watershed models, such as the widely applied Soil and Water Assessment Tool (SWAT). SWAT’s default parameterization generally produces unrealistic forest growth predictions, which we address here through an improved representation of forest dynamics using species-specific re-parameterizations. We applied this methodology to the two dominant pine species in the southeastern U.S., loblolly pine (Pinus taeda L.) and slash pine (Pinus elliotti). Specifically, we replaced unrealistic parameter values related to tree growth with physically meaningful parameters derived from publicly available remote-sensing products, field measurements, published literature, and expert knowledge. Outputs of the default and re-parameterized models were compared at four pine plantation sites across a range of management, soil, and climate conditions. Results were validated against MODIS-derived leaf area index (LAI) and evapotranspiration (ET), as well as field observations of total biomass. The re-parameterized model outperformed the default model in simulating LAI, biomass accumulation, and ET at all sites. The two parametrizations also resulted in substantially different mean annual water budgets for all sites, with reductions in water yield ranging from 13 to 45% under the new parameterization, highlighting the importance of properly parameterizing forest dynamics in watershed models. Importantly, our re-parameterization methodology does not require alteration to the SWAT code, allowing it to be readily adapted and applied in ongoing and future watershed modeling studies.

A catchment water balance assessment of an abrupt shift in evapotranspiration at the Hubbard Brook Experimental Forest, New Hampshire, USA

AbstractSmall catchments have served as sentinels of forest ecosystem responses to changes in air quality and climate. The Hubbard Brook Experimental Forest in New Hampshire has been tracking catchment water budgets and their controls – meteorology and vegetation – since 1956. Water budgets in four reference catchments indicated an approximately 30% increase in the evapotranspiration (ET) as estimated by the difference between precipitation (P) and runoff (RO) starting in 2010 and continuing through 2019. We analyzed the annual water budgets, cumulative deviations of the daily P, RO and water budget residual (WBR = P − RO), potential ET (PET) and indicators of subsurface storage to gain greater insight into this shift in the water budgets. The PET and the subsurface storage indicators suggest that this change in WBR was primarily due to increasing ET. While multiple long‐term hydrological and micrometeorological data sets were used to detect and investigate this change in ET, additional measurements of groundwater storage and soil moisture would enable better estimation of ET within the catchment water balance. Increasing the breadth of long‐term measurements across small gauged catchments allows them to serve as more effective sentinels of substantial hydrologic changes like the ET increase that we observed.

Quantifying floodwater impacts on a lake water budget via volume-dependent transient stable isotope mass balance

Abstract. Isotope mass balance models have undergone significant developments in the last decade, demonstrating their utility for assessing the spatial and temporal variability in hydrological processes and revealing significant value for baseline assessment in remote and/or flood-affected settings where direct measurement of surface water fluxes to lakes (i.e. stream gauging) are difficult to perform. In this study, we demonstrate that isotopic mass balance modelling can be used to provide evidence of the relative importance of direct floodwater inputs and temporary subsurface storage of floodwater at ungauged lake systems. A volume-dependent transient isotopic mass balance model was developed for an artificial lake (named lake A) in southern Quebec (Canada). This lake typically receives substantial floodwater inputs during the spring freshet period as an ephemeral hydraulic connection with a 150 000 km2 large watershed is established. First-order water flux estimates to lake A allow for impacts of floodwater inputs to be highlighted within the annual water budget. The isotopic mass balance model has revealed that groundwater and surface water inputs account for 60 %–71 % and 39 %–28 % of the total annual water inputs to lake A, respectively, which demonstrates an inherent dependence of the lake on groundwater. However, when considering the potential temporary subsurface storage of floodwater, the partitioning between groundwater and surface water inputs tends to equalize, and the lake A water budget is found to be more resilient to groundwater quantity and quality changes. Our findings suggest not only that floodwater fluxes to lake A have an impact on its dynamics during springtime but also significantly influence its long-term water balance and help to inform, understand, and predict future water quality variations. From a global perspective, this knowledge is useful for establishing regional-scale management strategies for maintaining water quality at flood-affected lakes, for predicting the response of artificial recharge systems in such settings, and for mitigating impacts due to land use and climate changes.

Egypt’s water budget deficit and suggested mitigation policies for the Grand Ethiopian Renaissance Dam filling scenarios

The Nile River is a unique environmental system and essential water resource for its basin riparian nations. Population growth, changes in precipitation patterns, damming and usage rights disputes present extreme challenges in utilizing and managing the basin’s primary water resource. These stress factors are of particular concern for highly populated Egypt, the furthest downstream recipient of the Nile’s water flow. Previously, colonial agreements had granted Egypt and Sudan the majority of water use rights on the Nile without neighboring Ethiopia receiving any specific allocation. Today, Ethiopia plans to increase its energy production through its Nile-powered Grand Ethiopian Renaissance Dam (GERD). While the 74-billion cubic meter (BCM) dam presents promising development opportunities for Ethiopia, the Nile’s altered flow will increase the existing water deficit for Egypt—the quantification and mitigation of which are still largely unconstrained and under intense debate. To address this deficiency, we estimate that the median total annual water budget deficit for Egypt during the filling period, considering seepage into the fractured rocks below and around the GERD reservoir, as well as the intrinsic water deficit and assuming no possible mitigation efforts by Egyptian authorities, will be ∼31 BCM yr−1, which would surpass one third of Egypt’s current total water budget. Additionally, we provide a feasibility index for the different proposed solutions to mitigate the above deficit and assess their economic impact on the GDP per capita. Our results suggest that the unmet annual deficit during the filling period can be partially addressed by adjusting the Aswan High Dam (AHD) operation, expanding groundwater extraction and by adopting new policies for cultivation of crops. If no prompt mitigation is performed, the short-term three-year filling scenario would generate a deficit that is equivalent to losses to the present cultivated area by up to 72% resulting in a total loss of the agricultural GDP by $51 billion during the above-mentioned filling period. Such figures are equivalent to a decrease in the total national GDP per capita by ∼8%, augmenting existing unemployment rates by 11%, potentially leading to severe socioeconomic instability.

DRAINMOD-P: A Model for Simulating Phosphorus Dynamics and Transport in Drained Agricultural Lands: II. Model Testing

HighlightsDRAINMOD-P was tested using a dataset from a drained field with desiccation cracks.Surface and subsurface phosphorus losses were mainly in the particulate form.Surface runoff was a major pathway for phosphorus loss in this field.The model performance in predicting edge-of-field phosphorus loss is promising.Abstract. The recently developed phosphorus (P) model DRAINMOD-P was tested using a four-year dataset from a subsurface-drained field in northwest Ohio with significant potential for desiccation cracking or preferential flow. The model satisfactorily predicted subsurface drainage discharge, with a monthly Nash-Sutcliffe efficiency (NSE) of 0.59 and index of agreement (IOA) of 0.89. Lack of annual water budget closure was reported and was likely caused by uncertainty in measured surface runoff and/or modeling approaches representing macropore flow. More than 80% of predicted surface and subsurface P losses were in the particulate form. Surface runoff was the major pathway for P loss, contributing 78% of predicted total P (TP) load. On average, predicted macropore flow represented about 15% of drainage discharge and contributed 21% of DRP loss via subsurface drains. The performance of DRAINMOD-P in predicting monthly dissolved reactive P and TP losses through subsurface drains can be rated as poor (NSE = 0.33 and IOA = 0.60) and very good (NSE = 0.81 and IOA = 0.95), respectively. DRAINMOD-P demonstrated potential for simulating P fate and transport in drained cropland. More testing is needed to further examine newly incorporated hydrological and biogeochemical components of the model. Keywords: Agricultural drainage, Edge-of-field phosphorus load, Macropore flow, Phosphorus model, Sediment yield, Water quality modeling.

Effects of shifting snowmelt regimes on the hydrology of non-alpine temperate landscapes

Across the globe, increasing winter temperatures resulting from climate change are affecting how and when snow melts. Despite the significance of snowfall to the annual water budget in mid-latitude, non-alpine regions, there has been limited research into shifting snowmelt dynamics in these areas. This study examines snowmelt changes resulting from warmer winters in the Great Lakes region of North America. Using multiple metrics calculated from station temperature data, recent years (2003–2017) were categorized as “warm” or “cool”, then snowmelt and regional hydrologic patterns within those year types were examined. Systematic differences were observed, including warmer years having less and earlier snowmelt than cooler years. Those changes have, in turn, led to both lower and earlier spring peak flows in streams and decreased net groundwater recharge in the northern regions. Additionally, we show that differences between warm and cool year types become more significant along a north–south gradient; differences between warm and cool years are more pronounced in northern regions with regard to streamflow, net recharge, melt amount and timing. These results provide a framework to examine linked changes between snowmelt and hydrology across high-latitude temperate regions of the world.

Analysis of floodplain forest sensitivity to drought.

Floodplain forests are very complex, productive ecosystems, capable of storing huge amounts of soil carbon. With the increasing occurrence of extreme events, they are today among the most threatened ecosystems. Our study's main goal was to assess the productivity of a floodplain forest located at Lanžhot in the Czech Republic from two perspectives: carbon uptake (using an eddy covariance method) and stem radius variations (using dendrometers). We aimed to determine which conditions allow for high ecosystem production and what role drought plays in reducing such production potential. Additionally, we were interested to determine the relative soil water content threshold indicating the onset and duration of this event. We hypothesized that summer drought in 2018 had the most significant negative effects on the overall annual carbon and water budgets. In contrast with our original hypothesis, we found that an exceptionally warm spring in 2018 caused a positive gross primary production (GPP) and evapotranspiration (ET) anomaly that consequently led in 2018 to the highest seasonal total GPP and ET from all of the investigated years (2015–2018). The results showed ring-porous species to be the most drought resistant. Relative soil water content threshold of approximately 0.45 was determined as indicating the onset of drought stress.This article is part of the theme issue ‘Impacts of the 2018 severe drought and heatwave in Europe: from site to continental scale’.

DRAINMOD Simulation of macropore flow at subsurface drained agricultural fields: Model modification and field testing

Macropores are critical pathways through which water and pollutants can bypass the soil matrix and be rapidly transported to subsurface drains and freshwater bodies. We modified the DRAINMOD model to simulate macropore flow using a simple approach as part of developing the DRAINMOD-P model to simulate phosphorus dynamics in artificially drained agricultural lands. The Hagen-Poiseuille’s law was used to estimate the flow capacity of macropores. When ponding depths on the soil surface are greater than Kirkham’s depth, water is assumed to flow through macropores directly to tile drains without interaction with the soil matrix. In the modified model, macropore size is adjusted based on wet or dry conditions while connectivity is altered by tillage. The model was tested using a 4-year data set from a subsurface drained field in northwest Ohio. The soils at the field are classified as very poorly drained and are prone to desiccation cracking. The modified model predicted the daily and monthly subsurface drainage with average Nash-Sutcliffe efficiency (NSE) values of 0.48 and 0.59, respectively. The cumulative drainage over the 4-year simulation period was under-predicted by 8%. The new macropore component was able to capture about 75% of 60 peak drainage flow events. However, surface runoff was over-predicted for the entire study period. Annual water budgets using measured data (precipitation, subsurface drainage, and surface runoff) and model predictions (evapotranspiration, vertical seepage, and change in storage) were not balanced with an average annual imbalance of 6.4 cm. The lack of closure in the water balance suggests that errors may have occurred in field measurements, particularly, surface runoff. Overall, incorporating macropore flow into DRAINMOD improved predictions of daily drainage peaks and enabled the model to predict subsurface drainage flux contributed by macropore flow, which is critical for expanding DRAINMOD to simulate phosphorus transport in subsurface drained agricultural land.

Watershed and streambank erosion modeling in a coldwater stream using the GWLF-E model: application and evaluation

While many watershed models estimate overland flow-related soil loss, streambank erosion is often overlooked even though it can be the dominant source of sediment in a catchment. We used the enhanced generalized watershed loading functions (GWLF-E) model to simulate the water budget, field erosion from the landscape, and streambank erosion along a gradient of agricultural to urban land cover from 1997 to 2015 in the Indian Mill Creek watershed of Michigan, USA. We evaluated discharge simulations using in-streamflow measurements and evaluated streambank erosion simulations using field-collected erosion pin measurements from a previous study. Annual water budget results suggest the creek is primarily groundwater fed, but that a per-subbasin average of 6–15% of precipitation becomes runoff. Field erosion contributed a per-subbasin average of 0.5–2.5 Mg ha−1 year−1 of sediment, while streambank erosion accounted for 0.2–50.1% of the subbasins’ total sediment yields. Average lateral erosion rate of streambanks in subbasins ranged from 0.04 to 7.37 cm year−1, with four subbasins exceeding 1.0 cm year−1. Evaluation with field-collected discharge data suggest that GWLF-E may have been overestimating discharge in the subbasins. Evaluation with streambank erosion data found that GWLF-E may have difficulty capturing the complexity of erosion rates, including instances where sediment is deposited on streambanks from upstream sources, and that it underestimates streambank erosion in headwater catchments sometimes over an order of magnitude. Our findings are important for watershed managers to evaluate modeling approaches and highlight the importance of modeling both field and streambank erosion.

Increased water yield and altered water partitioning follow wildfire in a forested catchment in the western United States

AbstractAs wildfires in much of the western United States increase in size, frequency, and severity, understanding the impact of these fires on water yield from forested headwater basins is essential to successful management of water resources. The current study examines the changes in partitioning of the hydrologic cycle in the Mill Creek Basin that follow the Chippy Creek Fire in Montana, USA, due to alterations to the vegetative regime. The analysis utilizes remote sensing‐based vegetation indices and evapotranspiration, a model‐interpolated precipitation product, and discharge data to assess annual water budgets and vegetative regimes in the Mill Creek Basin. After being almost 90% burned in the Chippy Creek Fire, vegetation in the catchment has shifted from almost exclusively mixed conifer forest to sagebrush scrubs and grasses. This shift in vegetation is accompanied by abrupt shifts in partitioning of the water budget, resulting in an altered ecohydrologic regime. Post‐fire, evapotranspiration decreases annually by 250 mm (46%), and evaporative fraction decreases by 0.53. However, evapotranspiration product biases may overestimate this decrease from pre‐ to post‐fire. This decrease in evapotranspiration results in an annual increase in streamflow of 136 mm, a 21% increase in the runoff ratio, and a 140% increase in water yield. These changes to the water budget are consistent for 10 years post‐fire and show no trend towards pre‐fire values during the study period. Results will help inform planning and management of water resources downstream of forested catchments that have been impacted by wildfire.

Partitioning snowmelt and rainfall in the critical zone: effects of climate type and soil properties

Abstract. Streamflow generation and deep groundwater recharge may be vulnerable to loss of snow, making it important to quantify how snowmelt is partitioned between soil storage, deep drainage, evapotranspiration, and runoff. Based on previous findings, we hypothesize that snowmelt produces greater streamflow and deep drainage than rainfall and that this effect is greatest in dry climates. To test this hypothesis we examine how snowmelt and rainfall partitioning vary with climate and soil properties using a physically based variably saturated subsurface flow model, HYDRUS-1D. We developed model experiments using observed climate from mountain regions and artificial climate inputs that convert all precipitation to rain, and then evaluated how climate variability affects partitioning in soils with different hydraulic properties and depths. Results indicate that event-scale runoff is higher for snowmelt than for rainfall due to higher antecedent moisture and input rates in both wet and dry climates. Annual runoff also increases with snowmelt fraction, whereas deep drainage is not correlated with snowmelt fraction. Deep drainage is less affected by changes from snowmelt to rainfall because it is controlled by deep soil moisture changes over longer timescales. Soil texture modifies daily wetting and drying patterns but has limited effect on annual water budget partitioning, whereas increases in soil depth lead to lower runoff and greater deep drainage. Overall these results indicate that runoff may be substantially reduced with seasonal snowpack decline in all climates, whereas the effects of snowpack decline on deep drainage are less consistent. These mechanisms help explain recent observations of streamflow sensitivity to changing snowpack and highlight the importance of developing strategies to plan for changes in water budgets in areas most at risk for shifts from snow to rain.

Energy and water mass balance of Lake Untersee and its perennial ice cover, East Antarctica

AbstractLake Untersee is one of the largest perennially ice-covered lakes in Dronning Maud Land. We investigated the energy and water mass balance of Lake Untersee to understand its state of equilibrium. The thickness of the ice cover is strongly correlated with sublimation rates; variations in sublimation rates across the ice cover are largely determined by wind-driven turbulent heat fluxes and the number of snow-covered days. Lake extent and water level have remained stable for the past 20 years, indicating that the water mass balance is in equilibrium. The lake is damned by the Anuchin Glacier and mass balance calculation suggest that subaqueous melting of terminus ice contributes 40–45% of the annual water budget; since there is no evidence of streams flowing into the lake, the lake must be connected to a groundwater system that contributes 55–60% in order to maintain the lake budget in balance. The groundwater likely flows at a rate of ~8.8 × 10−2m3s−1, a reasonable estimate given the range of subglacial water flux in the region. The fate of its well-sealed ice cover is likely tied to changes in wind regime, whereas changes in water budget are more closely linked to the response of surrounding glaciers to climate change.