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.     

The brecciated texture of polymict eucrites: Petrographic investigations of unequilibrated meteorites from the Antarctic Yamato collection

We report on the petrography and mineralogy of five Yamato polymict eucrites to better constrain the formation and alteration of crustal material on differentiated asteroids. Each sample consists of different lithic clasts that altogether form four dominant textures and therefore appear to originate from closely related petrological areas within Vesta′s crust. The textures range from subophitic to brecciated, porphyritic, and quench‐textured, that differ from section to section. Comparison with literature data for these samples is therefore difficult, which stresses that polymict eucrites are extremely complex in their petrography and investigation of only one thick section may not be representative for the host rock. We also show that sample Y‐793548 consists of more than one lithic unit and must therefore be classified as polymict instead of monomict. The variety and nature of lithic textures in the investigated Yamato meteorites indicate shock events, intense post‐magmatic thermal annealing, and secondary alteration. These postmagmatic features occur in different intensities, varying from clast to clast or among coexisting mineral fragments on a small, local scale. Several clasts within the eucrites studied have been modified by late‐stage alteration processes that caused deposition of Fe‐rich olivine and Fe enrichment along cracks crosscutting pyroxene crystals. However, formation of these secondary phases seems to be independent of the degree of thermal metamorphism observed within every type of clast, which would support a late‐stage metasomatism model for their formation.     

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.     

Symplectic Mappings of Co-orbital Motion in the Restricted Problem of Three Bodies

The dynamics of co-orbital motion in the restricted three-body problem are investigated by symplectic mappings. Analytical and semi-numerical mappings have been developed and studied in detail. The mappings have been tested by numerical integration of the equations of motion. These mappings have been proved to be useful for a quick determination of the phase space structure reflecting the main characteristics of the dynamics of the co-orbital problem.     

Scanning electron microscopy,cathodoluminescence, and Raman spectroscopy of experimentally shock‐metamorphosed quartzite

Abstract— We studied unshocked and experimentally (at 12, 25, and 28 GPa, with 25, 100, 450, and 750°C pre‐shock temperatures) shock‐metamorphosed Hospital Hill quartzite from South Africa using cathodoluminescence (CL) images and spectroscopy and Raman spectroscopy to document systematic pressure or temperature‐related effects that could be used in shock barometry. In general, CL images of all samples show CL‐bright luminescent patchy areas and bands in otherwise nonluminescent quartz, as well as CL‐dark irregular fractures. Fluid inclusions appear dominant in CL images of the 25 GPa sample shocked at 750°C and of the 28 GPa sample shocked at 450°C. Only the optical image of our 28 GPa sample shocked at 25°C exhibits distinct planar deformation features (PDFs). Cathodoluminescence spectra of unshocked and experimentally shocked samples show broad bands in the near‐ultraviolet range and the visible light range at all shock stages, indicating the presence of defect centers on, e.g., SiO₄ groups. No systematic change in the appearance of the CL images was obvious, but the CL spectra do show changes between the shock stages. The Raman spectra are characteristic for quartz in the unshocked and 12 GPa samples. In the 25 and 28 GPa samples, broad bands indicate the presence of glassy SiO₂, while high‐pressure polymorphs are not detected. Apparently, some of the CL and Raman spectral properties can be used in shock barometry.     

Optical properties of solid and liquid sulfur at visible and infrared wavelengths

Knowledge of the optical constants of elemental sulfur has potential applications to Venus, Jupiter, Io, Amalthea, and the Earth. The real part, n, of the index of refraction of liquid sulfur (at 133°C) and of solid orthorhombic sulfur (at 25°C) for the wavelength range 0.4–2.0 μm were measured ellipsometrically. The imaginary part, k, of the refractive index of liquid sulfur was obtained by transmittance measurements at the same temperature and wavelength range. The reflectance of semi-infinite slabs of solid and liquid sulfur is calculated using the measured n and k values. We confirm that sulfur melts on Io would be classified as “black” by the Voyager imaging system.     

Recent extreme slope failures in glacial environments: effects of thermal perturbation

This paper describes recent exceptional slope failures in high-mountain, glacial environments: the 2002 Kolka–Karmadon rock–ice avalanche in the Caucasus, a series of ice–rock avalanches on Iliamna Volcano, Alaska, the 2005 Mt. Steller rock–ice avalanche in Alaska, and ice and rock avalanches at Monte Rosa, Italy in 2005 and 2007. Deposit volumes range from 10⁶ to 10⁸ m³ and include rock, ice and snow. Here we focus on thermal aspects of these failures reflecting the involvement of glacier ice and permafrost at all sites, suggesting that thermal perturbations likely contributed to the slope failures. We use surface and troposphere air temperatures, near-surface rock temperatures, satellite thermal data, and recent 2D and 3D thermal modeling studies to document thermal conditions at the landslide sites. We distinguish between thermal perturbations of volcanic-geothermal and climatic origin, and thermal perturbations related to glacier–permafrost interaction. The data and analysis support the view that recent, current and future climatic change increases the likelihood of large slope failures in steep glacierized and permafrost terrain. However, some important aspects of these settings such as the geology and tectonic environment remain poorly understood, making the identification of future sites of large slope instabilities difficult. In view of the potentially large natural disasters that can be caused by such slope failures, improved data and understanding are needed.