National Floodplain Mapping: Datasets and Methods – 160,000 km in 12 months

This paper describes two projects requiring production of national floodplain maps for England and Wales – some 80,000 km of river. The novel solutions developed have brought together a national Digital Elevation Model (DEM), automatically-generated peak flow estimates at intervals along the watercourses and two alternative methods of calculating the outlines: normal depth calculation; and a purpose-built 2-dimensional raster-based floodplain model, JFLOW. The DEM was derived using Interferometric Synthetic Aperture Radar (IFSAR) techniques and has a vertical precision of ±0.5 m–1.0 m (RMSE) and a 5 m horizontal resolution. The flow estimates were derived by automating Flood Estimation Handbook (FEH) techniques. The normal depth calculations are applied at a number of discrete cross-sections with linear interpolation between to form a 3-dimensional water surface. This is overlain on the DEM to produce the flood outline. Careful manual checking is required at a number of stages. The JFLOW model is based on a discretised form of the 2-dimensional diffusive wave equation and directly simulates the flood outline in a series of overlapping short (1 km) reaches. Flood outlines from the overlapping reaches are merged to produce the overall flood envelope. The model has been written to work as a screen-saver, allowing distributed processing across all computers in an office and manual intervention is minimal. In simple valley situations both methods give similar results, but show differences in more complex areas. Each has advantages and disadvantages, but both have been shown to be a practicable solution to allow production of 160,000 km of flood outline in 12 months.     

Products of an effusive subglacial rhyolite eruption: Bláhnúkur, Torfajökull, Iceland

We present field observations from Bláhnúkur, a small volume (<0.1 km³) subglacial rhyolite edifice at the Torfajökull central volcano, south-central Iceland. Bláhnúkur was probably emplaced during the last glacial period (ca. 115–11 ka). The characteristics of the deposits suggest that they were formed by an effusive eruption in an exclusively subglacial environment, beneath a glacier >400 m thick. Lithofacies associations attest to complex patterns of volcano-ice interaction. Erosive channels at the base of the subglacial sequence are filled by both eruption-derived material and subglacial till, which show evidence for deposition by flowing meltwater. This suggests that meltwater was able to drain away from the vent area during the eruption. Much of the subglacial volcanic deposits consist of conical-to-irregularly shaped lava lobes typically 5–10 m long, set in poorly sorted breccias with an ash-grade matrix. A gradational lavabreccia contact at the base of lava lobes represents a fossilised fragmentation interface, driven by magma-water interaction as the lava flowed over poorly consolidated, waterlogged debris. Sets of columnar joints on the upper surfaces of lobes are interpreted as ice-contact features. The morphology of the lobes suggests that they chilled within conically shaped subglacial cavities 2–5 m high. Avalanche deposits mantling the flanks of Bláhnúkur appear to have been generated by the collapse of lava lobes and surrounding breccia. A variety of deposit characteristics suggests that this occurred both prior to and after quenching of the lava lobes. Collapse events may have occurred when the supporting ice walls were melted back from around the cooling lava lobes and breccias. Much larger lava flows were emplaced in the latter stages of the eruption. Columnar joint patterns suggest that these flowed and chilled within subglacial cavities 20 m high and 100–200 m in length. There is little evidence for magma-water interaction at lava flow margins which suggests that these larger cavities were drained of meltwater. As rhyolite magma rose to the base of the glacier, the nature of the subglacial cavity system played an important role in governing the style of eruption and the volcanic facies generated. We present evidence that the cavity system evolved during the eruption, reflecting variations in both melting rate and edifice growth that are best explained by a fluctuating eruption rate.     

Electrical resistivity ground imaging (ERGI): a new tool for mapping the lithology and geometry of channel-belts and valley-fills

