Changes in streamflow components following logging and regeneration in the southern forest of Western Australia

Two small experimental catchments were established in the south-west of Western Australia to study the effects of logging and subsequent regeneration on the mechanism of streamflow generation. Following a six year pre-treatment calibration period (1976–1981), one catchment (March Road) was logged and reforested in 1982 and the other (April Road South) remained as a control. Logging resulted in an increase in groundwater levels and subsequently groundwater discharge area. The deep, permanent groundwater levels in the valley and upslope areas rose until 1986 and then began to decline. The maximum rise was 5 m in the upslope areas. The duration of shallow, intermittent groundwater system, perched on underlying clay, was extended from 2–3 months in winter before logging to 5–6 months after logging. The shallow groundwater level rose in the valley and began to discharge at the ground surface in 1986. Logging resulted in an increase in streamflow. The maximum increase (≈18% of annual rainfall) was in 1983, one year after logging. The increase in streamflow was due to a substantial decrease in interception and evapotranspiration, increased recharge to the shallow groundwater system, decreased soil moisture deficit and consequently an increase in throughflow. The increase in base flow was about twice that of quick flow. The changes in streamflow and its components in the subsequent years were closely related to the groundwater discharge area. Most of the quick flow was generated as saturation excess overland flow from the groundwater discharge area in the valley. The expansion of the groundwater discharge area, increased soil moisture content, higher groundwater level and the presence of the shallow groundwater system for the extended periods were responsible for the process of streamflow generation.     

The role of run-on for overland flow and the characteristics of runoff generation in the Loess Plateau,China

Abstract

The Loess Plateau in China is overlain by deep and loose soil. As in other semi-arid regions, convective precipitation produces storms, typically of short duration, relatively high intensity and limited areal extent. Infiltration excess (Hortonian mechanism) of precipitation is conventionally assumed to be more prominent than saturation excess (Dunne mechanism) for storm runoff generation. This assumption is true at a point during the storm. However, the runoff generation mechanism is altered when the runoff is conditioned by a lateral redistribution movement of water, i.e. run-on, as the spatial scale increases. In the Loess Plateau, the effects of run-on may be significant, because of the deep and loose surface soil layer. In this study, the role of run-on for overland flow in the Upper Wei River basin, located in the Loess Plateau, is evaluated by means of a simple numerical model at the hillslope scale. The results show that almost all the Hortonian overland flow infiltrates into the soil along the flat hillslope and dry gully before it reaches the river channel. Most of the runoff is generated from the saturated soil near the river channel and from the subsurface. The run-on process takes much longer than the infiltration, facilitating rainfall–runoff modelling at a daily time step. A hydrological model is employed to investigate the characteristics of runoff generation in the Upper Wei River basin. The analysis shows that the subsurface flow contribution to total streamflow is more than 53% from October to March, while the overland flow contribution exceeds 72% from April to September.

Editor D. Koutsoyiannis; Associate editor Dawen Yang

Citation Liu, D.F., Tian, F.Q., Hu, H.C., and Hu, H.P., 2012. The role of run-on for overland flow and the characteristics of runoff generation in the Loess Plateau, China. Hydrological Sciences Journal, 57 (6), 1107–1117.     

The role of topography in controlling throughflow generation

Previous studies of the downslope movement of soil moisture have been largely confined to two-dimensional slope profiles—the influence of contour curvature (three-dimensional) on moisture movement has received little attention. An automatic tensiometer system has been used to monitor soil moisture status continuously in a single hollow and adjacent spurs. Maximum saturated hillslope flow is shown to coincide with the stream discharge ‘throughflow peak’. Relative to the hollow, little saturation is evident at the base of the spurs. This subcatchment model is applied to a major portion of the catchment, showing the general importance of topography in the control of stream discharge.     

Snowmelt runoff and soil moisture dynamics on steep subalpine hillslopes

