How does afforestation affect the hydrology of a blanket peatland? A modelling study

Over the last century, afforestation in Ireland has increased from 1% of the land area to 10%, with most plantations on upland drained blanket peatlands. This land use change is considered to have altered the hydrological response and water balance of upland catchments with implications for water resources. Because of the difficulty of observing these long‐term changes in the field, the aim of this study was to utilize a hydrological model to simulate the rainfall runoff processes of an existing pristine blanket peatland and then to simulate the hydrology of the peatland if it were drained and afforested. The hydrological rainfall runoff model (GEOtop) was calibrated and validated for an existing small (76 ha) pristine blanket peatland in the southwest of Ireland for the 2‐year period, 2007–2008. The current hydrological response of the pristine blanket peatland catchment with regard to streamflow and water table (WT) levels was captured well in the simulations. Two land use change scenarios of afforestation were also examined, (A) a young 10‐year‐old and (B) a semi‐mature 15‐year‐old Sitka Spruce forest. Scenario A produced similar streamflow dynamics to the pristine peatland, whereas total annual streamflow from Scenario B was 20% lower. For Scenarios A and B, on an annual average basis, the WT was drawn down by 16 and 20 cm below that observed in the pristine peatland, respectively. The maximum WT draw down in Scenario B was 61 cm and occurred in the summer months, resulting in a significant decrease in summer streamflow. Occasionally in the winter (following rainfall), the WT for Scenario B was just 2 cm lower than the pristine peatland, which when coupled with the drainage networks associated with afforestation led to higher peak streamflows. Copyright © 2012 John Wiley & Sons, Ltd.     

Estimating external magnetic field differences at high geomagnetic latitudes from a single station

Providing an accurate estimate of the magnetic field on the Earth's surface at a location distant from an observatory has useful scientific and commercial applications, such as in repeat station data reduction, space weather nowcasting or aeromagnetic surveying. While the correlation of measurements between nearby magnetic observatories at low and mid‐latitudes is good, at high geomagnetic latitudes () the external field differences between observatories increase rapidly with distance, even during relatively low magnetic activity. Thus, it is of interest to describe how the differences (or errors) in external magnetic field extrapolation from a single observatory grow with distance from its location. These differences are modulated by local time, seasonal and solar cycle variations, as well as geomagnetic activity, giving a complex temporal and spatial relationship. A straightforward way to describe the differences are via confidence intervals for the extrapolated values with respect to distance. To compute the confidence intervals associated with extrapolation of the external field at varying distances from an observatory, we used 695 station‐years of overlapping minute‐mean data from 37 observatories and variometers at high latitudes from which we removed the main and crustal fields to isolate unmodelled signals. From this data set, the pairwise differences were analysed to quantify the variation during a range of time epochs and separation distances. We estimate the 68.3%, 95.4% and 99.7% confidence levels (equivalent to the 1σ, 2σ and 3σ Gaussian error bounds) from these differences for all components. We find that there is always a small non‐zero bias that we ascribe to instrumentation and local crustal field induction effects. The computed confidence intervals are typically twice as large in the north–south direction compared to the east‐west direction and smaller during the solstice months compared to the equinoxes.     

Spatial variability of hydraulic conductivity and bulk density along a blanket peatland hillslope