Efforts to map the lithology and geometry of sand and gravel channel‐belts and valley‐fills are limited by an inability to easily obtain information about the shallow subsurface. Until recently, boreholes were the only method available to obtain this information; however, borehole programmes are costly, time consuming and always leave in doubt the stratigraphic connection between and beyond the boreholes. Although standard shallow geophysical techniques such as ground‐penetrating radar (GPR) and shallow seismic can rapidly obtain subsurface data with high horizontal resolution, they only function well under select conditions. Electrical resistivity ground imaging (ERGI) is a recently developed shallow geophysical technique that rapidly produces high‐resolution profiles of the shallow subsurface under most field conditions. ERGI uses measurements of the ground's resistance to an electrical current to develop a two‐dimensional model of the shallow subsurface (<200 m) called an ERGI profile. ERGI measurements work equally well in resistive sediments (‘clean’ sand and gravel) and in conductive sediments (silt and clay). This paper tests the effectiveness of ERGI in mapping the lithology and geometry of buried fluvial deposits. ERGI surveys are presented from two channel‐fills and two valley‐fills. ERGI profiles are compared with lithostratigraphic profiles from borehole logs, sediment cores, wireline logs or GPR. Depth, width and lithology of sand and gravel channel‐fills and adjacent sediments can be accurately detected and delineated from the ERGI profiles, even when buried beneath 1–20 m of silt/clay.     

Role of initial depth at basin margins in sequence architecture: field examples and computer models

A delay in the onset of sedimentation during fault‐related subsidence at a basin margin can occur in both extensional settings, where footwall tilting may cause a diversion of drainage patterns, and in strike‐slip basins, where a source area may be translated along the basin margin. The ‘initial depth’ created by this delay acts as pre‐depositional accommodation and is a partly independent variable. It controls the geometry of the first stratal units deposited at the basin margin and thus modifies the response of the depositional system to subsequent, syndepositional changes in accommodation. In systems with a sharp break in the depositional profile, such as the topset edge in coarse‐grained deltas, the initial depth controls the foreset height and therefore the progradational distance of the topset edge. The topset length, in turn, influences topset accommodation during cyclical base level variations and therefore is reflected in the resulting stacking patterns at both long‐ and short‐term time scales. In the simplified cases modelled in this study, it is the relationship between the initial depth and the net increase in depth over the interval of a relative sea‐level cycle (ΔH) that governs long‐ and short‐term stacking patterns. In situations where the initial depth is significantly larger than ΔH, the topset accommodation of the first delta is insufficient to contain the volume of sediment of younger sequences formed during subsequent relative sea‐level cycles. Therefore, the depositional system tends to prograde over a number of relative sea‐level cycles before the topset area increases so that the long‐term stacking pattern changes to aggradation. Stacking patterns of high‐frequency sequences are influenced by a combination of topset accommodation available and position of the short‐term relative sea‐level cycles on the rising or falling limb of a long‐term sea‐level curve. This determines whether deposits of short‐term cycles are accommodated in delta topsets or foresets, or in both. Variations in stacking pattern caused by different initial depths may be misinterpreted as due to relative sea level or sediment supply changes and it is necessary to consider initial bathymetry in modelling and interpretation of stacking patterns, especially in fault‐bounded basins.     

Petrophysical characterization and natural fracturing in an olivine-dolerite aquifer

As part of a study investigating the naturally-occurring fractures in mafic rocks, two holes were drilled 450 m apart through the Palisades dolerite sill in New York. Well-2 is 229 m deep and Well-3 was drilled to 305 m, both penetrating through the sill and into the underlying Triassic sediments of the Newark Basin. Both holes were logged with downhole geophysical tools, including the BHTV, which acoustically images fractures intersecting the well. Understanding the fracture pattern, density, and porosity in the sill is essential for identifying possible zones of active fluid flow and high permeability. Using the BHTV logs, 96 and 203 fractures were digitally mapped within the sill in Well-2 and Well-3, respectively. Most fractures appear to dip steeply (76-78°). There is a shift in fracture orientation, however, and these fractures may or may not be continuous over the short lateral distance between Well-2 and Well-3. The lithology of the sill as identified by drill chips is nevertheless continuous between the holes. Both intersect a 7 m thick olivine-rich layer about 15 m above the bottom of the sill. Several fractures identified in Well-2 have large apparent aperture (>6cm) which correspond to high porosity zones (6-14%) observed in the logs. Resistivity logs were used to compute porosity using Archie's law and match well with the neutron porosity log in Well-2. We use the relationship between porosity and fracture aperture within the sill at Well-2 to infer the porosity in Well-3. High-porosity, large-aperture zones, including the target olivine layer, are identified in both holes. Changes in the temperature gradient log indicate active fluid flow in the sill, although flow appears to be most active in the sediments. Direct field measurements of bulk permeability, hydrologic modeling of fluid flow and calibration of fracture and log porosity will be undertaken in the future.     

