Nonlinear baseflow recession analysis in watersheds with intermittent streamflow

Abstract

Discharge in most rivers consists mainly of baseflow exfiltrating from shallow groundwater reservoirs, while surface or other direct flows cease soon after rain storms or snowmelt. Analysis of observed baseflow recessions of two rivers in Turkey with intermittent flows and different geographical and climatic characteristics yielded nonlinear storage–outflow relationships of the highly seasonal aquifers. Baseflow separation was carried out using a nonlinear reservoir algorithm. Baseflow seasonality is related to the hydro-climatic conditions influencing groundwater recharge and evapotranspiration of groundwater. As intermittent streams generally have zero flows in the dry season, calibration of recession parameters is in many cases a complicated task.

Citation Aksoy, H. & Wittenberg, H. (2011) Nonlinear baseflow recession analysis in watersheds with intermittent streamflow. Hydrol. Sci. J. 56(2), 226–237.     

Focused Groundwater Controlled Feedbacks into the Hyporheic Zone During Baseflow Recession

Groundwater surface water interaction in the hyporheic zone remains an important challenge for water resources management and ecosystem restoration. In heterogeneous stratified glacial sediments, reach‐scale environments contain an uneven distribution of focused groundwater flow occurring simultaneously with diffusely discharging groundwater. This results in a variation of stream‐aquifer interactions, where focused flow systems are able to temporally dominate exchange processes. The research presented here investigates the direct and indirect influences focused groundwater discharge exerts on the hyporheic zone during baseflow recession. Field results demonstrate that as diffuse sources of groundwater deplete during baseflow recession, focused groundwater discharge remains constant. During baseflow recession the hyporheic zone is unable to expand, while the high nitrate concentration from focused discharge changes the chemistry of the stream. The final result is a higher concentration of nitrate in the hyporheic zone as this altered surface water infiltrates into the subsurface. This indirect coupling of focused groundwater discharge and the hyporheic zone is unaccounted for in hyporheic studies at this time. Results indicate important implications for the potential reduction of agricultural degradation of water quality.     

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.     

Using periodic hydrologic and geochemical sampling with limited continuous monitoring to characterize remote karst aquifers in the Kaweah River Basin,California, USA

Hydrogeologic field work in remote settings is often challenging: assessing spring behaviour and aquifer characteristics can be expensive in both time commitment and resources needed to assess these systems. In this study, we document the hydrology and geochemistry of 47 perennial karst springs in the Kaweah River, a mountain river basin in the Sierra Nevada, California. After preliminary hydrogeochemical characterization and grouping, selected springs were continuously monitored to further assess aquifer characteristics in each group. Later, in areas without previous dye‐tracing work, traces were conducted to establish connections between large sinking streams and springs. The springs have a wide range of inter‐spring and intra‐spring variability in discharge and geochemistry. We assessed this variability by performing statistical comparisons with spring chemistry and principal components analysis of all measured variables. Results show that springs can be divided into two distinct groups: high elevation springs of the Mineral King Valley and lower elevation springs throughout the rest of the basin. Continuous discharge, temperature and specific conductivity data from four springs (two from each group) were then used to characterize the hydrograph recession behaviour of springs in each group. Both groups showed statistically similar baseflow recession slopes, suggesting that both groups contain baseflow storage compartments with similar hydrogeologic properties. The biggest difference between each group is the variability in amount of water remaining in the aquifer during baseflow conditions. High elevation springs have lower baseflow discharges, relative to peak flow, than lower elevation springs, despite the fact that more precipitation falls at higher elevation. This is likely caused by differences in the amount of soil and epikarst storage, which are related to recent geomorphic events: high elevation aquifers were glaciated as recent as 41 thousand years ago (kya), while there is no evidence that low elevation aquifers were glaciated. As a result, lower elevations have developed thicker soils, weathered bedrock and epikarst. Copyright © 2016 John Wiley & Sons, Ltd.     

