Divergence and flow direction as indicators of subsurface heterogeneity and stage‐dependent storage in alluvial floodplains

Assuming homogeneity in alluvial aquifers is convenient, but limits our ability to accurately predict stream‐aquifer interactions. Research is needed on (i) identifying the presence of focused, as opposed to diffuse, groundwater discharge/recharge to streams and (ii) the magnitude and role of large‐scale bank and transient storage in alluvial floodplains relative to changes in stream stage. The objective of this research was to document and quantify the effect of stage‐dependent aquifer heterogeneity and bank storage relative to changes in stream stage using groundwater flow divergence and direction. Monitoring was performed in alluvial floodplains adjacent to the Barren Fork Creek and Honey Creek in northeastern Oklahoma. Based on results from subsurface electrical resistivity mapping, observation wells were installed in high and low electrical resistivity subsoils. Water levels in the wells were recorded real time using pressure transducers (August to October 2009). Divergence was used to quantify heterogeneity (i.e. variation in hydraulic conductivity, porosity, and/or aquifer thickness), and flow direction was used to assess the potential for large‐scale (100 m) bank or transient storage. Areas of localized heterogeneity appeared to act as divergence zones allowing stream water to quickly enter the groundwater system, or as flow convergence zones draining a large groundwater area. Maximum divergence or convergence occurred with maximum rates of change in flow rates or stream stage. Flow directions in the groundwater changed considerably between base and high flows, suggesting that the floodplains acted as large‐scale bank storage zones, rapidly storing and releasing water during passage of a storm hydrograph. During storm events at both sites, the average groundwater direction changed by at least 90° from the average groundwater direction during baseflow. Aquifer heterogeneity in floodplains yields hyporheic flows that are more responsive and spatially and temporally complex than would be expected compared to more common assumptions of homogeneity. Copyright © 2012 John Wiley & Sons, Ltd.     

Exploring the interaction of surface roughness and slope gradient in controlling rates of soil loss from sloping farmland on the Loess Plateau of China

Surface roughness and slope gradient are two important factors influencing soil erosion. The objective of this study was to investigate the interaction of surface roughness and slope gradient in controlling soil loss from sloping farmland due to water erosion on the Loess Plateau, China. Following the surface features of sloping farmland in the plateau region, we manually prepared rough surfaces using four tillage practices (contour drilling, artificial digging, manual hoeing, and contour plowing), with a smooth surface as the control measure. Five slope gradients (3°, 5°, 10°, 15°, and 20°) and two rainfall intensities (60 and 90 mm/hr) were considered in the artificial rainfall simulation experiment. The results showed that the runoff volume and sediment yield increased with increasing slope gradient under the same tillage treatment. At gentle slope gradients (e.g., 3° and 5°), the increase in surface roughness prevented the runoff and sediment production, that is, the surface roughness reduced the positive effect of slope gradient on the runoff volume and sediment yield to a certain extent. At steep slope gradients, however, the enhancing effect of slope gradient on soil erosion gradually increased and surpassed the reduction effect of surface roughness. This study reveals the existence of a critical slope gradient that influences the interaction of surface roughness and slope gradient in controlling soil erosion on sloping farmland. If the slope gradient is equal to or less than the critical value, an increase in surface roughness would decrease soil erosion. Otherwise, the increase in surface roughness would be ineffective for preventing soil erosion. The critical slope gradient would be smaller under higher rainfall intensity. These findings are helpful for us to understand the process of soil erosion and relevant for supporting soil and water conservation in the Loess Plateau region of China.     

Dissolved organic matter in surface runoff in the Loess Plateau of China: The role of rainfall events and land-use

Exploring the chemical characterization of dissolved organic matter (DOM) is important for understanding the fate of laterally transported organic matter in watersheds. We hypothesized that differences in water-extractable organic matter (WEOM) in soils of varying land uses and rainfall events may significantly affect the quality and the quantity of stream DOM. To test our hypotheses, characteristics of rainfall-runoff DOM and WEOM of source materials (topsoil from different land uses and gullies, as well as typical vegetation) were investigated at two adjacent catchments in the Loess Plateau of China, using ultraviolet–visible absorbance and excitation emission matrix fluorescence with parallel factor analysis (PARAFAC). Results indicated that land-use types may significantly affect the chemical composition of soil WEOM, including its aromaticity, molecular weight, and degree of humification. The PARAFAC analysis demonstrated that the soils and stream water were dominated by terrestrial/allochthonous humic-like substances and microbial transformable humic-like fluorophores. Shifts in the fluorescence properties of stream DOM suggested a pronounced change in the relative proportion of allochthonous versus autochthonous material under different rainfall patterns and land uses. For example, high proportions of forestland could provide more allochthonous DOM input. This study highlights the relevance of soils and hydrological dynamics on the composition and fluxes of DOM issuing from watersheds. The composition of DOM in soils was influenced by land-use type. Precipitation patterns influenced the proportion of terrestrial versus microbial origins of DOM in surface runoff. Contributions of allochthonous, terrestrially derived DOM inputs were highest from forested landscapes.     