This article presents the results of a field investigation of saturated hydraulic conductivity K_sat and bulk density (ρ_bd) in an Atlantic blanket bog in the southwest of Ireland. Starting at a peatland stream and moving along an uphill transect toward the peatland interior, ρ_bd and K_sat were examined at regular intervals. Saturated horizontal hydraulic conductivity (Kh_sat) and vertical (Kv_sat) was estimated at two depths: 10–20 and 30–40 cm below the peat surface, whereas ρ_bd was estimated for the full profile. We consider two separate zones, one a riparian zone extending 10 m from the stream and a second zone in the bog interior. We found that the K_sat was higher (~10^–5 m s^–1) in the bog interior than that in the riparian zone (~10^–6 m s^–1), whereas the converse applied to bulk density, with lowest density (~0.055 g cm^–3) at the interior and highest (~0.11 g cm^–3) at the riparian zone. In general, we found Kh_sat to be approximately twice the Kv_sat. These results support the idea that the lower K_sat at the margins control the hydrology of blanket peatlands. It is therefore important that the spatial variation of these two key properties be accommodated in hydrological models if the correct rainfall runoff characteristics are to be correctly modelled. Stream flow analysis over 3 years at the peatland catchment outlet showed that the stream runoff was composed of 8% base flow and 92% flood flow, suggesting that this blanket peatland is a source rather than a sink for floodwaters. Copyright © 2011 John Wiley & Sons, Ltd.     

Downstream hydraulic geometry of rivers in Victoria, Australia

Downstream variation of hydraulic geometry in rivers, characterized by fine textured banks and low width to depth ratios (7–25), is investigated in Victoria, Australia, with the aim of developing predictive models of channel geometry for large-scale spatial modeling applications. A one-dimensional hydraulic model is used to determine the mean bank-full geometry and discharge (Q_bf) for 93 sites which are investigated in relation to discharge of fixed average recurrence interval (ARI). The median bank-full ARI is estimated at 0.8 years with 75% of sites between 0.5 and 2.5 years. Exponents in the downstream hydraulic geometry relations for width, depth and velocity are respectively 0.43, 0.40 and 0.18 (Q = Q_bf) and 0.44, 0.38 and 0.03 (Q = Q₂, i.e., 2-year ARI), falling near the mode of global values. Q₂ and slope explain 66% of variance in Q_bf, while Q₂ explains 73% and 69% of the variance in width and depth relations, respectively: Q₂ provides a reliable substitute for Q_bf in spatial modeling applications. Spatial variation in hydraulic geometry relations within and between river basins remains largely unexplained. The W/D ratio characteristically decreases with increasing distance along the lower reaches of most rivers and this has contributed to the lower than expected value for the width exponent.     

Sediment partitioning,continental slopes and base‐of‐slope systems

Deciphering the role slope topography plays in partitioning sediment on siliciclastic continental slope and base‐of‐slope systems helps our understanding of slope depositional processes in significant ways: (1) by validation of large‐scale depositional process models for continental margins, (2) by validation of numerical basin‐scale stratigraphic forward models used to test and deploy source‐to‐sink (S2S) concepts and (3) by creating models for setting reservoir presence and quality expectations in frontier areas poorly constrained by wells and seismic. A global database consisting of >700 km of drilled stratigraphy provide empirical rock data lacking from most S2S studies. Analysis of calibrated seismic stratigraphic units characterised using the contextual framework laid out in this paper show that both gross depositional environments (GDEs) and sand content occur across slope profiles in systematic ways. The challenge in using these observations to quantify reservoir risk and uncertainty lies with relating the observations to depositional processes that can be used to characterise frontier basins that lack calibration. Depositional process‐based understanding encoded in 3D stratigraphic forward models (SFM) can simulate both lithologies and GDEs providing broad predictions for exploration at the scale of an entire basin or slope system. Stratigraphic forward models allow the integration of S2S understanding and provide a framework for testing sediment‐partitioning hypotheses in frontier settings. Valid S2S models must balance sediment yield from the source catchments with sinks, and be consistent with basin specific observations. The proportions of GDEs across the slope provide additional validation criteria to ensure the models are plausible.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 2. Controls on weathering and permeability profiles,drainage hydraulics,and solute export pathways