Snowmelt water supplies streamflow and growing season soil moisture in mountain regions, yet pathways of snowmelt water and their effects on moisture patterns are still largely unknown. This study examined how flow processes during snowmelt runoff affected spatial patterns of soil moisture on two steep sub‐alpine hillslope transects in Rocky Mountain National Park, CO, USA. The transects have northeast‐facing and east‐facing aspects, and both extend from high‐elevation bedrock outcrops down to streams in valley bottoms. Spatial patterns of both snow depth and near‐surface soil moisture were surveyed along these transects in the snowmelt and summer seasons of 2008–2010. To link these patterns to flow processes, soil moisture was measured continuously on both transects and compared with the timing of discharge in nearby streams. Results indicate that both slopes generated shallow lateral subsurface flow during snowmelt through near‐surface soil, colluvium and bedrock fractures. On the northeast‐facing transect, this shallow subsurface flow emerged through mid‐slope seepage zones, in some cases producing saturation overland flow, whereas the east‐facing slope had no seepage zones or overland flow. At the hillslope scale, earlier snowmelt timing on the east‐facing slope led to drier average soil moisture conditions than on the northeast‐facing slope, but within hillslopes, snow patterns had little relation to soil moisture patterns except in areas with persistent snow drifts. Results suggest that lateral flow and exfiltration processes are key controls on soil moisture spatial patterns in this steep sub‐alpine location. Copyright © 2014 John Wiley & Sons, Ltd.     

Interactions between perched and saprolite aquifers on a small,salt-affected and deeply weathered hillslope

A small hillslope was chosen to investigate the role of throughflow as a mechanism responsible for the movement of soil water and solutes towards a saline seep and as a source of recharge to a permanent, regional aquifer at depth. The hydraulic properties, chemical characteristics and physical responses of both systems were studied on a deeply weathered, salt-affected hillslope. Additional data were also obtained from other sites in south-western Australia. Regional groundwater flow occurred in a variably textured, deeply weathered material in which the hydraulic conductivity varied from < 0·001 to 0·14m day^?1. Perched groundwater flow (throughflow) occurred in the higher permeability (? 1·5 m day^?1), near-surface soil materials. Throughflow occurred throughout winter, contributing approximately 530 m³ of fresh (? 160 mg l^?1 Cl) water to a saline seep. By contrast, the deep aquifer discharged approximately 1100 m³ of waters with salt concentrations of 2000–6000 mg l^?1 Cl. Recharge and discharge rates to and from the deep aquifer, were estimated to be of the order of 5–20 mm a^?1 and 50–300 mm a^?1 respectively. Saturated conditions existed throughout winter within the seep and the immediately adjacent non-saline area, with up to 60 per cent of the hillslope soils becoming saturated after major rainfall events ( > 20 mm day^?1). In the mid-slopes, in particular along a central depression, saturation of the shallow soils caused macropore channel recharge to take place through the clay-textured subsoils. Water-level responses suggest that approximately 25–30 per cent of annual recharge occurred from one storm studied in September 1984. Recharge through macropore channels is a significant mechanism in the concave slope areas on the hillslope. Throughflow was found to be a major source of water, but not salt, contributing to the saline seep. In general, the contribution of throughflow was found to decrease further inland at other sites studied. However, at inland sites where perennial, perched aquifers have developed in deep sands, saline areas have been caused by throughflow and not by deep aquifer discharge.     

Stormflow generation and flowpath characteristics in an Amazonian rainforest catchment

The Amazon basin covers an area of roughly 7 × 10⁶ km² and encompasses diverse soil – landscape types with potentially differing hydrological behaviour. This study was conducted in the Ultisol landscape of the western Amazon basin in Peru. Processes of stormflow generation were investigated on an event basis in a first‐order rainforest catchment to establish a causal link between soil physical and precipitation characteristics, hillslope flowpaths and stormflow hydrograph attributes. A sharp decrease in soil hydraulic conductivity with depth and high rainfall intensity and frequency favour rapid near‐surface flowpaths, mainly in the form of saturation‐excess overland flow and return flow. The latter results in an almost random occurrence of overland flow, with no obvious topographic control. Hillslope flowpaths do not vary much with respect to the hydrograph attributes time of rise, response time, lag time and centroid lag time. They have the same response time as streamflow, but a somewhat lower time of rise and significantly shorter lag times. The recession constant for hillslope hydrographs is about 10 min, in contrast to the streamflow recession constants of 28, 75 and 149 min. Stormflow generation in this Ultisol rainforest catchment differs strongly from that reported for Oxisol rainforest catchments. These two soilscapes may define a spectrum of possible catchment hydrological behaviour in the Amazon basin. Copyright © 2000 John Wiley & Sons, Ltd.     