Groundwater flow and storage processes in an inactive rock glacier

Groundwater flow through coarse blocky landforms contributes to streamflow in mountain watersheds, yet its role in the alpine hydrologic cycle has received relatively little attention. This study examines the internal structure and hydrogeological characteristics of an inactive rock glacier in the Canadian Rockies using geophysical imaging techniques, analysis of the discharge hydrograph of the spring draining the rock glacier, and chemical and stable isotopic compositions of source waters. The results show that the coarse blocky sediments forming the rock glacier allow the rapid infiltration of snowmelt and rain water to an unconfined aquifer above the bedrock surface. The water flowing through the aquifer is eventually routed via an internal channel parallel to the front of the rock glacier to a spring, which provides baseflow to a headwater stream designated as a critical habitat for an at‐risk cold‐water fish species. Discharge from the rock glacier spring contributes up to 50% of basin streamflow during summer baseflow periods and up to 100% of basin streamflow over winter, despite draining less than 20% of the watershed area. The rock glacier contains patches of ground ice even though it may have been inactive for thousands of years, suggesting the resiliency of the ground thermal regime under a warming climate.     

Effects of subsurface drainage tiles on streamflow in Iowa agricultural watersheds: Exploratory hydrograph analysis

Flow from artificial subsurface (tile) drainage systems may be contributing to increasing baseflow in Midwestern rivers and increased losses of nitrate‐nitrogen. Standard hydrograph analysis techniques were applied to model simulation output and field monitoring from tile‐drained landscapes to explore how flow from drainage tiles affects stream baseflow and streamflow recession characteristics. DRAINMOD was used to simulate hydrologic response from drained (24 m tile spacing) and undrained agricultural systems. Hydrograph analysis was conducted using programs PART and RECESS. Field monitoring data were obtained from several monitoring sites in Iowa typical of heavily drained and less‐drained regions. Results indicate that flow from tile drainage primarily affects the baseflow portion of a hydrograph, increasing annual baseflow in streams with seasonal increases primarily occurring in the late spring and early summer months. Master recession curves from tile‐drained watersheds appear to be more linear than less‐tiled watersheds although comparative results of the recession index k were inconsistent. Considering the magnitude of non‐point source pollutant loads coming from tile‐drained landscapes, it is critical that more in‐depth research and analysis be done to assess the effects of tile drainage on watershed hydrology if water quality solutions are to be properly evaluated. Copyright © 2008 John Wiley & Sons, Ltd.     

Quantifying the role of karstic groundwater in a snowmelt-dominated hydrologic system

River discharge in mountainous regions of the world is often dominated by snowmelt, but base flows are sustained primarily by groundwater storage and discharge. Although numerous recent studies have focused on base-flow discharge in mountain systems, almost no work has explicitly investigated the role of karst groundwater in these systems across a full range of flow conditions. We directly measured groundwater discharge from 48 karst springs in the Kaweah River and its five forks in the Sierra Nevada mountains, California, United States. Relationships between spring and river discharge showed that karst aquifers and springs provide significant storage and delayed discharge to the river. Regression models showed that, of all potential seasonal groundwater storage compartments in the river basin, the area of karst (0.1–4.4%) present provides the best explanation of base-flow recession in each fork of the Kaweah River (directly measured contributions from karst springs ranged from 3.5 to 16% during high-flow to 20 to 65% during base-flow conditions). These results show that, even in settings where karst represents a small portion of basin area, it may play an over-sized role in seasonal storage and water resources in mountain systems. Karst aquifers are the single most important non-snow storage component in the Kaweah River basin, and likely provide similar water storage capacities and higher base flows in other mountain river systems with karst when compared with systems without karst.     

A method to assess annual average renewable groundwater reserves for large regions in Spain

Abstract

