Hydrology and related changes after harvesting native forest catchments and establishing pinus radiata plantations. Part 2. The native forest water balance and changes in streamflow after harvesting

Six small, steep, south-west facing catchments (1.63–4.62 ha) have been monitored in Westland, New Zealand since 1974. Two catchments were retained in native mixed evergreen forest and the rest were subjected to various harvesting and land preparation techniques before being planted with Pinus radiata between 1977 and 1980. The 11-year water balance for the native forest catchments was: rain = streamflow + interception loss + transpiration + seepage (2370mm = 1290mm + 620mm + 360mm + 100mm). In the year after treatment streamflow generally increased by 200–250 mm, except for one treatment (clearfelling, herbicide application, no riparian reserve) where the increase was 550 mm. The catchments were planted with Pinus radiata, but rapid colonization by bracken (Pteridium esculentum) and Himalayan honeysuckle (Leycesteria formosa) led to a rapid decline in streamflow, which returned to pre-treatment levels after an average of about five years. Streamflow yields then continued to decline for another two to three years before stabilizing at a level about 250mm yr^?1 lower than pre-treatment levels. At this time the catchments had a dense bracken/honeysuckle understorey beneath 5 m tall pine trees.     

The identification of runoff-production mechanisms using environmental isotopes in a tussock grassland catchment,eastern otago,New Zealand

A previous hydrometric study of runoff production in tussock grassland drainage basins in Otago (45°50′S, 169°45′E), New Zealand, revealed a marked change of slope in storm hydrograph recessions. An environmental isotope study was initiated to investigate the runoff mechanisms operating and to test specific hypotheses to explain this break in the hydrograph recession. The results indicated that for quickflow volumes in excess of 10mm, the first part of the storm hydrograph can be attributed to two separate sources, namely, ‘old’ water from a shallow, unconfined groundwater reservoir and ‘new’ water from saturation overland flow on the lower wetlands of concave slopes. Despite the extensive area of wetlands, ‘old’ water runoff from the unconfined groundwater reservoir is delivered more rapidly to the stream than ‘new’ water from saturation overland flow. Substantial surface storage in the wetlands has first to be exceeded before rain becomes a significant part of stream discharge. For quickflow volumes less than 10mm, only ‘old’ water from groundwater contributes to the first part of the hydrograph recession. This means that only the largest 7 per cent of storms (in terms of quickflow volume) generate quickflow containing significant amounts of ‘new water’. The second part of the recession of the storm hydrograph consists of ‘old’ water derived from a remarkably well-mixed shallow unconfined groundwater body.     

The spectacular Supernova SN1987A: Grain evaporation and condensation

Non-similarity solutions of the equations governing the motion of a perfect gas behind a cylindrical shock wave of variable strength have been obtained. These solutions are applicable to both the weak and the strong shocks. The nature of flow and field variables are illustrated through graphs. The total energy of the wave is taken to be constant.     

Complex P waveforms from a Gulf of Aden earthquake

Summary. It is shown that complex teleseismic P waveforms from a shallow earthquake in a tectonic area can be interpreted using a simple source model embedded in a plane layer velocity structure (with sea layer) whose details are based upon independent evidence. This gives hope that structural complexity in tectonic regions may not always make distant P -wave seismograms impossible to interpret, and that, instead, source complexity may be responsible for some of the many complex waveforms observed, even for earthquakes below magnitude m _b 5.5.     

Fault plane solutions using relative amplitudes of P and pP

Summary. One way of finding the fault plane orientations of small shallow earthquakes is by the generation of theoretical P -wave seismograms to match those observed at several distant stations. Here, a technique for determining the uniqueness of fault plane solutions computed using the modelling method of Douglas et al . is described. Relative amplitudes of P and pP , and their polarities if unambiguous, are measured on the observed seismograms to be modelled, and appropriate confidence limits are assigned to each measurement. A systematic search is then made for all fault plane orientations which satisfy these observations.
Examples show that if P and pP are not severely contaminated by other arrivals, a well-defined and unique fault plane orientation can often be computed using as few as three stations well distributed in azimuth. Further, even if pP is not identifiable on a particular seismogram, then an upper bound on its amplitude – deduced from the observed coda – still places a significantly greater constraint on the fault plane orientation than would be provided by a P onset polarity alone. Modelling takes account of all such information, and is able to further eliminate incompatible solutions (e.g. by the correct simulation of sP ). It follows that if solutions can be found which satisfy many observed seismograms, this places high significance on the validity of the assumed double-couple source mechanism.
This relative amplitude technique is contrasted with the familiar first motion method of fault plane determination which requires many polarity readings, whose reliabilities are difficult to quantify. It is also shown that fault plane orientations can be determined for earthquakes below the magnitude at which first motion solutions become unreliable or impossible.     

