Impact of water source dynamics on dissolved reactive phosphorus loadings in heterogenous karst agroecosystems with phosphatic limestones

Karst landscapes underlain with phosphatic limestones are now recognized to be an important contributor of fluvial phosphorus (P) to coastal waters. Specifically, karst agroecosystems may be a hotspot for dissolved reactive P (DRP) due to chronic over-application of organic and inorganic fertilizers that create legacy P accumulation in surface soils. Nevertheless, few studies have assessed the hydrologic controls on DRP transport in these systems at the watershed scale, which is the focus of this study. We analysed soil moisture, soil water extractable P, and storm event hydrologic and water quality data from a small heterogenous karst watershed (10.7 km²) in the Inner-Bluegrass Region of Central Kentucky, USA. Four storm events were sampled in winter, 2020 and were analysed for flow pathways using hydrograph recession analysis and water source connectivity using a tracer-based unmixing model. Based on hydrograph separation results, multiple linear regression analysis was performed to assess drivers of DRP concentrations and loadings. Soil water extractable P results showed stark vertical gradients with greater concentrations at both the surface and deeper soil zones, and minimum concentrations in the root zone. Results for the storm event analysis showed that water source connectivity provided superior prediction of DRP concentrations over the flow pathway analysis, which reflected the heterogeneity of karst maturity masking intermediate flow pathways. Findings from the MLR and loading analysis suggest waters sourced from the soil/epikarst produced significantly higher loadings compared with phreatic and precipitation water source in the three largest events, although concentrations fell between the phreatic (low) and precipitation (high) sources. Findings highlight variable activation of matrix-macropore exchange at different depths throughout the event. Collectively these results suggest existing models and approaches to assess karst hydrology need revision to improve management strategies in this critical landscape.     

Assessing urban rainfall-runoff response to stormwater management extent

We report an empirical analysis of the hydrologic response of three small, highly impervious urban watersheds to pulse rainfall events, to assess how traditional stormwater management (SWM) alters urban hydrographs. The watersheds vary in SWM coverage from 3% to 61% and in impervious cover from 45% to 67%. By selecting a set of storm events that involved a single rainfall pulse with >96% of total precipitation delivered in 60 min, we reduced the effect of differences between storms on hydrograph response to isolate characteristic responses attributable to watershed properties. Watershed-average radar rainfall data were used to generate local storm hyetographs for each event in each watershed, thus compensating for the extreme spatial and temporal heterogeneity of short-duration, intense rainfall events. By normalizing discharge values to the discharge peak and centring each hydrograph on the time of peak we were able to visualize the envelope of hydrographs for each group and to generate representative composite hydrographs for comparison across the three watersheds. Despite dramatic differences in the fraction of watershed area draining to SWM features across these three headwater tributaries, we did not find strong evidence that SWM causes significant attenuation of the hydrograph peak. Hydrograph response for the three watersheds is remarkably uniform despite contrasts in SWM, impervious cover and spatial patterns of land cover type. The primary difference in hydrograph response is observed on the recession limb of the hydrograph, and that change appears to be associated with higher storm-total runoff in the watersheds with more area draining to SWM. Our findings contribute more evidence to the work of previous authors suggesting that SWM is less effective at attenuating urban hydrographs than is commonly assumed. Our findings also are consistent with previous work concluding that percent impervious cover may have greater influence on runoff volume than percent SWM coverage.     

Hydrologic connectivity at the hillslope scale through intra-snowpack flow paths during snowmelt

The seasonal snowmelt period is a critical component of the hydrologic cycle for many mountainous areas. Changes in the timing and rate of snowmelt as a result of physical hydrologic flow paths, such as longitudinal intra-snowpack flow paths, can have strong implications on the partitioning of meltwater amongst streamflow, groundwater recharge, and soil moisture storage. However, intra-snowpack flow paths are highly spatially and temporally variable and thus difficult to observe. This study utilizes new methods to non-destructively observe spatio-temporal changes in the liquid water content of snow in combination with plot experiments to address the research question: What is the scale of influence that intra-snowpack flow paths have on the downslope movement of liquid water during snowmelt across an elevational gradient? This research took place in northern Colorado with study plots spanning from the rain-snow transition zone up to the high alpine. Results indicate an increasing scale of influence from intra-snowpack flow paths with elevation, showing higher hillslope connectivity producing larger intra-snowpack contributing areas for meltwater accumulation, quantified as the upslope contributing area required to produce observed changes in liquid water content from melt rate estimates. The total effective intra-snowpack contributing area of accumulating liquid water was found to be 17, 6, and 0 m² for the above tree line, near tree line, and below tree line plots, respectively. Dye tracer experiments show capillary and permeability barriers result in increased number and thickness of intra-snowpack flow paths at higher elevations. We additionally utilized aerial photogrammetry in combination with ground penetrating radar surveys to investigate the role of this hydrologic process at the small watershed scale. Results here indicate that intra-snowpack flow paths have influence beyond the plot scale, impacting the storage and transmission of liquid water within the snowpack at the small watershed scale.     