Metal accumulation in four molluscs

Specimens of mussels (Mytilus edulis), winkles (Littorina littorea), dogwhelks (Nucella lapillus) and limpets (Patella vulgata) were collected from the north-east coast of England in 1980–1981 in order to analyse the dried whole soft tissue for zinc, copper and iron. In each case, only a narrow size range of organisms was collected from a small environmentally homogeneous site and sex was taken into account. Most of the frequency distributions of metal concentrations obtained were found to show significant positive skewness. This skewness could not be explained either by the degree of mobility or type of feeding mechanism of the organisms or by any other known parameter, e.g. size, sex, habitat variation. Some individual mussels were found to be extremely good regulators of their zinc concentrations. On the other hand, other individual mussels had the ability to accumulate unusually high levels of zinc.     

The stable-beam seismic modelling method1

The stable-beam method for forward modelling of seismic data is introduced. The method is applicable to geometries which may be approximated by a series of single-valued (in depth) interfaces separating constant-velocity layers. For models of this restricted type, the results are of similar accuracy to those from waveequation-based methods whilst the run times are similar to, or better than, those for simple ray-tracing approaches. The basis of the method is to approximate interfaces by a series of straight-line segments. This allows very rapid and stable ray tracing through the model. Pseudodiffractions are then added from all of the interface discontinuities formed between adjacent segments. These pseudodiffractions have the effect of correcting for the phase, amplitude and wavefront continuity errors introduced by the interface approximation. Comparison of the stable-beam results to analytical, Kirchhoff, finite-difference and physical model results confirm the accuracy ofthe technique.     

Microstrain stability of Peninsular India 1864–1994

We report the results of the South Indian Strain Measuring Experiment (SISME) designed to determine whether strain related to microseismicity in the past century may have deformed the networks of the 19th century Great Trigonometrical Survey of India (GTS). More than a dozen GTS points were measured between Mangalore, Madras, and Kanyakumari in southernmost India using GPS geodesy to determine regional deformation. Detailed measurements were made near two of the original baselines of the survey to determine the reliability of dilatational strain data for the network. The regional measurements revealed negligible regional dilatational (+ 11.2 + 10 microstrain) and shear strain changes (0.66± 1.2μradians) in the southernmost 530 km of India. In addition to these measurements, we determined the rate of northward and eastward motion of a point in Bangalore (1991–1994) in the ITRF92 reference frame to be 39 ± 3.5 mm/year, and 51 ± 11 mm/year respectively. This is consistent with NUVEL-1A plate motion estimate for India. Simultaneous measurements to a point near Kathmandu reveal that the Indian plate and the Southern Himalaya are moving approximately in unison, placing an upper limit on the rate of creep processes beneath the lesser Himalaya of ≈6 mm/year, and suggesting relatively rigid behavior of the Indian plate north of Bangalore. The stability of the Indian plate is confirmed by the absence of significant changes in the lengths of the two baselines at Bangalore and Cape Comorin, which, within the limits of experimental error have not changed since 1869. The measurements place an upper limit for recent deformation in the southern peninsula, and hence a lower limit for the renewal time for intraplate earthquakes in the region of approximately 10,000 years, assuming shear failure strain of approximately 100 μradians. This, in turn, implies that recurrence intervals for Peninsular Earthquakes far exceed the length of the written historic record, suggesting that the characterisation of seismic recurrence intervals from historical studies is likely to be fruitless. In contrast, the SISME experiment demonstrates that the noise level of geodetic studies based on 19th century GTS data is less than 0.02 μstrain/year, providing considerable scope for delineating regions of anomalously high seismogenic strain, by GPS measurements at all available trig points of the 19th century GTS survey.