Predicting quantiles of water quality from catchment characteristics

Water quality is often highly variable both in space and time, which poses challenges for modelling the more extreme concentrations. This study developed an alternative approach to predicting water quality quantiles at individual locations. We focused on river water quality data that were collected over 25 years, at 102 catchments across the State of Victoria, Australia. We analysed and modelled spatial patterns of the 10th, 25th, 50th, 75th and 90th percentiles of the concentrations of sediments, nutrients and salt, with six common constituents: total suspended solids (TSS), total phosphorus (TP), filterable reactive phosphorus (FRP), total Kjeldahl nitrogen (TKN), nitrate-nitrite (NO_x), and electrical conductivity (EC). To predict the spatial variation of each quantile for each constituent, we developed statistical regression models and exhaustively searched through 50 catchment characteristics to identify the best set of predictors for that quantile. The models predict the spatial variation in individual quantiles of TSS, TKN and EC well (66%–96% spatial variation explained), while those for TP, FRP and NO_x have lower performance (37%–73% spatial variation explained). The most common factors that influence the spatial variations of the different constituents and quantiles are: annual temperature, percentage of cropping land area in catchment and channel slope. The statistical models developed can be used to predict how low- and high-concentration quantiles change with landscape characteristics, and thus provide a useful tool for catchment managers to inform planning and policy making with changing climate and land use conditions.     

Determination of subcatchment and watershed boundaries in a complex and highly urbanized landscape

Urban development significantly alters the landscape by introducing widespread impervious surfaces, which quickly convey surface run‐off to streams via stormwater sewer networks, resulting in “flashy” hydrological responses. Here, we present the inadequacies of using raster‐based digital elevation models and flow‐direction algorithms to delineate large and highly urbanized watersheds and propose an alternative approach that accounts for the influence of anthropogenically modified land cover. We use a semi‐automated approach that incorporates conventional drainage networks into overland flow paths and define the maximal run‐off contributing area. In this approach, stormwater pipes are clustered according to their slope attributes, which define flow direction. Land areas drained by each cluster and contributing (or exporting) flow to a topographically delineated catchment were determined. These land masses were subsequently added or removed from the catchment, modifying both the shape and the size. Our results in a highly urbanized Toronto, Canada, area watershed indicate a moderate net increase in the directly connected watershed area by 3% relative to a topographically forced method; however, differences across three smaller scale subcatchments are greater. Compared to topographic delineation, the directly connected watershed areas of both the upper and middle subcatchments decrease by 5% and 8%, respectively, whereas the lower subcatchment area increases by 15%. This is directly related to subsurface storm sewer pipes that cross topographic boundaries. When directly connected subcatchment area is plotted against total streamflow and flashiness indices using this method, the coefficients of variation are greater (0.93 to 0.97) compared to the use of digital elevation model‐derived subcatchment areas (0.78 to 0.85). The accurate identification of watershed and subcatchment boundaries should incorporate ancillary data such as stormwater sewer networks and retention basin drainage areas to reduce water budget errors in urban systems.     

Dynamic modelling of retrogressive landslides with emphasis on the role of clay sensitivity

This paper presents a detailed numerical study of the retrogressive failure of landslides in sensitive clays. The dynamic modelling of the landslides is carried out using a novel continuum approach, the particle finite element method, complemented with an elastoviscoplastic constitutive model. The multiwedge failure mode in the collapse is captured successfully, and the multiple retrogressive failures that have been widely observed in landslides in sensitive clays are reproduced with the failure mechanism, the kinematics, and the deposition being discussed in detail. Special attention has been paid to the role of the clay sensitivity on each retrogressive failure as well as on the final retrogression distance and the final run‐out distance via parametric studies. Moreover, the effects of the viscosity of sensitive clays on the failure are also investigated for different clay sensitivities.     

Wetlands and low-gradient topography are associated with longer hydrologic transit times in Precambrian Shield headwater catchments

The estimation of hydrologic transit times in a catchment provides insights into the integrated effects of water storage, mixing dynamics, and runoff generation processes. There has been limited effort to estimate transit times in southern boreal Precambrian Shield landscapes, which are characteristically heterogeneous with surface cover including till, thin soils, bedrock outcrops, and depressional wetland features that play contrasting hydrologic roles. This study presents approximately 3.5 years of precipitation and streamflow water isotope data and estimates mean transit times (MTTs) and the young water fraction (p_y) across six small catchments in the Muskoka-Haliburton region of south-central Ontario. The main objectives were to define a typical range of MTTs for headwater catchments in this region and to identify landscape variables that best explain differences in MTTs/p_y using airborne light detection and ranging and digital terrain analysis. Of the transit time distributions, the two parallel linear reservoir and gamma distributions best describe the hydrology of these catchments, particularly because of their ability to capture more extreme changes related to events such as snowmelt. The estimated MTTs, regardless of the modelling approach or distribution used, are positively associated with the percent wetland area and negatively with mean slope in the catchments. In this landscape, low-gradient features such as wetlands increase catchment scale water storage when antecedent conditions are dryer and decrease transit times when there is a moisture surplus, which plausibly explains the increases in MTTs and mean annual runoff from catchments with significant coverage of these landscape features.     