This paper proposes a method for assessing the renewable groundwater reserves of large regions for an average year, based on the integration of the recession curves for their basins' springs or the natural baseflow of their rivers. In this method, the hydrodynamic volumes (or renewable reserves), were estimated from the baseflow. It was assumed that the flow was the same as the natural recharge, and the recession coefficients were derived from the hydrogeological parameters and geometric characteristics of the aquifers and adjusted to fit the recession curves at gauging stations. The method was applied to all the aquifers of Spain, which have a total renewable groundwater reserve of 86 118 hm³—four times the mean annual recharge. However, the spatial distribution of these reserves is highly variable; 18.6% of the country's aquifers contain 94.7% of the entire reserve.

Citation Sanz, E. & Recio, B. (2010) A method to assess annual average renewable groundwater reserves for large regions in Spain. Hydrol. Sci. J. 56(1), 99–107.     

Estimation of baseflow residence times in watersheds from the runoff hydrograph recession: method and application in the Neversink watershed,Catskill Mountains,New York

A method for estimation of mean baseflow residence time in watersheds from hydrograph runoff recession characteristics was developed. Runoff recession characteristics were computed for the period 1993–96 in the 2 km² Winnisook watershed, Catskill Mountains, southeastern New York, and were used to derive mean values of subsurface hydraulic conductivity and the storage coefficient. These values were then used to estimate the mean baseflow residence time from an expression of the soil contact time, based on watershed soil and topographic characteristics. For comparison, mean baseflow residence times were calculated for the same period of time through the traditional convolution integral approach, which relates rainfall δ¹⁸O to δ¹⁸O values in streamflow. Our computed mean baseflow residence time was 9 months by both methods. These results indicate that baseflow residence time can be calculated accurately using recession analysis, and the method is less expensive than using environmental and/or artificial tracers. Published in 2002 by John Wiley & Sons, Ltd.     

Modification of the SWAT model to simulate regional groundwater flow using a multicell aquifer

The soil and water assessment tool (SWAT) has been widely used and thoroughly tested in many places in the world. The application of the SWAT model has pointed out that 2 of the major weaknesses of SWAT are related to the nonspatial reference of the hydrologic response unit concept and to the simplified groundwater concept, which contribute to its low performance in baseflow simulation and its inability to simulate regional groundwater flow. This study modified the groundwater module of SWAT to overcome the above limitations. The modified groundwater module has 2 aquifers. The local aquifer, which is the shallow aquifer in the original SWAT, represents a local groundwater flow system. The regional aquifer, which replaces the deep aquifer of the original SWAT, represents intermediate and regional groundwater flow systems. Groundwater recharge is partitioned into local and regional aquifer recharges. The regional aquifer is represented by a multicell aquifer (MCA) model. The regional aquifer is discretized into cells using the Thiessen polygon method, where centres of the cells are locations of groundwater observation wells. Groundwater flow between cells is modelled using Darcy's law. Return flow from cell to stream is conceptualized using a non‐linear storage–discharge relationship. The SWAT model with the modified aquifer module, the so‐called SWAT‐MCA, was tested in 2 basins (Wipperau and Neetze) with porous aquifers in a lowland area in Lower Saxony, Germany. Results from the Wipperau basin show that the SWAT‐MCA model is able (a) to simulate baseflow in a lowland area (where baseflow is a dominant source of streamflow) better than the original model and (b) to simulate regional groundwater flow, shown by the simulated groundwater levels in cells, quite well.     

Normal-mode analysis of the structure of baseflow-recession curves

It has been observed for several years that semi-logarithmic plots of the baseflow recession for many streams which partially penetrate aquifers in the United Kingdom are not single straight-line plots. Simple modelling studies indicate that this behaviour, for typical values of aquifer parameters, occurs for even the simplest of groundwater systems. In general such studies show that the baseflow component arising from even a structurally uncomplicated aquifer should be viewed as a superposition of many, distinct exponential terms. The origin of these terms does not necessarily lie in a varied hydrogeological structure, but can arise purely from the dynamics of groundwater flow. Such a form for the baseflow recession can be used to perform empirical fits to observed data. Moreover, where more than one exponential term is contributing to the total baseflow, great care must be taken in the interpretation of associated data.     