Preliminary evaluation of the runoff processes in a remote montane cloud forest basin using Mixing Model Analysis and Mean Transit Time

In this study, the Mean Transit Time and Mixing Model Analysis methods are combined to unravel the runoff generation process of the San Francisco River basin (73.5 km²) situated on the Amazonian side of the Cordillera Real in the southernmost Andes of Ecuador. The montane basin is covered with cloud forest, sub‐páramo, pasture and ferns. Nested sampling was applied for the collection of streamwater samples and discharge measurements in the main tributaries and outlet of the basin, and for the collection of soil and rock water samples. Weekly to biweekly water grab samples were taken at all stations in the period April 2007–November 2008. Hydrometric data, Mean Transit Time and Mixing Model Analysis allowed preliminary evaluation of the processes controlling the runoff in the San Francisco River basin. Results suggest that flow during dry conditions mainly consists of lateral flow through the C‐horizon and cracks in the top weathered bedrock layer, and that all subcatchments have an important contribution of this deep water to runoff, no matter whether pristine or deforested. During normal to low precipitation intensities, when antecedent soil moisture conditions favour water infiltration, vertical flow paths to deeper soil horizons with subsequent lateral subsurface flow contribute most to streamflow. Under wet conditions in forested catchments, streamflow is controlled by near surface lateral flow through the organic horizon. Exceptionally, saturation excess overland flow occurs. By absence of the litter layer in pasture, streamflow under wet conditions originates from the A horizon, and overland flow. Copyright © 2011 John Wiley & Sons, Ltd.     

PIPEFLOW CONTRIBUTING AREAS AND RUNOFF RESPONSE

Links that have been established between the characteristic responses of hillslope drainage processes and the size of drainage catchments enable valuable extrapolations and comparisons to be made, both within and beyond the range of the field data. To date, flow through natural soil pipes has not been included in these investigations. A methodology is offered for achieving this and the first such analyses of field data on pipeflow. The results indicate that pipeflow responses follow trends which are broadly similar to other hillslope drainage processes and are intermediate between throughflow and saturation overland flow. © 1997 by John Wiley & Sons, Ltd.     

Interactions and connectivity between runoff generation processes of different spatial scales

Monitoring runoff generation processes in the field is a prerequisite for developing conceptual hydrological models and theories. At the same time, our perception of hydrological processes strongly depends on the spatial and temporal scale of observation. Therefore, the aim of this study is to investigate interactions between runoff generation processes of different spatial scales (plot scale, hillslope scale, and headwater scale). Different runoff generation processes of three hillslopes with similar topography, geology and soil properties, but differences in vegetation cover (grassland, coniferous forest, and mixed forest) within a small v‐shaped headwater were measured: water table dynamics in wells with high spatial and temporal resolution, subsurface flow (SSF) of three 10 m wide trenches at the bottom of the hillslopes subdivided into two trench sections each, overland flow at the plot scale, and catchment runoff. Bachmair et al. ( 2012 ) found a high spatial variability of water table dynamics at the plot scale. In this study, we investigate the representativity of SSF observations at the plot scale versus the hillslope scale and vice versa, and the linkage between hillslope dynamics (SSF and overland flow) and streamflow. Distinct differences in total SSF within each 10 m wide trench confirm the high spatial variability of the water table dynamics. The representativity of plot scale observations for hillslope scale SSF strongly depends on whether or not wells capture spatially variable flowpaths. At the grassland hillslope, subsurface flowpaths are not captured by our relatively densely spaced wells (3 m), despite a similar trench flow response to the coniferous forest hillslope. Regarding the linkage between hillslope dynamics and catchment runoff, we found an intermediate to high correlation between streamflow and hillslope hydrological dynamics (trench flow and overland flow), which highlights the importance of hillslope processes in this small watershed. Although the total contribution of SSF to total event catchment runoff is rather small, the contribution during peak flow is moderate to substantial. Additionally, there is process synchronicity between spatially discontiguous measurement points across scales, potentially indicating subsurface flowpath connectivity. Our findings stress the need for (i) a combination of observations at different spatial scales, and (ii) a consideration of the high spatial variability of SSF at the plot and hillslope scale when designing monitoring networks and assessing hydrological connectivity. Copyright © 2013 John Wiley & Sons, Ltd.     

