Artificial pumping errors in the Kool-Parker scaling model of soil moisture hysteresis

Hysteresis models that eliminate the artificial pumping errors associated with the Kool-Parker (KP) soil moisture hysteresis model, such as the Parker-Lenhard (PL) method, can be computationally demanding in unsaturated transport models since they need to retain the wetting-drying history of the system. The pumping errors in these models need to be eliminated for correct simulation of cyclical systems (e.g. transport above a tidally forced watertable, infiltration and redistribution under periodic irrigation) if the soils exhibit significant hysteresis. A modification is made here to the PL method that allows it to be more readily applied to numerical models by eliminating the need to store a large number of soil moisture reversal points. The modified-PL method largely eliminates any artificial pumping error and so essentially retains the accuracy of the original PL approach. The modified-PL method is implemented in HYDRUS-1D (version 2.0), which is then used to simulate cyclic capillary fringe dynamics to show the influence of removing artificial pumping errors and to demonstrate the ease of implementation. Artificial pumping errors are shown to be significant for the soils and system characteristics used here in numerical experiments of transport above a fluctuating watertable. (c) 2005 Elsevier B.V. All rights reserved.

Efficient accommodation of local minima in watershed model calibration

The Gauss-Marquardt-Levenberg (GML) method of computer-based parameter estimation, in common with other gradient-based approaches, suffers from the drawback that it may become trapped in local objective function minima, and thus report optimized parameter values that are not, in fact, optimized at all. This can seriously degrade its utility in the calibration of watershed models where local optima abound. Nevertheless, the method also has advantages, chief among these being its model-run efficiency, and its ability to report useful information on parameter sensitivities and covariances as a by-product of its use. It is also easily adapted to maintain this efficiency in the face of potential numerical problems (that adversely affect all parameter estimation methodologies) caused by parameter insensitivity and/or parameter correlation. The present paper presents two algorithmic enhancements to the GML method that retain its strengths, but which overcome its weaknesses in the face of local optima. Using the first of these methods an intelligent search for better parameter sets is conducted in parameter subspaces of decreasing dimensionality when progress of the parameter estimation process is slowed either by numerical instability incurred through problem ill-posedness, or when a local objective function minimum is encountered. The second methodology minimizes the chance of successive GML parameter estimation runs finding the same objective function minimum by starting successive runs at points that are maximally removed from previous parameter trajectories. As well as enhancing the ability of a GML-based method to find the global objective function minimum, the latter technique can also be used to find the locations of many non-global optima (should they exist) in parameter space. This can provide a useful means of inquiring into the well-posedness of a parameter estimation problem, and for detecting the presence of bimodal parameter and predictive probability distributions. The new methodologies are demonstrated by calibrating a Hydrological Simulation Program-FORTRAN (HSPF) model against a time series of daily flows. Comparison with the SCE-UA method in this calibration context demonstrates a high level of comparative model run efficiency for the new method. (c) 2006 Elsevier B.V. All rights reserved.

A high-precision record of mid-late Holocene sea-level events from emergent coral pavements in the Houtman Abrolhos Islands, southwest Australia

Early work on sea-levels in southwest Australia claimed to recognise a Holocene sea-level highstand which was not seen in better known sea-level records elsewhere at the time. More recent work has confirmed that a mid-Holocene highstand Occurred about 6 kyr ago. As new data on oscillating sea-levels from the region have recently been published, a high continuity, precisely dated and accurately surveyed record was obtained from emergent coral pavements in the leeward Houtman Abrolhos Islands (Serventy Island), a tectonically stable region from where good-quality Holocene sea-level data have been previously obtained from corals. From the mid-Holocene highstand ca. 7 U/Th kyr ago, sea-level declined linearly during the remainder of the Holocene as the carbonate platform prograded leewards. Hydro-isostatic controls are probably significant in the record. (c) 2005 Elsevier Ltd and INQUA. All rights reserved.

Characterisation of Chinese Tang sancai from Gongxian and Yaozhou kilns using ICP-MS trace element and TIMS Sr–Nd isotopic analysis

Tang sancai is one of the most important types of Chinese ceramics. To determine the provenance of Tang sancai is important to study ancient trade and other issues. In this paper we compare ICP-MS trace elements and TIMS Sr-Nd isotopes of visibly similar Tang sancai from two major production centres Gongxian and Yaozhou. The variation in contents/ratios of many of > 40 trace elements is small for samples from Gongxian, yet is considerably bigger for that from Yaozhou. However, the variation in Sr-87/Sr-86 and Nd-143/Nd-144 isotopic ratios is very small for samples from both places. Gongxian and Yaozhou samples have distinctive Sr-Nd isotopic and trace element features despite their similarity in major elements, and these analysis data can be interpreted with geochemistry, indicating that Sr and Nd isotopes have great potential in ceramic provenance studies. The distinct characterisation of these samples provides valuable criteria for identifying provenance of Tang sancai of uncertain origin. Two modern fakes are also analysed, and they can as well be distinguished from antique Tang sancai using above criteria. (c) 2005 Published by Elsevier Ltd.

