Construction dynamics of a lava channel

We use a kinematic GPS and laser range finder survey of a 200 m-long section of the Muliwai a Pele lava channel (Mauna Ulu, Kilauea) to examine the construction processes and flow dynamics responsible for the channel–levee structure. The levees comprise three packages. The basal package comprises an 80–150 m wide ′a′a flow in which a ∼2 m deep and ∼11 m wide channel became centred. This is capped by a second package of thin (<45 cm thick) sheets of pahoehoe extending no more than 50 m from the channel. The upper-most package comprises localised ′a′a overflows. The channel itself contains two blockages located 130 m apart and composed of levee chunks veneered with overflow lava. The channel was emplaced over 50 h, spanning 30 May–2 June, 1974, with the flow front arriving at our section (4.4 km from the vent) 8 h after the eruption began. The basal ′a′a flow thickness yields effusion rates of 35 m³ s⁻¹ for the opening phase, with the initial flow advancing across the mapped section at ∼10 m/min. Short-lived overflows of fluid pahoehoe then built the levee cap, increasing the apparent channel depth to 4.8 m. There were at least six pulses at 90–420 m³ s⁻¹, causing overflow of limited extent lasting no more than 5 min. Brim-full flow conditions were thus extremely short-lived. During a dominant period of below-bank flow, flow depth was ∼2 m with an effusion rate of ∼35 m³ s⁻¹, consistent with the mean output rate (obtained from the total flow bulk volume) of 23–54 m³ s⁻¹. During pulses, levee chunks were plucked and floated down channel to form blockages. In a final low effusion rate phase, lava ponded behind the lower blockage to form a syn-channel pond that fed ′a′a overflow. After the end of the eruption the roofed-over pond continued to drain through the lower blockage, causing the roof to founder. Drainage emplaced inflated flows on the channel floor below the lower blockage for a further ∼10 h. The complex processes involved in levee–channel construction of this short-lived case show that care must be taken when using channel dimensions to infer flow dynamics. In our case, the full channel depth is not exposed. Instead the channel floor morphology reflects late stage pond filling and drainage rather than true channel-contained flow. Components of the compound levee relate to different flow regimes operating at different times during the eruption and associated with different effusion rates, flow dynamics and time scales. For example, although high effusion rate, brim-full flow was maintained for a small fraction of the channel lifetime, it emplaced a pile of pahoehoe overflow units that account for 60% of the total levee height. We show how time-varying volume flux is an important parameter in controlling channel construction dynamics. Because the complex history of lava delivery to a channel system is recorded by the final channel morphology, time-varying flow dynamics can be determined from the channel morphology. Developing methods for quantifying detailed flux histories for effusive events from the evidence in outcrop is therefore highly valuable. We here achieve this by using high-resolution spatial data for a channel system at Kilauea. This study not only indicates those physical and dynamic characteristics that are typical for basaltic lava flows on Hawaiian volcanoes, but also a methodology that can be widely applied to effusive basaltic eruptions.     

Does quantifying antecedent flow conditions improve stream phosphorus export estimation?

A reliable and economical method for the estimation of nutrient export (e.g. phosphorus) in stream flow from catchments is necessary to quantify the impact of land use or land use change upon aquatic systems. The transport of phosphorus (P) from soil to water is known to impact negatively on water quality. A key observation from studies is that most P export occurs during high stream flow. However, it is not yet clear how flood-antecedent conditions affect the P export during flood events. In this study, the P loss from soil to water as represented by soluble reactive phosphorus (SRP) in stream waters from three different catchments, varying in land use, scale and location in Ireland was monitored over 1 year. This study examined the role of antecedent stream flow conditions on SRP export and identifies a catchment-specific relationship between SRP flood event load (EL) and a flow ratio (FR). The FR is defined as the ratio of the flood event volume (EV) to the pre-event volume (PEV). The latter is the cumulative flow volume for a number of days preceding the event. This PEV period was found to be longer (average 81 days) in the grassland catchments which were known to be saturated with soil P than in the forested catchments (average 21 days) with minimal soil P. This FR ratio is a measure of the antecedent hydrological state (wet or dry) of the catchment. For SRP for each catchment, a specific relationship between SRP EL and FR was identified. The annual SRP export was estimated, using this ratio and compared with the concentration/discharge (C/Q) method. The new flow ratio method was used with data from 12 flood events during the year to estimate an annual export of SRP. For the two grassland catchments in the study, using the FR method, we estimated an SRP export of 1.77 and 0.41 kg ha⁻¹ yr⁻¹. Using the C/Q method, for the same sites, our estimate of SRP export was 1.70 and 0.50 kg ha⁻¹ yr⁻¹ respectively. The C/Q method used SRP concentrations covering 40% of the year while the FR method used only 12 flood events covering less than 2% of the year. This new method which takes account of the antecedent flow state of the river is an alternative to and may be more promising than the traditional C/Q method, particularly when short duration or flood sampling of water quality is carried out.     

Sm-Nd studies of Archaean metasediments and metavolcanics from West Greenland and their implications for the Earth's early history

