Spatial and temporal variability of ground water recharge in central Australia: a tracer approach

Two environmental tracer methods are applied to the Ti-Tree Basin in central Australia to shed light on the importance of recharge from floodouts of ephemeral rivers in this arid environment. Ground water carbon-14 concentrations from boreholes are used to estimate the average recharge rate over the interval between where the ground water sample first entered the saturated zone and the bore. Environmental chloride concentrations in ground water samples provide estimates of the recharge rate at the exact point in the landscape where the sample entered the saturated zone. The results of the two tracer approaches indicate that recharge rates around one of the rivers and an extensive floodplain are generally higher than rates of diffuse recharge that occurs in areas of lower topographic relief. Ground water 2H/1H and 18O/16O compositions are all depleted in the heavier isotopes (delta2H = -67 per thousand to -50 per thousand; delta18O = -9.2 per thousand to -5.7%o) compared with the long-term, amount-weighted mean isotopic composition of rainfall in the area (delta2H = -33.8 per thousand; delta18O = -6.3 per thousand). This indicates that recharge throughout the basin occurs only after intense rainfall events of at least 150 to 200 mm/month. Finally, a recharge map is developed to highlight the spatial extent of the two recharge mechanisms. Floodout recharge to the freshest ground water (TDS <1,000 mg/L) is approximately 1.9 mm/year compared with a mean recharge rate of approximately 0.2 mm/year to the remainder of the basin. These findings have important implications for management of the ground water resource.     

Groundwater flow and solute transport at the Mourquong saline-water disposal basin, Murray Basin, southeastern Australia

Saline groundwater and drainage effluent from irrigation are commonly stored in some 200 natural and artificial saline-water disposal basins throughout the Murray-Darling Basin of Australia. Their impact on underlying aquifers and the River Murray, one of Australia's major water supplies, is of serious concern. In one such scheme, saline groundwater is pumped into Lake Mourquong, a natural groundwater discharge complex. The disposal basin is hydrodynamically restricted by low-permeability lacustrine clays, but there are vulnerable areas in the southeast where the clay is apparently missing. The extent of vertical and lateral leakage of basin brines and the processes controlling their migration are examined using (1) analyses of chloride and stable isotopes of water (²H/¹H and ¹⁸O/¹⁶O) to infer mixing between regional groundwater and lake water, and (2) the variable-density groundwater flow and solute-transport code SUTRA. Hydrochemical results indicate that evaporated disposal water has moved at least 100 m in an easterly direction and that there is negligible movement of brines in a southerly direction towards the River Murray. The model is used to consider various management scenarios. Salt-load movement to the River Murray was highest in a "worst-case" scenario with irrigation employed between the basin and the River Murray. Present-day operating conditions lead to little, if any, direct movement of brine from the basin into the river. Electronic Publication     

The termination of the last major phase of aeolian sand movement,coastal dunefields,Denmark

Optically stimulated luminescence (OSL) dating of sand samples from stabilized (or inactive) coastal dunes in Denmark provides information on the age of the termination phase of the last major aeolian activity period. A total of 26 sand samples were taken from four different coastal dunefields around the North Sea, Skagerrak and Kattegat coasts of Denmark. The OSL dates indicate that the last major phase of aeolian activity terminated between ad 1860 and 1905. Most of the dunes examined in this study were active around 1820, during a period documented to have been very stormy. A dune management scheme started around 1792, and this no doubt was a major cause of dunefield stabilization, but an overall decline of storminess, particularly spring and summer storminess, around the end of the 19th century must also have contributed to the increasing inactivity of coastal dunes. The new OSL dates on aeolian sand movement agree well with historical data and data from topographic maps on dune movement. This agreement supports the observation from earlier work that OSL dating of recent aeolian sand movement is accurate over the last few decades to centuries. Copyright © 2005 John Wiley & Sons, Ltd.     

The effects of rock uplift and rock resistance on river morphology in a subduction zone forearc,Oregon, USA