Computer analysis of regional groundwater flow and boundary conditions in the basin of Mexico

Mexico City is situated in the Valley of Mexico on the extensive lacustrine clays that overlay highly productive aquifers of both volcanic and sedimentary origin. The Valley is closed by volcanic mountains. The natural hydraulic boundary conditions associated withe mountain ranges and their relationship to the important aquifers were studied using a two-dimensional, steady-state finite-element model in cross section. Four cross sections were analysed under hydrologic conditions existing prior to the large scale pumping of the aquifers. Factors such as bulk hydraulic conductivities and regional infiltration rates were obtained from field observations and the literature to assess location of the associated groundwater divides, and the water-table in the mountains. The modeled flow patterns are consistent with the historical hydrologic records piezometric characteristics and observed surface features of the groundwater in the Basin of Mexico. From the modeling results, the groundwater recharge in the mountains is 30–50% of the mean average precipitation. Higher and lower rates result in a flow regime that is not compatible with field observations. In general the location of the divides in the mountains is displaced towards the Valley of Mexico, which influences the groundwater budget of the Valley. The water table in places is several hundred metres below ground surface, in accordance with field observations of a very thick unsaturated zone. Before major aquifer exploitation began about 50 years ago, 40–50% of the total discharge into the Valley was by upward flow through the lacustrine deposits. The best results were obtained using a subsurface distribution of hydrostratigraphic units based on recently published geological interpretations.     

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.     

Tile drainage causes flashy streamflow response in Ohio watersheds

Artificial subsurface (tile) drainage is used to increase trafficability and crop yield in much of the Midwest due to soils with naturally poor drainage. Tile drainage has been researched extensively at the field scale, but knowledge gaps remain on how tile drainage influences the streamflow response at the watershed scale. The purpose of this study is to analyse the effect of tile drainage on the streamflow response for 59 Ohio watersheds with varying percentages of tile drainage and explore patterns between the Western Lake Erie Bloom Severity Index to streamflow response in heavily tile-drained watersheds. Daily streamflow was downloaded from 2010 to 2019 and used to calculated mean annual peak daily runoff, mean annual runoff ratio, the percent of observations in which daily runoff exceeded mean annual runoff (T_Qmean), baseflow versus stormflow percentages, and the streamflow recession constant. Heavily-drained watersheds (>40% of watershed area) consistently reported flashier streamflow behaviour compared to watersheds with low percentages of tile drainage (<15% of watershed area) as indicated by significantly lower baseflow percentages, T_Qmean, and streamflow recession constants. The mean baseflow percent for watersheds with high percentages of tile drainage was 20.9% compared to 40.3% for watersheds with low percentages of tile drainage. These results are in contrast to similar research regionally indicating greater baseflow proportions and less flashy hydrographs (higher T_Qmean) for heavily-drained watersheds. Stormflow runoff metrics in heavily-drained watersheds were significantly positively correlated to western Lake Erie algal bloom severity. Given the recent trend in more frequent large rain events and warmer temperatures in the Midwest, increased harmful algal bloom severity will continue to be an ecological and economic problem for the region if management efforts are not addressed at the source. Management practices that reduce the streamflow response time to storm events, such as buffer strips, wetland restoration, or drainage water management, are likely to improve the aquatic health conditions of downstream communities by limiting the transport of nutrients following storm events.     

Homogenization of spatial patterns of hydrologic response in artificially drained agricultural catchments

