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.     

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.     

A WKBJ spectral method for computation of SV synthetic seismograms in a cylindrically symmetric medium

Summary. A method of synthetic seismogram computation for teleseismic SV -waves is developed in order to treat quantitatively SV -waves in problems of body wave source inversion and source—receiver structure studies. The method employs WKBJ theory for a generalized ray in a vertically inhomogeneous half-space and the propagator matrix technique for waves in near-surface homogeneous layers. Wavenumber integration is done along the real axis of the wavenumber plane and anelasticity is included by using complex velocity in all regions of the earth model. The near-surface source structure is taken into account in the computation for the case of the shallow source by allowing a point source to be located in the homogeneous layers. Source and receiver area structures are also allowed to differ. A general moment tensor point source is considered.     

The mass ratio of Nova Muscae 1991

We report intermediate resolution H spectroscopy of the black hole candidate Nova Muscae 1991 during quiescence. We classify the companion star as a K3-K4V which contributes 85±6 percent to the total flux from the binary. The photospheric absorption lines are broadened by 106±13 kms⁻¹ with respect to template field stars, leading to a system mass ratio of q =M₁/M₂ = 7.8_−2.0^+3.4. Doppler imaging of the H line shows strong emission coming from the secondary star (EW=3.1±0.6Å) which we associate with chromospheric activity. However, the hot-spot is not detected and this may indicate a lower mass transfer rate than in other X-ray transients of comparable orbital periods. The surface brightness distribution of the accretion disk in H follows a relation I∝R^−1.1, less steep than typically observed in cataclysmic variables. Updated system parameters are also presented.     

Near-Infrared Imaging Spectrophotometry Of Comet C/1995 O1 (Hale-Bopp) Near Perihelion

We present 1- to 5-μm broadband and CVF images of comet Hale-Bopp taken 1997 February 10.5 UT, 50 days before perihelion. All the images exhibit a nonspherical coma with a bright “ridge” in the direction of the dust tail approximately 10″ from the coma. Synthetic aperture spectrophotometry implies that the optically important grains are of a radius ≤0.4 μm; smallest radius for any comet seen to date. The variation of the integrated surface brightness with radial distance from the coma (ρ) in all the images closely follows the “steady state” ρ⁻¹ model for comet dust ablation (Gehrz and Ney, 1992). The near-infrared colors taken along the dust tail are not constant implying the dust grain properties vary with coma distance. This revised version was published online in July 2006 with corrections to the Cover Date.     

The D-CIXS X-ray mapping spectrometer on SMART-1

The D-CIXS Compact X-ray Spectrometer will provide high quality spectroscopic mapping of the Moon, the primary science target of the ESA SMART-1 mission. D-CIXS consists of a high throughput spectrometer, which will perform spatially localised X-ray fluorescence spectroscopy. It will also carry a solar monitor, to provide the direct calibration needed to produce a global map of absolute lunar elemental abundances, the first time this has been done. Thus it will achieve ground breaking science within a resource envelope far smaller than previously thought possible for this type of instrument, by exploiting two new technologies, swept charge devices and micro-structure collimators. The new technology does not require cold running, with its associated overheads to the spacecraft. At the same time it will demonstrate a radically novel approach to building a type of instrument essential for the BepiColombo mission and potential future planetary science targets.