Relationships between riverbed morphology, concavity, rock type and rock uplift rate are examined to independently unravel the contribution of along-strike variations in lithology and rates of vertical deformation to the topographic relief of the Oregon coastal mountains. Lithologic control on river profile form is reflected by convexities and knickpoints in a number of longitudinal profiles and by general trends of concavity as a function of lithology. Volcanic and sedimentary rocks are the principal rock types underlying the northern Oregon Coast Ranges (between 46°30′ and 45°N) where mixed bedrock–alluvial channels dominate. Average concavity, θ, is 0·57 in this region. In the alluviated central Oregon Coast Ranges (between 45° and 44°N) values of concavity are, on average, the highest (θ = 0·82). South of 44°N, however, bedrock channels are common and θ = 0·73. Mixed bedrock–alluvial channels characterize rivers in the Klamath Mountains (from 43°N south; θ = 0·64). Rock uplift rates of ≥0·5 mm a⁻¹, mixed bedrock–alluvial channels, and concavities of 0·53–0·70 occur within the northernmost Coast Ranges and Klamath Mountains. For rivers flowing over volcanic rocks θ = 0·53, and θ = 0·72 for reaches crossing sedimentary rocks. Whereas channel type and concavity generally co-vary with lithology along much of the range, rivers between 44·5° and 43°N do not follow these trends. Concavities are generally greater than 0·70, alluvial channels are common, and river profiles lack knickpoints between 44·5° and 44°N, despite the fact that lithology is arguably invariant. Moreover, rock uplift rates in this region vary from low, ≤0·5 mm a⁻¹, to subsidence (<0 mm a⁻¹). These observations are consistent with models of transient river response to a decrease in uplift rate. Conversely, the rivers between 44° and 43°N have similar concavities and flow on the same mapped bedrock unit as the central region, but have bedrock channels and irregular longitudinal profiles, suggesting the river profiles reflect a transient response to an increase in uplift rate. If changes in rock uplift rate explain the differences in river profile form and morphology, it is unlikely that rock uplift and erosion are in steady state in the Oregon coastal mountains. Copyright © 2006 John Wiley & Sons, Ltd.     

Which way to the Moon?

Andrew Ball and Ian Crawford report on NASA's planning for a return to the Moon, with some thoughts on opportunities open to the UK.     

Remotely sensed estimation of forest canopy density: A comparison of the performance of four methods

In recent years, a number of alternative methods have been proposed to predict forest canopy density from remotely sensed data. To date, however, it remains difficult to decide which method to use, since their relative performance has never been evaluated. In this study the performance of: (1) an artificial neural network, (2) a multiple linear regression, (3) the forest canopy density mapper and (4) a maximum likelihood classification method was compared for prediction of forest canopy density using a Landsat ETM+ image. Comparison of confusion matrices revealed that the regression model performed significantly worse than the three other methods. These results were based on a z-test for comparison of weighted kappa statistics, which is an appropriate statistic for analysis of ranked categories. About 89% of the variance of the observed canopy density was explained by the artificial neural networks, which outperformed the other three methods in this respect. Moreover, the artificial neural networks gave an unbiased prediction, while other methods systematically under or over predicted forest canopy density. The choice of biased method could have a high impact on canopy density inventories.     

The relationships between tropical Atlantic sea level variability and major climate indices

Relationships were examined between variability in tropical Atlantic sea level and major climate indices with the use of TOPEX/POSEIDON altimeter and island tide gauge data with the aim of learning more about the external influences on the variability of the tropical Atlantic ocean. Possible important connections were found between indices related to the El Niño–Southern Oscillation (ENSO) and the sea levels in all three tropical regions (north, equatorial, and south), although the existence of only one major ENSO event within the decade of available altimetry means that a more complete investigation of the ENSO-dependence of Atlantic sea level changes has to await for the compilation of longer data sets. An additional link was found with the Indian Ocean Dipole (IOD) in the equatorial region, this perhaps surprising observation is probably an artifact of the similarity between IOD and ENSO time series in the 1990s. No evidence was obtained for significant correlations between tropical Atlantic sea level and North Atlantic Oscillation or Antarctic Oscillation Index. The most intriguing relationship observed was between the Quasi-Biennial Oscillation and sea level in a band centered approximately on 10°S. A plausible explanation for the relationship is lacking, but possibilities for further research are suggested.     