Changes in hillslope hydrology in a perched,shallow soil system due to clearcutting and residual biomass removal

Sustainable fuels legislation and volatility in energy prices have put additional pressures on the forestry sector to intensify the harvesting of biomass for “advanced biofuel” production. To better understand how residual biomass removal after harvest affects forest hydrology in relatively low slope terrain, a Before-After-Control-Impact (BACI) study was conducted in the USDA Forest Service's Marcell Experimental Forest, Minnesota, USA. Hydrological measurements were made from 2010–2013 on a forested hillslope that was divided into three treatment blocks, where one block was harvested and residual biomass removed (Biomass Removed), the second was harvested and residual biomass left (Biomass Left), and the last block was left as an Unharvested Control. The pre-harvest period (2 years) was 2010–11 and post-harvest (2 years) was 2012–13. Water table elevation at the upslope and downslope position, subsurface runoff, and soil moisture were measured between May–November. Mixed effect statistical models were used to compare both the before-after and “control” treatment ratios (ratios between harvested hillslopes and the Unharvested Control hillslope). Subsurface runoff significantly increased (p < .05) at both harvested hillslopes but to a greater degree on the Biomass Left hillslope. Greater subsurface runoff volumes at both harvested hillslopes were driven by substantial increases during fall, with additional significant increases during summer on the Biomass Left hillslope. The hydrological connectivity, inferred from event runoff ratios, increased due to harvesting at both hillslopes but only significantly on the Biomass Left hillslope. The winter harvest minimized soil disturbance, resulting in no change to the effective hydraulic conductivity distribution with depth. Thus, the observed hydrological changes were driven by increased effective precipitation and decreased evapotranspiration, increasing the duration that both harvested hillslopes were hydrologically active. The harvesting of residual biomass appears to lessen hydrological connectivity relative to leaving residual biomass on the hillslope, potentially decreasing downstream hydrological impacts of similar forestry operations.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

Geostatistical features of streambed vertical hydraulic conductivities in Frenchman Creek Watershed in Western Nebraska

This study evaluated the spatial variability of streambed vertical hydraulic conductivity (K_v) in different stream morphologies in the Frenchman Creek Watershed, Western Nebraska, using different variogram models. Streambed K_v values were determined in situ using permeameter tests at 10 sites in Frenchman, Stinking Water and Spring Creeks during the dry season at baseflow conditions. Measurements were taken both in straight and meandering stream channels during a 5 day period at similar flow conditions. Each test site comprised of at least three transects and each transect comprised of at least three K_v measurements. Linear, Gaussian, exponential and spherical variogram models were used with Kriging gridding method for the 10 sites. As a goodness-of-fit statistic for the variogram models, cross-validation results showed differences in the median absolute deviation and the standard deviation of the cross-validation residuals. Results show that using the geometric means of the 10 sites for gridding performs better than using either all the K_v values from the 93 permeameter tests or 10 K_v values from the middle transects and centre permeameters. Incorporating both the spatial variability and the uncertainty involved in the measurement at a reach segment can yield more accurate grid results that can be useful in calibrating K_v at watershed or sub-watershed scales in distributed hydrological models.     

A land-use mix allocation model considering adjacency,intensity, and proximity

ABSTRACT

Land-Use Mix (LUM) refers to the strategy of integrating complementary functions within a building or area. While LUM has become a dominant approach in urban planning, its actual benefits and vision for spatial planning remain unclear. To clarify this issue, this study discerns the spatial features of land-use patterns depending on the compatibilities among land-use categories. Accordingly, this study introduces three LUM measures – adjacency, intensity, and proximity – to identify differences in the spatial distribution of land-use categories. Based on these measures, a land-use allocation model is developed to specify spatial patterns satisfying the given compatibilities. This model is tested by applying the concept of the neighborhood unit on a case study of normative land-use patterns subject to specified compatibilities. The results describe spatial features of four compatibility sets, including a set exhibiting a compatibility conflict between the same land-use pair and LUM measures when, for example, a given land-use pair is compatible in terms of intensity but incompatible in terms of proximity. Understanding the spatial features of a normative land-use pattern that satisfies various possible compatibilities will facilitate the incorporation of the LUM approach into local planning guidance and zoning ordinances.     

Identification of mercurian volcanism: Resolution effects and implications for MESSENGER

Abstract— The possibility of volcanism on Mercury has been a topic of discussion since Mariner 10 returned images of half the planet's surface showing widespread plains material. These plains could be volcanic or lobate crater ejecta. An assessment of the mechanics of the ascent and eruption of magma shows that it is possible to have widespread volcanism, no volcanism on the surface whatsoever, or some range in between. It is difficult to distinguish between a lava flow and lobate crater ejecta based on morphology and morphometry. No definite volcanic features have been identified on Mercury. However, known lunar volcanic features cannot be identified in images with similar resolutions and viewing geometries as the Mariner 10 dataset. Examination of high‐resolution, low Sun angle Mariner 10 images reveals several features which are interpreted to be flow fronts; it is unclear if these are volcanic flows or ejecta flows. This analysis implies that a clear assessment of volcanism on Mercury must wait for better data. MESSENGER (MErcury: Surface, Space ENvironment, GEochemistry, Ranging) will take images with viewing geometries and resolutions appropriate for the identification of such features.