Oil and fishing — conflict in the North Sea

A programme of research is underway at Aberdeen University on the conflicts arising in the North Sea between the fishing and offshore oil and gas industries. This research note summarizes the work that has been undertaken to date and the main aspects on which future work will concentrate.     

Plate boundary deformation and continuing deflation of the Askja volcano,North Iceland,determined with GPS, 1987–1993

GPS geodetic measurements were conducted around the Askja central volcano located at the divergent plate boundary in north Iceland in 1987, 1990, 1992 and 1993. The accuracy of the 1987 and 1990 measurements is in the range of 10 mm for horizontal components; the accuracy of the 1992 and 1993 measurements is about 4 mm in the horizontal plane. Regional deformation in the Askja region is dominated by extension. Points located outside a 30–45 km wide plate boundary deformation zone indicate a displacement of 2.4±0.5 cm/a in the direction N 99°E±12° of the Eurasian plate relative to the North American plate in the period 1987–1990. Within the plate boundary deformation zone extensional strain accumulates at a rate of 0.8 strain/a. Displacement of control points next to Askja (>7 km from the caldera center) in the periods 1990–1993 and 1992–1993 show deflation and contraction towards the caldera. These results are in accordance with the results obtained by other geodetic methods in the area, which indicate that the deflation at Askja occurs in response to a pressure decrease at about 2.8 km depth, located close to the center of the main Askja caldera. A Mogi point source was fixed at this location and the GPS data used to solve for the source strength. A central subsidence of 11±2.5 cm in the period 1990–1993 is indicated, and 5.5±1.5 cm in the period 1992–1993. The maximum tensional strain rate, according to the point source model, occurs at a horizontal distance of 2.5–6 km from the source, at the same location as the main caldera boundary. Discrepancies between the observed displacements and predicted displacements from the Mogi model near the Askja caldera can be attributed to the regional eastwest extension that occurs at Askja.     

Geochemical fingerprinting of oceanic basalts with applications to ophiolite classification and the search for Archean oceanic crust

Two geochemical proxies are particularly important for the identification and classification of oceanic basalts: the Th–Nb proxy for crustal input and hence for demonstrating an oceanic, non-subduction setting; and the Ti–Yb proxy for melting depth and hence for indicating mantle temperature and thickness of the conductive lithosphere. For the Th–Nb proxy, a Th/Yb–Nb/Yb projection demonstrates that almost all oceanic basalts lie within a diagonal MORB–OIB array with a principal axis of dispersion along the array. However, basalts erupted at continental margins and in subduction zones are commonly displaced above the MORB–OIB array and/or belong to suites with principal dispersion axes which are oblique to the array. Modelling of magma–crust interaction quantifies the sensitivity of the Th–Nb proxy to process and to magma and crustal compositions. For the Ti–Yb proxy, the equivalent Ti/Yb–Nb/Yb projection features a discriminant boundary between low Ti/Yb MORB and high Ti/Yb OIB that runs almost parallel to the Nb/Yb axis, reflecting the fact that OIB originate by melting beneath thicker lithosphere and hence by less melting and with residual garnet. In the case of volcanic-rifted margins and oceanic plume–ridge interactions (PRI), where hot mantle flows toward progressively thinner lithosphere (often becoming more depleted in the process), basalts follow diagonal trends from the OIB to the MORB field. Modelling of mantle melting quantifies the sensitivity of the Ti–Nb proxy to mantle potential temperature and lithospheric thickness and hence defines the petrogenetic basis by which magmas plot in the OIB or MORB fields. Oceanic plateau basalts lie mostly in the centre of the MORB part of that field, reflecting a high degree of melting of fertile mantle. Application of the proxies to some examples of MORB ophiolites helps them to be further classified as C (contaminated)-MORB, N (normal)-MORB, E (enriched)-MORB and P (plume)-MORB ophiolites, which may add a useful dimension to ophiolite classification. In the Archean, the hotter magmas, higher crustal geotherms and higher Th contents of contaminants all result in widespread crustal input that is easy to detect geochemically with the Th–Nb proxy. Application of this proxy to Archean greenstones demonstrates that almost all exhibit a crustal component even when reputedly oceanic. This indicates, either that some interpretations need to be re-examined or that intra-oceanic crustal input is important in the Archean making the proxy less effective in distinguishing oceanic from continental settings. The Ti–Yb proxy is not effective for fingerprinting Archean settings because higher mantle potential temperatures mean that lithospheric thickness is no longer the critical variable in determining the presence or absence of residual garnet.