The advance of a chemical weathering front into the bedrock of a hillslope is often limited by the rate weathering products that can be carried away, maintaining chemical disequilibrium. If the weathering front is within the saturated zone, groundwater flow downslope may affect the rate of transport and weathering—however, weathering also modifies the rock permeability and the subsurface potential gradient that drives lateral groundwater flow. This feedback may help explain why there tends to be neither “runaway weathering” to great depth nor exposed bedrock covering much of the earth and may provide a mechanism for weathering front advance to keep pace with incision of adjacent streams into bedrock. This is the second of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory. Here, we show how a simplified kinetic model of 1‐D rock weathering can be extended to consider lateral flow in a 2‐D hillslope. Exact and approximate analytical solutions for the location and thickness of weathering within the hillslope are obtained for a number of cases. A location for the weathering front can be found such that lateral flow is able to export weathering products at the rate required to keep pace with stream incision at steady state. Three pathways of solute export are identified: “diffusing up,” where solutes diffuse up and away from the weathering front into the laterally flowing aquifer; “draining down,” where solutes are advected primarily downward into the unweathered bedrock; and “draining along,” where solutes travel laterally within the weathering zone. For each pathway, a different subsurface topography and overall relief of unweathered bedrock within the hillslope is needed to remove solutes at steady state. The relief each pathway requires depends on the rate of stream incision raised to a different power, such that at a given incision rate, one pathway requires minimal relief and, therefore, likely determines the steady‐state hillslope profile.     

A low‐dimensional model of bedrock weathering and lateral flow coevolution in hillslopes: 1. Hydraulic theory of reactive transport

This is the first of a two‐part paper exploring the coevolution of bedrock weathering and lateral flow in hillslopes using a simple low‐dimensional model based on hydraulic groundwater theory (also known as Dupuit or Boussinesq theory). Here, we examine the effect of lateral flow on the downward fluxes of water and solutes through perched groundwater at steady state. We derive analytical expressions describing the decline in the downward flux rate with depth. Using these, we obtain analytical expressions for water age in a number of cases. The results show that when the permeability field is homogeneous, the spatial structure of water age depends qualitatively on a single dimensionless number, Hi. This number captures the relative contributions to the lateral hydraulic potential gradient of the relief of the lower‐most impermeable boundary (which may be below the weathering front within permeable or incipiently weathered bedrock) and the water table. A “scaled lateral symmetry” exists when Hi is low: age varies primarily in the vertical dimension, and variations in the horizontal dimension x almost disappear when the vertical dimension z is expressed as a fraction z/H(x) of the laterally flowing system thickness H(x). Taking advantage of this symmetry, we show how the lateral dimension of the advection–diffusion‐reaction equation can be collapsed, yielding a 1‐D vertical equation in which the advective flux downward declines with depth. The equation holds even when the permeability field is not homogeneous, as long as the variations in permeability have the same scaled lateral symmetry structure. This new 1‐D approximation is used in the accompanying paper to extend chemical weathering models derived for 1‐D columns to hillslope domains.     

The effects of disproportional load contributions on quantifying vegetated filter strip sediment trapping efficiencies

Vegetated filter strips (VFSs) are a best management practice (BMP) commonly implemented adjacent to row-cropped fields to trap overland transport of sediment and other constituents present in agricultural runoff. Although they have been widely adopted, insufficient data exist to understand their short and long-term effectiveness. High inter-event variability in performance has been observed, yet the majority of studies report average removal efficiencies over observed or simulated events, ignoring the disproportional effects of loads into and out of VFSs over longer periods of time. We argue that due to positively correlated sediment concentration-discharge relationships, disproportional contribution of runoff events transporting sediment over the course of a year (i.e., temporal inequality), decreased performance with increasing flow rates, and effects of antecedent moisture condition, VFS removal efficiencies over annual time scales may be significantly lower than reported per-event averages. By applying a stochastic approach, we investigated the extent of disparity between reporting average efficiencies from each runoff event over the course of 1 year versus the total annual load reduction. Additionally, we examined the effects of soil texture, concentration-discharge relationship, and VFS slope in contributing to this disparity, with the goal of revealing potential errors that may be incurred by ignoring the effects of temporal inequality in quantifying VFS performance. Simulation results suggest that ignoring temporal inequality can lead to overestimation of annual performance by as little as < 2% and to as much as > 20%, with the greatest disparities observed for soils with high clay content.