Parametric study of a physically-based,plot-scale hillslope hydrological model through virtual experiments

Abstract

A physically-based hillslope hydrological model with shallow overland flow and rapid subsurface stormflow components was developed and calibrated using field experiments conducted on a preferential path nested hillslope in northeast India. Virtual experiments were carried out to perform sensitivity analysis of the model using the automated parameter estimation (PEST) algorithm. Different physical parameters of the model were varied to study the resulting effects on overland flow and subsurface stormflow responses from the theoretical hillslopes. It was observed that topographical shapes had significant effects on overland flow hydrographs. The slope profiles, surface storage, relief, rainfall intensity and infiltration rates primarily controlled the overland flow response of the hillslopes. Prompt subsurface stormflow responses were mainly dominated by lateral preferential flow, as soil matrix flow rates were very slow. Rainfall intensity and soil macropore structures were the most influential parameters on subsurface stormflow. The number of connected soil macropores was a more sensitive parameter than the size of macropores. In hillslopes with highly active vertical and lateral preferential pathways, saturation excess overland flow was not evident. However, saturation excess overland flow was generated if the lateral macropores were disconnected. Under such conditions, rainfall intensity, duration and preferential flow rate governed the process of saturation excess overland flow generation from hillslopes.

Editor D. Koutsoyiannis; Associate editor C. Perrin     

Understanding plot‐scale hydrology of Lesser Himalayan watershed—A field study and HYDRUS‐2D modelling approach

Soil moisture dynamics have a significant effect on overland flow generation. Catchment aspect is one of the major controlling factors of overland flow and soil moisture behaviour. A few experimental studies have been carried out in the uneven topography of the Himalayas. This study presents plot‐scale experiments using portable rainfall simulator at an altitude of 1,230 m above mean sea level and modelling of overland flow using observed datasets. Two plots were selected in 2 different aspects of Aglar watershed of Lesser Himalaya; the agro‐forested (AF) plot was positioned at the north aspect whereas the degraded (DE) plot was located at the south aspect of the hillslope. HS flumes and rain gauges were installed to measure the runoff at the outlet of the plot and the rainfall depth during rainfall simulation experiments. Moreover, 10 soil moisture sensors were installed at upslope and downslope locations of both the plots at 5, 15, 25, 35, and 45 cm depth from ground level to capture the soil moisture dynamics. The tests were conducted at intensities of 79.8 and 75 mm/hr in AF plot and 82.2 and 72 mm/hr in the DE plot during Test 1 and Test 2, respectively. The observed data indicate the presence of reinfiltration process only in the AF plot. The high water holding capacity and the presence of reinfiltration process results in less runoff volume in the AF plot compared with the DE plot. The Hortonian overland flow mechanism was found to be the dominant overland flow mechanism as only a few layers of top soil get saturated during all of the rainfall–runoff experiments. The runoff, rainfall, and soil moisture data were subsequently used to calibrate the parameters of HYDRUS‐2D overland flow module to simulate the runoff hydrograph and soil moisture. The components of hydrograph were evaluated in terms of peak discharge, runoff volume and time of concentration, the results were found to be within the satisfactory range. The goodness of fit of simulated hydrographs were more than 0.85 and 0.95 for AF and DE plot, respectively. The model produced satisfactory simulation results of soil moisture for all of the rainfall–runoff experiments. The HYDRUS‐2D overland flow module was found promising to simulate the runoff hydrograph and soil moisture in plot‐scale research.     

Towards estimation of extreme floods: examination of the roles of runoff process changes and floodplain flows

This paper presents the development and application of a distributed rainfall-runoff model for extreme flood estimation, and its use to investigate potential changes in runoff processes, including changes to the ‘rating curve’ due to effects of over-bank flows, during the transition from ‘normal’ floods to ‘extreme’ floods. The model has two components: a hillslope runoff generation model based on a configuration of soil moisture stores in parallel and series, and a distributed flood routing model based on non-linear storage-discharge relationships for individual river reaches that includes the effects of floodplain geometries and roughnesses. The hillslope water balance model contains a number of parameters, which are measured or derived a priori from climate, soil and vegetation data or streamflow recession analyses. For reliable estimation of extreme discharges that may extend beyond recorded data, the parameters of the flood routing model are estimated from hydraulic properties, topographic data and vegetation cover of compound channels (main channel and floodplains). This includes the effects of the interactions between the main channel and floodplain sections, which tend to cause a change to the rating curve. The model is applied to the Collie River Basin, 2545 km², in Western Australia and used to estimate the probable maximum flood (PMF) from probable maximum precipitation estimates for this region. When moving from normal floods to the PMFs, application of the model demonstrates that the runoff generation process changes with a substantial increase of saturation excess overland flow through the expansion of saturated areas, and the dominant runoff process in the stream channel changes from in-bank to over-bank flows. The effects of floodplain inundation and floodplain vegetation can significantly reduce the magnitude of the estimated PMFs. This study has highlighted the need for the estimation of a number of critical parameters (e.g. cross-sectional geometry, floodplain vegetation, soil depths) through concerted field measurements or surveys, and targeted laboratory experiments.     