Anthropogenic modifications to the landscape, with agricultural activities being a primary driver, have resulted in significant alterations to the hydrologic cycle. Artificial drainage, including surface and subsurface drainage (tile drains), is one of the most extensive manipulations in agricultural landscapes and thus is expected to provide a distinct signature of anthropogenic modification. This study adopts a data synthesis approach in an effort to characterize the signature of artificial subsurface drainage. Daily discharge data from 24 basins across the state of Iowa, which encapsulate a range of anthropogenic modifications, are assessed using a variety of flow metrics. Results indicate that the presence of artificial subsurface drainage leads to a homogenization of landscape hydrologic response. Non‐tiled watersheds exhibit a decrease in the area‐normalized peak discharge and an increase in the baseflow ratio (baseflow/streamflow) with increases in the spatial scale, while scale invariance is apparent in tiled basins. Within‐basin variability in hydrograph recession coefficients also appears to decrease with increases in the proportion of the catchment that is artificially drained. Finally, the differences between tiled and non‐tiled landscapes disappear at scales greater than approximately 2200 km², indicating that this may be a threshold scale for studying the effects of tile drainage. This decrease in within‐basin variability and the scale invariance of hydrologic metrics in artificially drained watersheds are attributed to the creation of a bypass flow hydrologic pathway that bypasses the complexity of the catchment travel paths. Spatial homogeneity in responses implies that it may be possible to develop more parsimonious hydrologic models for these regions. Copyright © 2013 John Wiley & Sons, Ltd.     

Thermal Imagery of Groundwater Seeps: Possibilities and Limitations

Quantifying groundwater flow at seepage faces is crucial because seepage faces influence the hydroecology and water budgets of watersheds, lakes, rivers and oceans, and because measuring groundwater fluxes directly in aquifers is extremely difficult. Seepage faces provide a direct and measurable groundwater flux but there is no existing method to quantitatively image groundwater processes at this boundary. Our objective is to determine the possibilities and limitations of thermal imagery in quantifying groundwater discharge from discrete seeps. We developed a conceptual model of temperature below discrete seeps, observed 20 seeps spectacularly exposed in three dimensions at an unused limestone quarry and conducted field experiments to examine the role of diurnal changes and rock face heterogeneity on thermal imagery. The conceptual model suggests that convective air‐water heat exchange driven by temperature differences is the dominant heat transfer mechanism. Thermal imagery is effective at locating and characterizing the flux of groundwater seeps. Areas of active groundwater flow and ice growth can be identified from thermal images in the winter, and seepage rates can be differentiated in the summer. However, the application of thermal imagery is limited by diverse factors including technical issues of image acquisition, diurnal changes in radiation and temperature, and rock face heterogeneity. Groundwater discharge rates could not be directly quantified from thermal imagery using our observations but our conceptual model and experiments suggest that thermal imagery could quantify groundwater discharge when there are large temperature differences, simple cliff faces, non‐freezing conditions, and no solar radiation.     

Deduction of groundwater flow regime in a basaltic aquifer using geochemical and isotopic data: The Golan Heights, Israel case study

Groundwater flow-paths through shallow-perch and deep-regional basaltic aquifers at the Golan Heights, Israel, are reconstructed by using groundwater chemical and isotopic compositions. Groundwater chemical composition, which changes gradually along flow-paths due to mineral dissolution and water–rock interaction, is used to distinguish between shallow-perched and deep-regional aquifers. Groundwater replenishment areas of several springs are identified based on the regional depletion in rainwater δ¹⁸O values as a function of elevation (−0.25‰ per 100 m). Tritium concentrations assist in distinguishing between pre-bomb and post-bomb recharged rainwater.

It was found that waters emerging through the larger springs are lower in δ¹⁸O than surrounding meteoric water and poor in tritium; thus, they are inferred to originate in high-elevation regions up to 20 km away from their discharge points and at least several decades ago. These results verify the numerically simulated groundwater flow field proposed in a previous study, which considered the geological configuration, water mass balance and hydraulic head spatial distribution.     

Increasing non-linearity of the storage-discharge relationship in sub-Arctic catchments