Integration of ice-core, marine and terrestrial records for the Australian Last Glacial Maximum and Termination: a contribution from the OZ INTIMATE group

The degree to which Southern Hemisphere climatic changes during the end of the last glacial period and early Holocene (30-8 ka) were influenced or initiated by events occurring in the high latitudes of the Northern Hemisphere is a complex issue. There is conflicting evidence for the degree of hemispheric 'teleconnection' and an unresolved debate as to the principle forcing mechanism(s). The available hypotheses are difficult to test robustly, however, because the few detailed palaeoclimatic records in the Southern Hemisphere are widely dispersed and lack duplication. Here we present climatic and environmental reconstructions from across Australia, a key region of the Southern Hemisphere because of the range of environments it covers and the potentially important role regional atmospheric and oceanic controls play in global climate change. We identify a general scheme of events for the end of the last glacial period and early Holocene but a detailed reconstruction proved problematic. Significant progress in climate quantification and geochronological control is now urgently required to robustly investigate change through this period. Copyright (c) 2006 John Wiley & Sons, Ltd.

High-precision TIMS U-series and AMS C-14 dating of a coral reef lagoon sediment core from southern South China Sea

Accurate dating of lagoon sediments has been a difficult problem, although lagoon profiles, usually with high deposition rates, have a great potential for high-resolution climate reconstruction. We report 26 high-precision TIMS U-series dates (on 25 coral branches) and five AMS C-14 dates (on foraminifera) for a 15.4-m long lagoon core from Yongshu Reef, Nansha area, southern South China Sea. All the dates are in the correct stratigraphical sequence, providing the best chronology so far reported for lagoon deposits. The results reveal a similar to 4000-a continuous depositional history, with sedimentation rates varying from 0.8 to 24.6 mm a(-1), with an average of 3.85 mm a(-1), which corresponds to an average net carbonate accumulation rate of similar to 2700 g CaCO3 m(-2) a(-1), significantly higher than the mean value (800 +/- 400 g CaCO3 m(-2) a(-1)) used for lagoons in general in previous studies of global carbonate budget. Episodes of accelerated depositions within the last 1000 years correlate well with strong storm events identified by U-series dates of storm-transported coral blocks in the area. However, in the longer term, the sedimentation rates during the past 1000 years were much higher than earlier on, probably due to more vigorous wave-reef interaction as a result of relative sea-level fall since 500 AD and expansion of reef flat area, supplying more sediments. The coral TIMS U-series ages and foraminifera AMS 14C dates reveal intriguing apparent radiocarbon reservoir ages (R) from 572 to 1052 years, which are much higher than global mean values of similar to 400 years. (c) 2006 Elsevier Ltd. All rights reserved.

Climatic variability in the southwest Pacific during the Last Termination (20-10 kyr BP)

The degree to which palaeoclimatic changes in the Southern Hemisphere co-varied with events in the high latitude Northern Hemisphere during the Last Termination is a contentious issue, with conflicting evidence for the degree of 'teleconnection' between different regions of the Southern Hemisphere. The available hypotheses are difficult to test robustly, however, because there are few detailed palaeoclimatic records in the Southern Hemisphere. Here we present climatic reconstructions from the southwestern Pacific, a key region in the Southern Hemisphere because of the potentially important role it plays in global climate change. The reconstructions for the period 20-10 kyr BP were obtained from five sites along a transect from southern New Zealand, through Australia to Indonesia, supported by 125 calibrated C-14 ages. Two periods of significant climatic change can be identified across the region at around 17 and 14.2 cal kyr BP, most probably associated with the onset of warming in the West Pacific Warm Pool and the collapse of Antarctic ice during Meltwater Pulse-1A, respectively. The severe geochronological constraints that inherently afflict age models based on radiocarbon dating and the lack of quantified climatic parameters make more detailed interpretations problematic, however. There is an urgent need to address the geochronological limitations, and to develop more precise and quantified estimates of the pronounced climate variations that clearly affected this region during the Last Termination. (c) 2005 Elsevier Ltd. All rights reserved.

Geochemical changes recorded in Lynch's Crater, Northeastern Australia, over the past 50 ka