Sub‐grid scale parameterization of hillslope runoff and erosion processes for catchment‐scale models of semi‐arid landscapes

The processes of hillslope runoff and erosion are typically represented at coarse spatial resolution in catchment‐scale models due to computational limitations. Such representation typically fails to incorporate the important effects of topographic heterogeneity on runoff generation, overland flow, and soil erosion. These limitations currently undermine the application of distributed catchment models to understand the importance of thresholds and connectivity on hillslope and catchment‐scale runoff and erosion, particularly in semi‐arid environments. This paper presents a method for incorporating high‐resolution topographic data to improve sub‐grid scale parameterization of hillslope overland flow and erosion models. Results derived from simulations conducted using a kinematic wave overland flow model at 0.5 m spatial resolution are used to parameterize the depth–discharge relationship in the overland flow model when applied at 16 m resolution. The high‐resolution simulations are also used to derive a more realistic parameterization of excess flow shear stress for use in the 16 m resolution erosion model. Incorporating the sub‐grid scale parameterization in the coarse‐resolution model (16 m) leads to improved predictions of overland flow and erosion when evaluated using results derived from high‐resolution (0.5 m) model simulations. The improvement in performance is observed for a range of event magnitudes and is most notable for erosion estimates due to the non‐linear dependency between the rates of erosion and overland flow. Copyright © 2013 John Wiley & Sons, Ltd.     

Revegetation impacts on runoff generation processes in a typical loess catchment

Since 1999, large-scale ecosystem restoration has been implemented in the Loess Plateau, effectively increasing regional vegetation coverage. Vegetation restoration has significantly elevated the saturated hydraulic conductivity (Ks) of the near-surface soil layers and increased the vertical heterogeneity of the Ks profile. Many studies have examined the change of runoff due to revegetation, yet the impacts of Ks profile on the soil moisture distribution and runoff generation processes were less explored. In this study, numerical simulations were conducted to investigate how changes in the Ks profile caused by vegetation restoration influenced the hydrological responses at event scale. The numerical simulation results show that the increase of surface Ks caused by vegetation restoration can effectively reduce runoff at event scale. Moreover, the enhancement of vertical heterogeneity of Ks profiles can significantly change the vertical profile of soil water content, prompting more water to percolate into the deep soil layer. When rainfall exceeds a threshold, the accumulation of soil water above the relatively less permeable layer can cause short-term saturation in shallow soil layers, resulting in a transient perched water table. As a result, after the vegetation restoration in the Loess Plateau, though Horton overland flow is still the main runoff generation mechanism, there is a possibility of the emergence of Dunne overland flow under the high vegetation coverage (e.g., NDVI larger than 0.5). This emergence of new runoff generation mechanism, saturation excess runoff, in the Loess Plateau due to the vegetation restoration could provide scientific guidance for water and sediment movement, soil and water conservation practices, and desertification control in the Loess Plateau.     

A numerical investigation of bedrock groundwater recharge and exfiltration on soil mantled hillslopes

