Effects of glacial retreat on proglacial streams and riparian zones in the Coast and North Cascade Mountains

Proglacial stream development was studied in coastal British Columbia and Washington, focusing on reaches exposed by post‐Little Ice Age (LIA) glacier retreat, to address three principal questions: (i) Does the legacy of LIA glaciation influence the evolution of channel morphology? (ii) How long does it take for riparian forest to establish following glacier retreat? (iii) Can newly exposed proglacial streams provide suitable fish habitat? Channel morphologies were identified by field surveys of 69 reaches in 10 catchments. Riparian forest development and potential fish habitat were characterized in those reaches and an additional 22 catchments using GIS analysis. The landscape template imposed by the Quaternary glaciation appears to override most of the modern effects of the LIA in controlling channel‐reach morphology. Binary logistic regression analysis identified elevation and time since deglaciation as primary controls on the presence of riparian forest. At higher elevations, establishment of morphologically functional riparian forest could take several centuries, prolonged by channel instability associated with post‐LIA sediment inputs. Of the recently deglaciated streams included in this analysis, the majority (86%) of the total length was of suitable gradient for fish and could be accessed either by downstream populations or from adjacent lakes. Predicted maximum weekly average stream temperature (MWAT) indicated that the post‐LIA study streams were thermally suitable for cold‐water fish. A future scenario of glacier loss would cause a 14% decline in accessible cold‐water thermal habitat in post‐LIA streams. Decreased summer flows due to glacier retreat could further limit usable habitat by reducing stream depths and wetted perimeters. Copyright © 2013 John Wiley & Sons, Ltd.     

Landslide inventory in a rugged forested watershed: a comparison between air-photo and field survey data

Landslide inventories are routinely compiled by means of aerial photo interpretation (API). When examining photo pairs, the forest canopy (notably in old-growth forest) hides a population of “not visible” landslides. In the present study, we attempt to estimate how important is the contribution of landslides not detectable from aerial photographs to the global mass of sediment production from mass failures on forested terrain of the Capilano basin, coastal British Columbia. API was coupled with intensive fieldwork for identification and measurement of all landslides. A 30-year framework was adopted. We show that “not visible” landslides can represent up to 85% of the total number of failures and account for 30% of the volume of debris mobilised. Such percentages display high sub-basin variability with rates of sediment production varying by one order of magnitude between two sub-basins of the study area. This is explained qualitatively by GIS-based analysis of slope frequency distributions, drainage density, and spatial distribution of surficial materials. Such observations find further support in the definitions of transport-limited and supply-limited basins. As a practical consideration to land managers, we envisage that supplementary fieldwork for landslide identification is mandatory in transport-limited systems only. Fieldwork has demonstrated that gully-related failures have a greater importance than one could expect from API.     

Impact of wastewater discharge on the channel morphology of ephemeral streams

The impact of wastewater flow on the channel bed morphology was evaluated in four ephemeral streams in Israel and the Palestinian Territories: Nahal Og, Nahal Kidron, Nahal Qeult and Nahal Hebron. Channel changes before, during and after the halting of wastewater flow were monitored. The wastewater flow causes a shift from a dry ephemeral channel with intermittent floods to a continuous flow pattern similar to that of humid areas. Within a few months, nutrient‐rich wastewater flow leads to rapid development of vegetation along channel and bars. The colonization of part of the active channel by vegetation increases flow resistance as well as bank and bed stability, and limits sediment availability from bars and other sediment stores along the channels. In some cases the established vegetation covers the entire channel width and halts the transport of bed material along the channel. During low and medium size flood events, bars remain stable and the vegetation intact. Extreme events destroy the vegetation and activate the bars. The wastewater flow results in the development of new small bars, which are usually destroyed by flood flows. Due to the vegetation establishment, the active channel width decreases by up to 700 per cent. The deposition of fine sediment and organic material changed the sediment texture within the stable bar surface and the whole bed surface texture in Nahal Hebron. The recovery of Nahal Og after the halting of the wastewater flow was relatively fast; within two flood seasons the channel almost returned to pre‐wastewater characteristics. The results of the study could be used to indicate what would happen if wastewater flows were introduced along natural desert streams. Also, the results could be used to predict the consequences of vegetation removal as a result of human intervention within the active channel of humid streams. Copyright © 2001 John Wiley & Sons, Ltd.     

Extension of the Henry semi-analytical solution for saltwater intrusion in stratified domains