The Arctic is warming at an unprecedented rate. We hypothesis that as seasonally frozen soils thaw and recede in extent as a response to this warming, flow path diversity and thus hydrologic connectivity increases. This enhanced hydrologic connectivity then increases the non-linearity of the storage-discharge relationship in a catchment. The objective of this study is to test this hypothesis by quantifying trends and spatio-temporal differences in the degree of linearity in the storage-discharge relationships for 16 catchments within Northern Sweden from 1950 to 2018. We demonstrate a clear increase in non-linearity of the storage-discharge relationship over time for all catchments with 75% showing a statistically significant increase in non-linearity. Spring has significantly more linear storage-discharge relationships than summer for most catchments (75%) supporting the idea that seasonally frozen soils with a low degree of hydrological connectivity have a linear storage-discharge relationship. For the period considered, spring also showed greater change in storage-discharge relationship trends than summer signifying that changes in recessions are primarily occurring during the thawing period. Separate storage-discharge analyses combined with preceding winter conditions demonstrated that especially cold winters with little snow yielded springs and summers with more linear storage-discharge relationships. We show that streamflow recession analysis reflects ongoing hydrological change of an arctic landscape as well as offering new metrics for tracking change across arctic and sub-arctic landscapes.     

How to construct recursive digital filters for baseflow separation

Recursive digital filtering of hydrographs is a baseflow separation method that can easily be automated and has been recommended for providing reproducible results. In the past, different formulations of the most simple filter type, the so‐called one‐parameter filter, have been proposed. In this paper, a theoretical framework is developed for filter algorithms that were constructed under the assumption that the outflow from an aquifer is linearly proportional to its storage. It is shown that these one‐parameter filters describing an exponential baseflow recession are all special cases of a two‐parameter filter whose equation is specified. Its parameters are the recession constant—which can be objectively determined by a recession analysis—and BFI_max, the maximum value of the baseflow index that can be modelled by the algorithm. This introduces a subjective element into the baseflow calculation, since BFI_max is not measurable. A preliminary analysis based on the results of conventional separation techniques shows that it might be possible to find typical BFI_max values for classes of catchments that can be unequivocally distinguished by their hydrological and hydrogeological characteristics. Copyright © 2004 John Wiley & Sons, Ltd.     

Run‐off processes from mountains to foothills: The role of soil stratigraphy and structure in influencing run‐off characteristics across high to low relief landscapes

The critical zone features that control run‐off generation, specifically at the regional watershed scale, are not well understood. Here, we addressed this knowledge gap by quantitatively and conceptually linking regional watershed‐scale run‐off regimes with critical zone structure and climate gradients across two physiographic provinces in the Southeastern United States. We characterized long‐term (~20 years) discharge and precipitation regimes for 73 watersheds with United States Geological Survey in‐stream gaging stations across the Appalachian Mountain and Piedmont physiographic provinces of North Carolina. Watersheds included in this analysis had <10% developed land and ranged in size from 14.1–4,390 km². Thirty‐four watersheds were located in the Piedmont physiographic province, which is typically classified as a low relief landscape with deep, highly weathered soils and regolith. Thirty‐nine watersheds were located in the Appalachian Mountain physiographic province, which is typically classified as a steeper landscape with highly weathered, but shallower soils and regolith. From the United States Geological Survey daily mean run‐off time series, we calculated annual and seasonal baseflow indices (BFI), minimum, mean, and maximum daily run‐off, and Pearson's correlation coefficients between precipitation and baseflow. Our results showed that Appalachian Mountain watersheds systematically had higher minimum daily flows and BFI values. Piedmont watersheds displayed much larger deviations from mean annual BFI in response to year‐to‐year variability in precipitation. A series of linear regression models between 21 landscape metrics and annual BFIs showed non‐linear and complex terrestrial–hydrological relationships across the two provinces. From these results, we discuss how distinct features of critical zone architecture, with specific focus on soil depth and stratigraphy, may be dominating the regulation of hydrological processes and run‐off regimes across these provinces.