Quantitative estimation of submarine groundwater discharge using airborne thermal infrared data acquired at two different tidal heights

Submarine groundwater discharge (SGD) plays an important role in coastal biogeochemical processes and hydrological cycles, particularly off volcanic islands in oligotrophic oceans. However, the spatial and temporal variations of SGD are still poorly understood owing to difficulty in taking rapid SGD measurements over a large scale. In this study, we used four airborne thermal infrared surveys (twice each during high and low tides) to quantify the spatiotemporal variations of SGD over the entire coast of Jeju Island, Korea. On the basis of an analytical model, we found a linear positive correlation between the thermal anomaly and squares of the groundwater discharge velocity and a negative exponential correlation between the anomaly and water depth (including tide height and bathymetry). We then derived a new equation for quantitatively estimating the SGD flow rates from thermal anomalies acquired at two different tide heights. The proposed method was validated with the measured SGD flow rates using a current meter at Gongcheonpo Beach. We believe that the method can be effectively applied for rapid estimation of SGD over coastal areas, where fresh groundwater discharge is significant, using airborne thermal infrared surveys.     

Hydrologic responses to land cover change: the case of Jedeb mesoscale catchment,Abay/Upper Blue Nile basin,Ethiopia

The objective of this study was to quantify the impacts of land use/land cover (LULC) change on the hydrology of the Jedeb, an agricultural dominated mesoscale catchment, in the Abay/Upper Blue Nile basin, Ethiopia. Two methods have been used. First, the trends of certain daily flow variability parameters were evaluated to detect statistical significance of the change of the hydrologic response. Second, a conceptual monthly hydrological model was used to detect changes in the model parameters over different periods to infer LULC change. The results from the statistical analysis of the daily flows between 1973 and 2010 reveal a significant change in the response of the catchment. Peak flow is enhanced, i.e. response appears to be flashier. There is a significant increase in the rise and fall rates of the flow hydrograph, as well as the number of low‐flow pulses below a threshold level. The discharge pulses show a declining duration with time. The model result depicts a change in model parameters over different periods, which could be attributed to an LULC change. The model parameters representing soil moisture conditions indicated a gradual decreasing trend, implying limited storage capacity likely attributed to increasing agricultural farming practices in the catchment. This resulted in more surface runoff and less infiltration into the soil layers. The results of the monthly flow duration curve analysis indicated large changes of the flow regime. The high flow has increased by 45% between the 1990s and 2000s, whereas the reduction in low flows was larger: a 15% decrease between 1970s and 1980s, 39% between 1980s and 1990s and up to 71% between 1990s and 2000s. These results, could guide informed catchment management practices to reduce surface runoff and augment soil moisture level in the Jedeb catchment. Copyright © 2013 John Wiley & Sons, Ltd.     

Alteration of hydrological processes and streamflow with juniper (Juniperus virginiana) encroachment in a mesic grassland catchment

To understand the effect of woody plant encroachment on hydrological processes of mesic grasslands, we quantified infiltration capacity in situ, the temporal changes in soil water storage, and streamflow of a grassland catchment and a catchment heavily encroached by juniper (Juniperus virginiana, eastern redcedar) in previously cultivated, non‐karst substrate grasslands in north‐central Oklahoma for 3 years. The initial and steady‐state infiltration rates under the juniper canopy were nearly triple to that of the grassland catchment and were intermediate in the intercanopy spaces within the encroached catchment. Soil water content and soil water storage on the encroached catchment were generally lower than on the grassland catchment, especially when preceding the seasons of peak rainfall in spring and fall. Frequency and magnitude of streamflow events were reduced in the encroached catchment. Annual runoff coefficients for the encroached catchment averaged 2.1%, in contrast to 10.6% for the grassland catchment. Annual streamflow duration ranged from 80 to 250 h for the encroached catchment compared with 600 to 800 h for the grassland catchment. Our results showed that the encroachment of juniper into previously cultivated mesic grasslands fundamentally alters catchment hydrological function. Rapid transformation of mesic grassland to a woodland state with juniper encroachment, if not confined, has the potential to drastically reduce soil water, streamflow and flow duration of ephemeral streams in the Southern Great Plains. Copyright © 2013 John Wiley & Sons, Ltd.     

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

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