The Henry semi-analytical solution is developed for stratified aquifers with exponential and Gaussian permeability–depth relationships and small dispersion. The semi-analytical solution is sought by expanding the stream function and the concentration in an infinite Fourier series truncated at given orders. Due to the heterogeneity, the semi-analytical solution contains additional terms and quickly becomes impractical because of the high truncation orders required to yield stable results. An appropriate evaluation of the most expensive term, involving fairly complex summations, is proposed to render the computation of the semi-analytical solution affordable. Three configurations with the same transmissivity as that for the homogeneous Henry problem are investigated. The first configuration considers an exponential decay of the permeability (k) with depth. The second and third configurations consider an exponential increase and a Gaussian high-low-high distribution of k with depth, respectively. The three configurations are also investigated numerically using an accurate numerical code based on the method of lines and advanced spatial discretization schemes. An excellent agreement is obtained between the semi-analytical and the numerical solutions for all test cases. The results show that the amount of saltwater intrusion is strongly dependent on the heterogeneity. Further, the heterogeneous Henry problem is found to be more suitable for benchmarking density-driven numerical codes because the numerical solution is more sensitive to the grid size than for the homogeneous problem.     

Specimens size effect D/H on cyclic behaviour and liquefaction potential of clean sand

Acta Geotechnica - Owing to the lack of complementary shear stress on vertical boundaries of a direct simple shear (DSS) specimen, the stress/strain distributions in soil specimens are not uniform....     

Debris flow initiation and sediment recharge in gullies

Landslides that enter gullied low-order drainages can either initiate debris flow or stop, depositing sediment in the channel. This process is one of the most common ways that debris flows initiate, but little attention to date has been paid to evaluating the factors that affect whether or not the initial landslide will become a debris flow or deposit sediment in the channel. Statistically significant parameters that determine whether slope failures become debris flows or act to recharge in-channel sediment are channel gradient, angle of entry of failure into the channel, initial failure volume, and the amount of in-channel stored sediment. Steeper channels, low angles of entry, lower volumes of in-channel sediment, and larger initial failures were more likely to result in debris flows. This study found that as the volume of in-channel stored sediment increased, the volume of initial failure required to initiate a debris flow also increased. This result calls into question the simple supply-limited model of cyclical debris recharge and debris flow in low-order gullied drainages and suggests a negative feedback mechanism between debris accumulation and debris flow susceptibility.     

The structural inventory of a small complex impact crater: Jebel Waqf as Suwwan,Jordan

The investigation of terrestrial impact structures is crucial to gain an in‐depth understanding of impact cratering processes in the solar system. Here, we use the impact structure Jebel Waqf as Suwwan, Jordan, as a representative for crater formation into a layered sedimentary target with contrasting rheology. The complex crater is moderately eroded (300–420 m) with an apparent diameter of 6.1 km and an original rim fault diameter of 7 km. Based on extensive field work, IKONOS imagery, and geophysical surveying we present a novel geological map of the entire crater structure that provides the basis for structural analysis. Parametric scaling indicates that the structural uplift (250–350 m) and the depth of the ring syncline (<200 m) are anomalously low. The very shallow relief of the crater along with a NE vergence of the asymmetric central uplift and the enhanced deformations in the up‐range and down‐range sectors of the annular moat and crater rim suggest that the impact was most likely a very oblique one (~20°). One of the major consequences of the presence of the rheologically anisotropic target was that extensive strata buckling occurred during impact cratering both on the decameter as well as on the hundred‐meter scale. The crater rim is defined by a circumferential normal fault dipping mostly toward the crater. Footwall strata beneath the rim fault are bent‐up in the down‐range sector but appear unaffected in the up‐range sector. The hanging wall displays various synthetic and antithetic rotations in the down‐range sector but always shows antithetic block rotation in the up‐range sector. At greater depth reverse faulting or folding is indicated at the rim indicating that the rim fault was already formed during the excavation stage.     

Mapping groundwater discharge to a coastal lagoon using combined spatial airborne thermal imaging,radon (222Rn) and multiple physicochemical variables

Coastal lagoons are significant wetland environments found on coastlines throughout the world. Groundwater seepage may be a key component of lagoon water balances, though only a few studies have investigated large (>100 km²) coastal lagoons. In this study, we combined airborne thermal infrared imagery with continuous measurements of radon (²²²Rn—a natural groundwater tracer), conductivity, water temperature and dissolved oxygen to map groundwater seepage to a large coastal lagoon in New Zealand. We found evidence of seepage along the margins of the lagoon but not away from the margins. Our findings confirmed previously known seepage zones and identified new potential locations of groundwater inflow. Both point source and diffuse seepage occurred on the western and northwestern margins of the lagoon and parallel to the barrier between the lagoon and sea. These observations imply geologic controls on seepage. The combination of remote sensing and in-situ radon measurements allowed us to effectively map groundwater discharge areas across the entire lagoon. Combined, broad-scale qualitative methods built confidence in our interpretation of groundwater discharge locations in a large, dynamic coastal lagoon.