There are many geochemical reconstructions of environmental change in the mid and high latitudes but relatively few in the tropical latitudes, despite their considerable potential for reconstructing environmental processes that cannot be identified using more traditional proxies. Here we present one reconstruction of environmental change for the tropics. This reconstruction covers the past 50 ka using a suite of geochemical data from the high-resolution sequence of Lynch's Crater in northeast Queensland, Australia, a region highly sensitive to El Nino-Southern Oscillation (ENSO) activity. The 23 major oxides and trace elements measured Could be summarised by extracting three axes using principal components analysis (accounting for 72% of the variability). The data indicate that the greatest variability in the geochemical data accounted for erosional activity within the catchment that was associated with past changes in the frequency of ENSO activity (though this was less sensitive during wetter periods, probably as a result of buffering by high vegetation cover). The remaining variability was largely explained by elements that form complexes with organic compounds (e.g., humic acids) and those that are important nutrients for specific vegetation types (and therefore a measure of vegetation distribution). For more detailed reconstructions, further work is required to disentangle the complex controls of clements within sedimentary sequences. (c) 2005 Elsevier B.V. All rights reserved.

Hydrogen from coal: Production and utilisation technologies

Error condition detected Although coal may be viewed as a dirty fuel due to its high greenhouse emissions when combusted, a strong case can be made for coal to be a major world source of clean H-2 energy. Apart from the fact that resources of coal will outlast oil and natural gas by centuries, there is a shift towards developing environmentally benign coal technologies, which can lead to high energy conversion efficiencies and low air pollution emissions as compared to conventional coal fired power generation plant. There are currently several world research and industrial development projects in the areas of Integrated Gasification Combined Cycles (IGCC) and Integrated Gasification Fuel Cell (IGFC) systems. In such systems, there is a need to integrate complex unit operations including gasifiers, gas separation and cleaning units, water gas shift reactors, turbines, heat exchangers, steam generators and fuel cells. IGFC systems tested in the USA, Europe and Japan employing gasifiers (Texaco, Lurgi and Eagle) and fuel cells have resulted in energy conversions at efficiency of 47.5% (HHV) which is much higher than the 30-35% efficiency of conventional coal fired power generation. Solid oxide fuel cells (SOFC) and molten carbonate fuel cells (MCFC) are the front runners in energy production from coal gases. These fuel cells can operate at high temperatures and are robust to gas poisoning impurities. IGCC and IGFC technologies are expensive and currently economically uncompetitive as compared to established and mature power generation technology. However, further efficiency and technology improvements coupled with world pressures on limitation of greenhouse gases and other gaseous pollutants could make IGCC/IGFC technically and economically viable for hydrogen production and utilisation in clean and environmentally benign energy systems. (c) 2005 Elsevier B.V. All rights reserved.

AUSLEM (AUStralian Land Erodibility Model): A tool for identifying wind erosion hazard in Australia

We present AUSLEM (AUStralian Land Erodibility Model), a land erodibility modelling system that utilizes a rule-set of surficial and climatic thresholds applied through a Geographic Information System (GIs) modelling framework to predict landscape susceptibility to wind erosion. AUSLEM is distinctive in that it quantitatively assesses landscape susceptibility to wind erosion at a 5 x 5 km. spatial resolution on a monthly time-step across Australia. The system was implemented for representative wet (1984), dry (1994), and average rainfall (1997) years with corresponding low, high and moderate dust storm day frequencies. Results demonstrate that AUSLEM can identify landscape erodibility, and provide an interpretation of the physical nature and distribution of erodible landscapes in Australia. Further, results offer an assessment of the dynamic tendencies of erodibility in space and time in response to the El Nino Southern Oscillation (ENSO) and seasonal synoptic scale climate variability. A comparative analysis of AUSLEM output with independent national and international wind erosion, atmospheric aerosol and dust event records indicates a high level of model competency. (c) 2006 Elsevier B.V. All rights reserved.

Tidal impacts on riparian salinities near estuaries

Groundwater-dependent riparian biota is known to be sensitive to changes in soil and groundwater salinity in estuarine systems. The groundwater flow and salinity behaviour in a phreatic aquifer adjoining a partially penetrating, tidal. estuary is investigated through two-dimensional numerical experiments for a lateral cross-section, which explore the influence of factors, such as aquifer and soil materials, tidal amptitudes, and regional groundwater hydraulic gradients. The density contrast between estuarine water and the fresh groundwater drives saltwater penetration of the aquifer even in the case of a marked groundwater hydraulic gradient towards the estuary. We show that tidal fluctuations in estuaries can significantly affect the groundwater salinity distribution in adjacent density-stratified phreatic aquifers. This has consequences for the expected distribution of salinity-sensitive biota in the hyporheic zone as well as vegetation and fauna dependent on water in the riparian soil and aquifer. The shape of the dense saltwater wedge propagating into the adjacent groundwater system is also modified by the estuarine tidal signal, although this effect appears to have only minor influence on the maximum distance penetrated into the aquifer (i.e., location of the 'toe' of the wedge). Tide-induced changes to riparian groundwater salinity are advection-driven, as evidenced by the modified time-averaged groundwater flow dynamics. (c) 2006 Elsevier B.V. All rights reserved.