Dynamic river–groundwater exchange in the presence of a saline,semi‐confined aquifer

Understanding groundwater–surface water exchange in river banks is crucial for effective water management and a range of scientific disciplines. While there has been much research on bank storage, many studies assume idealized aquifer systems. This paper presents a field‐based study of the Tambo Catchment (southeast Australia) where the Tambo River interacts with both an unconfined aquifer containing relatively young and fresh groundwater (<500 μS/cm and <100 years old) and a semi‐confined artesian aquifer containing old and saline groundwater (electrical conductivity > 2500 μS/cm and >10 000 years old). Continuous groundwater elevation and electrical conductivity monitoring within the different aquifers and the river suggest that the degree of mixing between the two aquifers and the river varies significantly in response to changing hydrological conditions. Numerical modelling using MODFLOW and the solute transport package MT3DMS indicates that saline water in the river bank moves away from the river during flooding as hydraulic gradients reverse. This water then returns during flood recession as baseflow hydraulic gradients are re‐established. Modelling also indicates that the concentration of a simulated conservative groundwater solute can increase for up to ~34 days at distances of 20 and 40 m from the river in response to flood events approximately 10 m in height. For the same flood event, simulated solute concentrations within 10 m of the river increase for only ~15 days as the infiltrating low‐salinity river water drives groundwater dilution. Average groundwater fluxes to the river stretch estimated using Darcy's law were 7 m³/m/day compared with 26 and 3 m³/m/day for the same periods via mass balance using Radon (²²²Rn) and chloride (Cl), respectively. The study shows that by coupling numerical modelling with continuous groundwater–surface water monitoring, the transient nature of bank storage can be evaluated, leading to a better understanding of the hydrological system and better interpretation of hydrochemical data. Copyright © 2015 John Wiley & Sons, Ltd.     

The exponential decline in saturated hydraulic conductivity with depth: a novel method for exploring its effect on water flow paths and transit time distribution

The strong vertical gradient in soil and subsoil saturated hydraulic conductivity is characteristic feature of the hydrology of catchments. Despite the potential importance of these strong gradients, they have proven difficult to model using robust physically based schemes. This has hampered the testing of hypotheses about the implications of such vertical gradients for subsurface flow paths, residence times and transit time distribution. Here we present a general semi‐analytical solution for the simulation of 2D steady‐state saturated‐unsaturated flow in hillslopes with saturated hydraulic conductivity that declines exponentially with depth. The grid‐free solution satisfies mass balance exactly over the entire saturated and unsaturated zones. The new method provides continuous solutions for head, flow and velocity in both saturated and unsaturated zones without any interpolation process as is common in discrete numerical schemes. This solution efficiently generates flow pathlines and transit time distributions in hillslopes with the assumption of depth‐varying saturated hydraulic conductivity. The model outputs reveal the pronounced effect that changing the strength of the exponential decline in saturated hydraulic conductivity has on the flow pathlines, residence time and transit time distribution. This new steady‐state model may be useful to others for posing hypotheses about how different depth functions for hydraulic conductivity influence catchment hydrological response. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessing the ‘two water worlds’ hypothesis and water sources for native and exotic evergreen species in south‐central Chile

Recent studies using water‐stable isotopes (δ¹⁸O and δ²H) have suggested an ecohydrological separation of water flowing to streams or recharging groundwater and water used by trees, known as the ‘two water worlds’ (TWW) hypothesis. In this study, we measured water isotopic composition in precipitation [open field and throughfall, i.e. local meteoric water line (LMWL)] and the mobile water compartment (i.e. stream and soil solution), bulk soil water and xylem water over a period of 1.5 years in two headwater catchments: NF, covered with old growth native evergreen forest (Aetoxicon punctatum, Laureliopsis philippiana and Eucriphya cordifolia), and EP, covered with 4 and 16‐year‐old Eucalyptus nitens stands. Our results show that precipitation, stream and soil solution plot approximately along the LMWL, while xylem waters from all studied tree species plot below the LMWL, supporting the TWW hypothesis. However, we also found evidence of ecohydrological connectivity during the wet season, likely controlled by the amount of antecedent precipitation. These observations hold for all investigated tree species. On both sites, a different precipitation source for stream and xylem water was observed. However, in EP, bulk soil showed a similar precipitation source as xylem water from both E. nitens stands. This suggests that E. nitens may use water that is recharging the bulk soil compartment. We conclude that under a rainy temperate climate, the TWW hypothesis is temporal and does not apply during wet seasons. Copyright © 2016 John Wiley & Sons, Ltd.