The dynamics of a channel-fed lava flow on Pico Partido volcano, Lanzarote

A short length of channel on Pico Partido volcano, Lanzarote, provides us the opportunity to examine the dynamics of lava flowing in a channel that extends over a sudden break in slope. The 1–2-m-wide, 0.5–2-m-deep channel was built during the 1730–1736 eruptions on Lanzarote and exhibits a sinuous, well-formed channel over a steep (11° slope) 100-m-long proximal section. Over-flow units comprising smooth pahoehoe sheet flow, as well as evidence on the inner channel walls for multiple (at least 11) flow levels, attest to unsteady flow in the channel. In addition, superelevation is apparent at each of the six bends along the proximal channel section. Superelevation results from banking of the lava as it moves around the bend thus causing preferential construction of the outer bank. As a result, the channel profile at each bend is asymmetric with an outer bank that is higher than the inner bank. Analysis of superelevation indicates flow velocities of ~8 m s^–1. Our analysis of the superelevation features is based on an inertia-gravity balance, which we show is appropriate, even though the down-channel flow is in laminar flow. We use a viscosity-gravity balance model, together with the velocities calculated from superelevation, to obtain viscosities in the range 25–60 Pa s (assuming that the lava behaved as a Newtonian liquid). Estimated volume fluxes are in the range 7–12 m³ s^–1. An apparent down-flow increase in derived volume flux may have resulted from variable supply or bulking up of the flow due to vesiculation. Where the channel moves over a sharp break in slope and onto slopes of ~6°, the channel becomes less well defined and widens considerably. At the break of slope, an elongate ridge extends across the channel. We speculate that this ridge was formed as a result of a reduction in velocity immediately below the break of slope to allow deposition of entrained material or accretion of lava to the channel bed as a result of a change in flow regime or depth.     

Synergy of multiple geophysical approaches to unravel explosive eruption conduit and source dynamics - A case study from Stromboli

An effective approach to understanding the dynamics of explosive volcanic eruptions and the conduit systems that drive them is through synergy of multiple data sets. Three data sets that lend themselves to ease of integration are seismic, infrasonic and thermal. Although approaches involving these data have been used to record volcanological phenomena since 1862, 1955 and 1965, respectively, their integrated use has only developed since 1999. When combined, these three data sets allow constraint of shallow system geometry and the dynamics of the explosive events that occur within that system. Using Stromboli volcano (Italy) as a case study, we review the complete range of geochemical and geophysical studies that can be applied. In doing so, we aim to show how integration of these diverse studies allows insights into a plumbing system and the dynamics of the eruptive activity that the system feeds. When combined at Stromboli, these data provide constraint of multiple system parameters including chamber depths, gas and magma fluxes, shallow system magma residence times, explosion source depths, and the rise/ejection velocities of ascending gas slugs and ejecta. In turn, these results allow various conduit and eruption dynamic models to be applied and tested.The persistent and repeated mildly explosive events that characterize Stromboli have been modeled in terms of the coalescence of gas within the magma to form large gas slugs that ascend the remaining portion of the conduit to burst at the free surface. Our integrated seismic, infrasonic and thermal data sets indicate that gas coalescence occurs at a depth of ∼260 m, with a typical event frequency of ∼9/h. Infrasonic and thermal data show the explosion source to be located 20-220 m below the vent. Thermal data give emission velocities for the ejected fragments of 8-20 m/s, which converts to gas jet velocities of 23-39 m/s. Tracking these parameters in space and time shows that, although eruptions at Stromboli can be grouped into two characteristic types (simple and complex-each of which characterizes a particular crater, NE and SW, respectively), events within each type show significant short-term variability. The system does, however, appear robust, maintaining its characteristic strombolian eruption style after significant effusive phases and more energetic explosive events.