Partitioning of deformation within a subduction channel during exhumation of high-pressure rocks: a case study from the Western Alps

The metamorphic belt of the Western Alps was subjected to widespread extensional tectonism at the end of the Eocene (ca. 45-35 Ma). Extension was accommodated by hinterland-directed movements along gently inclined extensional shear zones, which facilitated rapid exhumation of high-pressure and ultra-high-pressure rocks. This deformation resulted in a normal metamorphic sequence. Extension in the inner parts of the Western Alps was coeval with shortening at the front of the belt (foreland-directed thrusts), which took place during decompression, and emplaced higher grade metamorphic units over lower grade metamorphic rocks, thus forming an inverse metamorphic sequence. Two mechanisms for this extensional episode are discussed: (1) collapse of an overthickened lithosphere, and (2) internal readjustments within the orogenic wedge due to subduction channel dynamics. We favour the latter mechanism because it can account for the development of the observed inverse and normal metamorphic sequences along foreland-directed thrusts and hinterland-directed detachments, respectively. This hypothesis is supported by published structural, metamorphic and geochronological data from four geological transects through the Western Alps. This study also emphasizes the importance of post-shearing deformation (e.g. horizontal buckling versus vertical flattening), which can modify the distribution of hinterland- and foreland-directed shear zones in orogenic belts. (c) 2006 Elsevier Ltd. All rights reserved.

Mapping structural parameters and species composition of riparian vegetation using IKONOS and Landsat ETM+ data in Australian tropical savannahs

Government agencies responsible for riparian environments are assessing the utility of remote sensing for mapping and monitoring vegetation structural parameters. The objective of this work was to evaluate Ikonos and Landsat-7 ETM+ imagery for mapping structural parameters and species composition of riparian vegetation in Australian tropical savannahs for a section of Keelbottom Creek, Queensland, Australia. Vegetation indices and image texture from Ikonos data were used for estimating leaf area index (R-2 = 0.13) and canopy percentage foliage cover (R-2 = 0.86). Pan-sharpened Ikonos data were used to map riparian species composition (overall accuracy = 55 percent) and riparian zone width (accuracy within +/- 3 m). Tree crowns could not be automatically delineated due to the lack of contrast between canopies and adjacent grass cover. The ETM+ imagery was suited for mapping the extent of riparian zones. Results presented demonstrate the capabilities of high and moderate spatial resolution imagery for mapping properties of riparian zones.

A worldwide correlation for exponential bed particle size variation in subaerial aqueous flows

The particle size of the bed sediments in or on many natural streams, alluvial fans, laboratory flumes, irrigation canals and mine waste deltas varies exponentially with distance along the stream. A plot of the available worldwide exponential bed particle size diminution coefficient data against stream length is presented which shows that all the data lie within a single narrow band extending over virtually the whole range of stream lengths and bed sediment particle sizes found on Earth. This correlation applies to both natural and artificial flows with both sand and gravel beds, irrespective of either the solids concentration or whether normal or reverse sorting occurs. This strongly suggests that there are common mechanisms underlying the exponential diminution of bed particles in subaerial aqueous flows of all kinds. Thus existing models of sorting and abrasion applicable to some such flows may be applicable to others. A comparison of exponential laboratory abrasion and field diminution coefficients suggests that abrasion is unlikely to be significant in gravel and sand bed streams shorter than about 10 km to 100 km, and about 500 km, respectively. Copyright (C) 1999 John Wiley & Sons, Ltd.

A numerical study of steady state evaporative conditions applied to mine tailings

The suction profile of a desiccating soil is dependent on the water table depth, the soil-water retention characteristics, and the climatic conditions. In this paper, an unsaturated flow model, which simulates both liquid and vapour flow, was used to investigate the effects of varying the water table depth and the evaporation rate on the evaporative fluxes from a desiccating tailings deposit under steady-state conditions. Results obtained showed that at a critical evaporation rate, beyond which evaporation is no longer dictated by climatic conditions, the matric suction profiles remain basically unchanged. The critical evaporation rate varies inversely with the water table depth. It is associated with the maximum evaporative flux that might be extracted from a soil at steady-state conditions. The time required to establish steady-state conditions is directly proportional to the water table depth, and it acquires a maximum value at the critical evaporation rate. A detailed investigation of the movement of the drying front demonstrated the significance of attaining a matric suction of about 3000 kPa on the contribution to flow in the vapour phase.