Here we use Richards Equation models of variably saturated soil and bedrock groundwater flow to investigate first-order patterns of the coupling between soil and bedrock flow systems. We utilize a Monte Carlo sensitivity analysis to identify important hillslope parameters controlling bedrock recharge and then model the transient response of bedrock and soil flow to seasonal precipitation. Our results suggest that hillslopes can be divided into three conceptual zones of groundwater interaction, (a) the zone of lateral unsaturated soil moisture accumulation (upper portion of hillslope), (b) the zone of soil saturation and bedrock recharge (middle of hillslope) and (c) the zone of saturated-soil lateral flow and bedrock groundwater exfiltration (bottom of hillslope). Zones of groundwater interaction expand upslope during periods of precipitation and drain downslope during dry periods. The amount of water partitioned to the bedrock groundwater system a can be predicted by the ratio of bedrock to soil saturated hydraulic conductivity across a variety of hillslope configurations. Our modelled processes are qualitatively consistent with observations of shallow subsurface saturation and groundwater fluctuation on hillslopes studied in our two experimental watersheds and support a conceptual model of tightly coupled shallow and deep subsurface circulation where groundwater recharge and discharge continuously stores and releases water from longer residence time storage.     

On the contribution of groundwater to streamflow in laterite catchments of the Darling Range,south-western Australia

In deeply weathered laterite catchments of the Darling Range in south-western Australia, the direct contribution (i.e., discharge) of permanent groundwater to streamflow has long been considered as minor. Instead, downslope shallow throughflow was thought to dominate, generating more than 90% of streamflow. We used a chemical hydrograph separation approach to estimate annual groundwater discharge for three catchments over periods of up to 39 years, and found that direct groundwater contributions to streamflow were far more variable across catchments and through time than has previously been acknowledged. The estimated proportion of annual streamflow sourced directly from groundwater ranged from 0 to 93% and was related linearly to the size of the groundwater discharge area in the catchment valley floor. In contrast, contributions from shallow sources including shallow throughflow varied primarily and linearly with annual rainfall. However, the response to rainfall was “amplified” in a predictable way by the size of the groundwater discharge area, consistent with the variable source area concept. We derived a functional relationship between catchment annual rainfall-runoff ratio and groundwater discharge area and successfully applied this to a further four catchments, inferring that the results were broadly applicable across the Darling Range. The implications for an improved understanding of streamflow generating processes in the study region, and for laterite catchments generally, are discussed.     

Seasonal changes in runoff generation in a small forested mountain catchment

This study aimed to investigate the seasonal variability of runoff generation processes, the sources of stream water, and the controls on the contribution of event water to streamflow for a small forested catchment in the Italian pre‐Alps. Hydrometric, isotopic, and electrical conductivity data collected between August 2012 and August 2013 revealed a marked seasonal variability in runoff responses. Noticeable differences in runoff coefficients and hydrological dynamics between summer and fall/spring rainfall events were related to antecedent moisture conditions and event size. Two‐component and three‐component hydrograph separation and end‐member mixing analysis showed an increase in event water contributions to streamflow with event size and average rainfall intensity. Event water fractions were larger during dry conditions in the summer, suggesting that stormflow generation in the summer consisted predominantly of direct channel precipitation and some saturated overland flow from the riparian zone. On the contrary, groundwater and hillslope soil water contributions dominated the streamflow response during wet conditions in fall. Seasonal differences were also noted between event water fractions computed based on isotopic and electrical conductivity data, likely because of the dilution effect during the wetter months. Copyright © 2014 John Wiley & Sons, Ltd.     

Water flux and sediment transport within a forested landscape: the role of connectivity,subsurface flow,and slope length scale on transport mechanism

The connectivity and upscaling of overland runoff and sediment transport are important issues in hillslope hydrology to identify water flux and sediment transport within landscape. These processes are highly variable in time and space with regard to their interactions with vegetation and soil surface conditions. The generation of overland runoff and its spatial connectivity were examined along a slope to determine the variations in the transport mechanism of runoff and soil particles by rain splash and overland runoff. Field experiments were conducted by erosion plots on a steep hillslope at lengths of 5, 10, and 15 m. The overland runoff connectivity and flow transport distance decreased with the slope length, while spatial variability of infiltration increased significantly with the slope length. Observation of subsurface flow revealed that surface soil and litter layer could have important role in water transport. However, the surface soil water content and water flux transport along the slope was highly variable for different storm events; the variability was related to the complexity of the system, mainly by way of the initial wetness conditions and infiltration characteristics. Only net rain‐splashed soil was measurable, but examination of the water flux, overland runoff and sediment transport connectivity, characteristics of sheetwash, and the variability in spatial infiltration indicated an increase in the contribution of the rain splash transport mechanism along the slope. Copyright © 2013 John Wiley & Sons, Ltd.