The assumption of uniform specific discharge: unsafe at any time?

Nearby catchments in the same landscape are often assumed to have similar specific discharge (runoff per unit catchment area). Five years of streamflow from 14 nested catchments in a 68 km² landscape was used to test this assumption, with the hypothesis that the spatial variability in specific discharge is smaller than the uncertainties in the measurement. The median spatial variability of specific discharge, defined as subcatchment deviation from the catchment outlet, was 33% at the daily scale. This declined to 24% at a monthly scale and 19% at an annual scale. These specific discharge differences are on the same order of magnitude as predicted for major land‐use conversions or a century of climate change. Spatial variability remained when considering uncertainties in specific discharge, and systematic seasonal patterns in specific discharge variation further provide confidence that these differences are more than just errors in the analysis of catchment area, rainfall variability or gauging. Assuming similar specific discharge in nearby catchments can thus lead to spurious conclusions about the effects of disturbance on hydrological and biogeochemical processes. Copyright © 2016 John Wiley & Sons, Ltd.     

Influences of vegetation disturbance on hydrogeomorphic response following wildfire

Quantifying the linkages between vegetation disturbance by fire and the changes in hydrologic processes leading to post‐fire erosional response remains a challenge. We measured the influence of fire severity, defined as vegetation disturbance (using a satellite‐derived vegetation disturbance index, VDI), landscape features that organize hydrologic flow pathways (relief and elongation ratios), and pre‐fire vegetation type on the probability of the occurrence of post‐fire gully rejuvenation (GR). We combined field surveys across 270 burned low‐order catchments (112 occurrences of GR) and geospatial analysis to generate a probabilistic model through logistic regression. VDI alone discriminated well between catchments where GR did and did not occur (area under the curve = 0.78, model accuracy = 0.72). The strong effect of vegetation disturbance on GR suggests that vegetation exerts a primary influence on the occurrence of infiltration excess run‐off and post‐fire erosion and that major gully erosion will not occur until fire consumes aboveground biomass. Other topographic and local factors also influenced GR response, including catchment elongation, per cent pre‐fire shrub, mid‐slope riparian vegetation, armoured headwaters, firehose effects, and concentration of severe burn in source areas. These factors highlight the need to consider vegetation effects in concert with local topography and site conditions to understand the propensity for flow accumulation leading to GR. We present a process‐based conceptual hydrologic model where vegetation loss from fire decreases rainfall attenuation and surface roughness, leading to accelerated flow accumulation and erosion; these effects are also influenced by interactions between fire severity and landscape structure. The VDI metric provides a continuous measure of vegetation disturbance and, when placed in a hydrologic context, may improve quantitative analysis of burned‐area susceptibility to erosive rainfall, hazard prediction, ecological effects of fire, landform evolution, and sensitivity to climate change. Copyright © 2015 John Wiley & Sons, Ltd.     

Hydrological impacts of seismic lines in the wetland‐dominated zone of thawing,discontinuous permafrost,Northwest Territories,Canada

Intensive seismic exploration in the Northwest Territories began in the late 1960s. Since that time, the legacy of seismic surveys – i.e. straight lines cutting through boreal forest and tundra – has remained visible throughout northern Canada and Alaska. The removal of trees and compaction of the ground surface alter the thermophysical properties of the active (i.e. seasonally thawed) layer to such an extent that the underlying permafrost seriously degrades or even disappears completely. Such a transformation along linear corridors that cut indiscriminately across different terrain types with contrasting hydrological functions has potentially serious implications to the redistribution of water and energy within and among landscape units with feedbacks to permafrost thaw, land cover change and run‐off generation. This paper characterizes the flow and storage of water and energy along a seismic cut line in the high boreal zone of discontinuous permafrost in order to improve the understanding of these processes, their interactions and hydrological implications. As such, this paper lays a conceptual foundation for the development of numerical models needed to predict the hydrological and thermal impact of seismic lines in this sensitive region. We used ground‐penetrating radar and multi‐year ground temperatures and water levels along a seismic line to estimate the degree of permafrost degradation below it. The seismic line studied extends from a permafrost‐free wetland (flat bog), over a permafrost body (peat plateau) and into another permafrost‐free wetland (channel fen). It was found that once thaw had lowered the permafrost table below the ground surface elevation of the flat bog and channel fen, the seismic line forms a hydrological connection between them. It was also shown that during the permafrost thaw process, seismic lines develop a perennially thawed layer (talik) between the overlying active layer and underlying permafrost and that the talik conveys water as a conduit throughout the year. The implications of such drainage through seismic lines and networks on basin drainage in peatland‐dominated regions with discontinuous permafrost are also discussed. Copyright © 2015 John Wiley & Sons, Ltd.     