Metasedimentary and metavolcanic rocks from the Archaean of West Greenland have been examined for evidence of crustal components greater than 3.8 Ga in age and for their compatibility with the presently adopted bulk Earth Sm-Nd parameters. Sm-Nd isotopic data have been obtained for the garbenschiefer metagabbro unit, metasediments from the Isua supracrustal belt, gneisses interior to the Isua belt and metasediments from the Malene supracrustal belt.Using estimates of emplacement age (T) of between 3.77 and 3.67 Ga for the parental volcanics to the garbenschiefer unit, initial¹⁴³Nd/¹⁴⁴Nd ratios yield positiveε_Nd^T values between +1.0 and +3.1 (relative to the CHUR parameters) for seven out of eight samples. Model Sm-Nd ages for the Isua gneisses and metasediments are only compatible with their estimated stratigraphic ages if their sources were ca.+2ε_Nd relative to CHUR at those times. Similarly, model Sm-Nd ages for the Malene samples are only compatible with stratigraphic age constraints when based on a source evolution with positiveε_Nd^T. Implications of these results for the early development of the Earth's mantle are discussed.     

Quantifying the spatial distribution of sediment transport in an experimental gully system using the morphological method

Whilst time-series of sediment transport in gullies in both laboratory experimental and field settings can be determined through instrumentation, quantifying the spatial distribution of transport rates remains challenging. The morphological method, which was proposed for estimating bed-material transport in both one- and two-dimensions in rivers, provides an alternative. Here, we developed this method for gully systems. A laboratory catchment was used to simulate gully erosion. High-resolution topographical data were acquired by close-range digital photogrammetry. Morphological changes were determined using high-resolution topographic data and an associated level of detection. Based on measured morphological changes, one-dimensional (1D) and two-dimensional (2D) sediment transport rates were calculated via cross-section by cross-section routing (1D) and cell by cell routing (2D). The 1D application provided a general trend of longitudinal variation of sediment transport for the whole gully system, increased gradually from zones of headward extension to a zone downstream where erosion and deposition were in balance, and sediment transport rates less variable in space. For the 2D application, hydrological and blended hydrological-hydraulic routing solutions were compared. We found that the level of negative transport was insensitive to whether or not a blended hydrological-hydraulic routing was used and that results from applying the hydrological routing throughout were not significantly degraded. We also found that consideration should be given to spatial and temporal resolution of the topographic data. The 2D application provided spatial patterns of sediment transport that vary with gully evolution. The main gully remained a high transport corridor but branch transport became more important through time. The framework we report provides an additional tool for both experimental and field quantification of the spatial patterns of sediment transport in gullies; and quantification of how these patterns change under different forcing factors.     

Large N-S convergence at the northern edge of the Tibetan plateau? New Early Cretaceous paleomagnetic data from Hexi Corridor, NW China

Nine Early Cretaceous paleomagnetic sites have been collected in the Yumen area of the Hexi Corridor (NW China). Magnetic directions isolated at lower temperatures fail the fold test, and lie close to the geocentric axial dipole field direction before tilt correction. High temperature components are carried by magnetite and/or hematite, all with normal polarity, and pass the fold test. The average paleomagnetic pole from the nine sites is at λ=75.5°N, φ=169.9°E (A₉₅=7.7°). These results are consistent with those from other areas of the North China block (NCB), but significantly different from those from the Qaidam Basin on the southern side of the Qilian Mountains. They suggest that: (1) the Yumen region behaved as a rigid part of the NCB since at least the Early Cretaceous; (2) 740±500 km of north-south directed convergence has taken place between the NCB and Qaidam, within the Qilian Mountains and (3) extrusion of Qaidam was accompanied by a 23±5° relative rotation with respect to North China. This is larger than implied by the maximum left lateral slip on the Altyn Tagh fault system. The same data imply some 1000±800 km of Cenozoic motion between the Tarim and NCB blocks, which were so far believed to have formed a rigid entity since at least the Jurassic. One interpretation could be that all Tarim and Qaidam Cretaceous paleomagnetic samples from red beds, but not those from Yumen and the NCB, suffered significant inclination shallowing, as observed in Cenozoic red beds from Central Asia. So far, we do not find support for this possibility. Possible tectonic interpretations include: (1) the existence of a large, as yet uncharted, tectonic discontinuity between Tarim and the NCB in the vicinity of the desert corridor near 95-100°E longitude; (2) the occurrence of significant deformation within southwestern Tarim, to the north of Yingjisha where paleomagnetic sites were obtained, or (3) persistent clockwise rotation of Tarim with respect to the NCB, for at least 20 Ma, at the rate found for current block kinematics.     

Strontium in rivers of the Baltic Basin

The rivers in the Baltic Basin drain a mixture of bedrocks ranging from Mesozoic-Paleozoic sediments in the south to Proterozoic-Archean intrusives in the north. The rivers in the sedimentary basin in the south have high concentrations of Sr, in the interval 100–500 µg l^–1 while the⁸⁷Sr/⁸⁶Sr ratio is close to that of seawater, i.e. 0.71. The northern rivers in the Precambrian shield area on the other hand have low Sr concentrations of 15–50 µg l^–1 with high⁸⁷Sr/⁸⁶Sr ratios of about to 0.73 (0.721–0.745). The riverine input of dissolved Sr to the brackish Baltic Sea approaches 60 tons year^–1, with a weighted mean concentration approaching 130 µg l^–1 and a weighted mean⁸⁷Sr/⁸⁶Sr ratio close to 0.712. Although the sedimentary area in the south supplies only about 43% of the total river discharge, it gives about 88% of the total Sr input. Because of this and the strong regional riverine variation in⁸⁷Sr/⁸⁶Sr ratio, Sr and its isotopes seem to be a convenient tool to unveil mixing relations of water masses in the northern Baltic Sea, provided high resolution analyses are applied. For an overall characterization of water mixing in the Baltic Sea, the Nd system will be superior to that of Sr.