Extreme precipitation time trends in Ontario, 1960–2010

The impact of climate change on the behaviour of intensity–duration–frequency curves is critical to the estimation of design storms, and thus to the safe design of drainage infrastructure. The present study develops a regional time trend methodology that detects the impact of climate change on extreme precipitation from 1960 to 2010. The regional time trend linear regression method is fitted to different durations of annual maximum precipitation intensities derived from multiple sites in Ontario, Canada. The results show the relationship between climate change and increased extreme precipitation in this province. The regional trend analysis demonstrates, under nonstationary conditions arising from climate change, that the intensity of extreme precipitation increased decennially between 1.25% for the 30‐min storm and 1.82% for the 24‐h storm. A comparison of the results with a regional Mann–Kendall test validates the found regional time‐trend results. The results are employed to extrapolate the intensity–duration–frequency curves temporally and spatially for future decades across the province. The results of the regional time trend assessment help with the establishment of new safety margins for infrastructure design in Ontario. Copyright © 2016 John Wiley & Sons, Ltd.     

Comparative analysis of baseflow characteristics of two Andean catchments,Ecuador

Baseflow in the Andes is commonly considered to be related with the release of water stored in páramos. Páramo is the predominant ecosystem above 3500 m a.s.l. and is characterized by a rainy and cold climate with low evapotranspiration. However, this baseflow concept is based on hydrological process studies in small Andean catchments of a few square kilometre with a homogeneous land cover. Middle‐sized Andean catchments, like the subcatchments of Tarqui and Yanuncay, Ecuador, are rarely homogeneous or uniformly covered by páramo. The objectives of this study are therefore to investigate baseflow characteristics in heterogeneous Andean catchments and to identify relationships between baseflow processes and physical characteristics such as storage and recharge. Hereby, the contribution to baseflow of páramo and other sources such as alluvial aquifers is quantified. This study uses nonlinear recession analysis, physically based filters and digital filters for comparison of baseflow of neighbouring but distinct subcatchments. The Yanuncay subcatchment shows a clearly different storage capacity and recession. The storage capacity of Yanuncay is 50% higher than for Tarqui because of its higher coverage of páramo. On the other hand, considerable storage capacity has also been found in the Tarqui subcatchment, which has a limited páramo area but a significant alluvial aquifer. It is shown that improved understanding of the specific baseflow characteristics such as storage and recharge and its relationships to the heterogeneity of the land cover in Andean catchments will lead to a better assessment of the water resources and give new insights for effective management actions. Copyright © 2014 John Wiley & Sons, Ltd.     

Spatial variation of wetlands and flux of dissolved organic carbon in boreal headwater streams

In order to investigate the relation between water chemistry and functional landscape elements, spatial data sets of characteristics for 68 small (0·2–1·5 km²) boreal forest catchments in western central Sweden were analysed in a geographical information system (GIS). The geographic data used were extracted from official topographic maps. Water sampled four times at different flow situations was analysed chemically. This paper focuses on one phenomenon that has an important influence on headwater quality in boreal, coniferous forest streams: generation and export of dissolved organic carbon (DOC). It is known that wetland cover (bogs and fens) in the catchment is a major source of DOC. In this study, a comparison was made between a large number of headwater catchments with varying spatial locations and areas of wetlands. How this variation, together with a number of other spatial variables, influences the DOC flux in the streamwater was analysed by statistical methods. There were significant, but not strong, correlations between the total percentages of wetland area and DOC flux measured at a medium flow situation, but not at high flow. Neither were there any significant correlations between the percentage of wetland area connected to streams, nor the percentage of wetland area within a zone 50 m from the stream and the DOC flux. There were, however, correlations between catchment mean slope and the DOC flux in all but one flow situations. This study showed that, considering geographical data retrieved from official sources, the topography of a catchment better explains the variation in DOC flux than the percentage and locations of distinct wetland areas. This emphasizes the need for high‐resolution elevation models accurate enough to reveal the sources of DOC found in headwater streams. Copyright © 2007 John Wiley & Sons, Ltd.