Double-diffusive convection with a non-linear equation of state: Part I. The accurate conservation of properties in a two-layer system

If the equation of state is nonlinear, a given flux of heat across a double diffusive interface causes different buoyancy fluxes in the upper and lower layers. This results in different convective activity in the two layers and can lead to preferential entrainment across the interface in one direction (i.e. a migration of the interface). In this paper we derive the conservation equations for properties (e.g. heat and a solute) across a double diffusive interface between two well-mixed layers. A nondimensional measure of the entrainment across an interface and the most suitable choice for the buoyancy flux ratio are presented. Some surprising facts emerge. First, even for a linear equation of state and in the absence of direct entrainment across the interface, the flux of water across a finger interface is shown to be important. Second, for the heat-solute system, the heat balance equations for each well-mixed layer contain terms proportional to the heat of solution of the solute and the partial specific enthalpy of pure water in a seawater solution. Third, the rate of change of gravitational potential energy of the two-layer system is shown to have several extra terms in addition to the two commonly quoted major terms.     

Effects of detailed soil spatial information on watershed modeling across different model scales

Hydro-ecological modelers often use spatial variation of soil information derived from conventional soil surveys in simulation of hydro-ecological processes over watersheds at mesoscale (10–100 km²). Conventional soil surveys are not designed to provide the same level of spatial detail as terrain and vegetation inputs derived from digital terrain analysis and remote sensing techniques. Soil property layers derived from conventional soil surveys are often incompatible with detailed terrain and remotely sensed data due to their difference in scales. The objective of this research is to examine the effect of scale incompatibility between soil information and the detailed digital terrain data and remotely sensed information by comparing simulations of watershed processes based on the conventional soil map and those simulations based on detailed soil information across different simulation scales. The detailed soil spatial information was derived using a GIS (geographical information system), expert knowledge, and fuzzy logic based predictive mapping approach (Soil Land Inference Model, SoLIM). The Regional Hydro-Ecological Simulation System (RHESSys) is used to simulate two watershed processes: net photosynthesis and stream flow. The difference between simulation based on the conventional soil map and that based on the detailed predictive soil map at a given simulation scale is perceived to be the effect of scale incompatibility between conventional soil data and the rest of the (more detailed) data layers at that scale. Two modeling approaches were taken in this study: the lumped parameter approach and the distributed parameter approach. The results over two small watersheds indicate that the effect does not necessarily always increase or decrease as the simulation scale becomes finer or coarser. For a given watershed there seems to be a fixed scale at which the effect is consistently low for the simulated processes with both the lumped parameter approach and the distributed parameter approach.     

Knickpoint recession rate and catchment area: the case of uplifted rivers in Eastern Scotland

Knickpoint behaviour is a key to understanding both the landscape responses to a base‐level fall and the corresponding sediment fluxes from rejuvenated catchments, and must be accommodated in numerical models of large‐scale landscape evolution. Knickpoint recession in streams draining to glacio‐isostatically uplifted shorelines in eastern Scotland is used to assess whether knickpoint recession is a function of discharge (here represented by its surrogate, catchment area). Knickpoints are identified using DS plots (log slope versus log downstream distance). A statistically significant power relationship is found between distance of headward recession and catchment area. Such knickpoint recession data may be used to determine the values of m and n in the stream power law, E = KA^mSⁿ. The data have too many uncertainties, however, to judge definitively whether they are consistent with m = n = 1 (bedrock erosion is proportional to stream power and KPs should be maintained and propagate headwards) or m = 0·3, n = 0·7 (bedrock incision is proportional to shear stress and KPs do not propagate but degrade in place by rotation or replacement). Nonetheless, the E Scotland m and n values point to the dominance of catchment area (discharge) in determining knickpoint retreat rates and are therefore more consistent with the stream power law formulation in which bedrock erosion is proportional to stream power. Copyright © 2005 John Wiley & Sons, Ltd.     

Basal sediment evacuation by subglacial meltwater: suspended sediment transport from Haut Glacier d'Arolla,Switzerland

Proglacial suspended sediment transport was monitored at Haut Glacier d'Arolla, Switzerland, during the 1998 melt season to investigate the mechanisms of basal sediment evacuation by subglacial meltwater. Sub‐seasonal changes in relationships between suspended sediment transport and discharge demonstrate that the structure and hydraulics of the subglacial drainage system critically influenced how basal sediment was accessed and entrained. Under hydraulically inefficient subglacial drainage at the start of the melt season, sediment availability was generally high but sediment transport increased relatively slowly with discharge. Later in the melt season, sediment transport increased more rapidly with discharge as subglacial meltwater became confined to a spatially limited network of channels following removal of the seasonal snowpack from the ablation area. Flow capacity is inferred to have increased more rapidly with discharge within subglacial channels because rapid changes in discharge during highly peaked diurnal runoff cycles are likely to have been accommodated largely by changes in flow velocity. Basal sediment availability declined during channelization but increased throughout the remainder of the monitored period, resulting in very efficient basal sediment evacuation over the peak of the melt season. Increased basal sediment availability during the summer appears to have been linked to high diurnal water pressure variation within subglacial channels inferred from the strong increase in flow velocity with discharge. Basal sediment availability therefore appears likely to have been increased by (1) enhanced local ice‐bed separation leading to extra‐channel flow excursions and[sol ]or (2) the deformation of basal sediment towards low‐pressure channels due to a strong diurnally reversing hydraulic gradient between channels and areas of hydraulically less‐efficient drainage. Copyright © 2005 John Wiley & Sons, Ltd.     

Spatial Concentration in Australian Regional Development,Exogenous Shocks and Regional Demographic Outcomes: a South Australian case study

ABSTRACT

As a tribute to the massive contribution of our friend and colleague Graeme Hugo to the population and settlement geography of Australian rural areas, this paper presents a longitudinal study from his home State. It forms part of a wider study of the long-term demographic relationships between Australia’s rapidly growing regional cities and their surrounding functional regions. Of particular interest is the question of what effect the accelerating concentration of population and economic activity into a given regional city will have for the longer term demographic sustainability of its functional region as a whole. Taking the case of Port Lincoln, regional capital of most of South Australia’s Eyre Peninsula, it examines the nature of change in the functional region over the period 1947–2011, and investigates the forces feeding, and partly counteracting, the population concentration process, informed by concepts of evolutionary economic geography. In particular it traces the demographic impact (particularly differential migration and ageing trends) of exogenous shocks to the region’s essentially primary productive economic base during the period of major change from 1981 to 2011.     

Temporal evolution and compositional signatures of two supervolcanic systems recorded in zircons from Mangakino volcanic centre,New Zealand

Mangakino, the oldest rhyolitic caldera centre delineated in the Taupo Volcanic Zone of New Zealand, generated two very large (super-sized) ignimbrite eruptions, the 1.21 ± 0.04 Ma >500 km³ Ongatiti and ~1.0 Ma ~1,200 km³ Kidnappers events, the latter of which was followed after a short period of erosion by the ~200 km³ Rocky Hill eruption. We present U/Pb ages and trace-element analyses on zircons from pumice clasts from these three eruptions by Secondary Ion Mass Spectrometry (SIMS) using SHRIMP-RG instruments to illustrate the evolution of the respective magmatic systems. U–Pb age spectra from the Ongatiti imply growth of the magmatic system over ~250 kyr, with a peak of crystallisation around 1.32 Ma, ~100 kyr prior to eruption. The zircons are inferred to have then remained stable in a mush with little crystallisation and/or dissolution before later rejuvenation of the system at the lead-in to eruption. The paired Kidnappers and Rocky Hill eruptions have U–Pb zircon ages and geochemical signatures that suggest they were products of a common system grown over ~200 kyr. The Kidnappers and Rocky Hill samples show similar weakly bimodal age spectra, with peaks at 1.1 and 1.0 Ma, suggesting that an inherited antecrystic population was augmented by crystals grown at ages within uncertainty of the eruption age. In the Kidnappers, this younger age peak is dominantly seen in needle-shaped low U grains with aspect ratios of up to 18. In all three deposits, zircon cores show larger ranges and higher absolute concentrations of trace elements than zircon rims, consistent with zircon crystallisation from evolving melts undergoing crystal fractionation involving plagioclase and amphibole. Abundances and ratios of many trace elements frequently show variations between different sectors within single grains, even where there is no visible sector zoning in cathodoluminescence (CL) imaging. Substitution mechanisms, as reflected in the molar (Sc + Y + REE³⁺)/P ratio, differ in the same growth zone between the sides (along a-axis and b-axis: values approaching 1.0) and tips (c-axis: values between 1.5 and 5.0) of single crystals. These observations have implications for the use of zircons for tracking magmatic processes, particularly in techniques where CL zonation within crystals is not assessed and small analytical spot sizes cannot be achieved. These observations also limit applicability of the widely used Ti-in-zircon thermometer. The age spectra for the Ongatiti and Kidnappers/Rocky Hill samples indicate that both magmatic systems were newly built in the time-breaks after respective previous large eruptions from Mangakino. Trace element variations defining three-component mixing suggest that zircons, sourced from multiple melts, contributed to the population in each system.     

Trace Element Partitioning and Accessory Phase Saturation during H2O-Saturated Melting of Basalt with Implications for Subduction Zone Chemical Fluxes

Experimental phase equilibrium and trace element partitioningdata are reported for H₂O-saturated mid-ocean ridge basalt at2·5 GPa, 750–900°C and oxygen fugacities atthe nickel–nickel oxide buffer. Garnet, omphacite andrutile are present at all temperatures. Amphibole and epidotedisappear as residual phases above 800°C; allanite appearsabove 750°C. The Na–Al-rich silicate glass presentin all run products is likely to have quenched from a supercriticalliquid. Trace element analyses of glasses demonstrate the importantcontrol exerted by residual minerals on liquid chemistry. Inaddition to garnet, which controls heavy rare earth elements(HREE) and Sc, and rutile, which controls Ti, Nb and Ta, allanitebuffers the light REE (LREE; La–Sm) contents of liquidsto relatively low levels and preferentially holds back Th relativeto U. In agreement with previous experimental and metamorphicstudies we propose that residual allanite plays a key role inselectively retaining trace elements in the slab during subduction.Experimental data and analyses of allanite-bearing volcanicrocks are used to derive a model for allanite solubility inliquids as a function of pressure, temperature, anhydrous liquidcomposition and LREE content. The large temperature dependenceof allanite solubility is very similar to that previously determinedfor monazite. Our model, fitted to 48 datapoints, retrievesLREE solubility (in ppm) to within a factor of 1· 40over a pressure range of 0–4 GPa, temperature range of700–1200°C and for liquids with anhydrous SiO₂ contentsof 50–84 wt %. This uncertainty in LREE content is equivalentto a temperature uncertainty of only ± 27°C at 1000K, indicating the potential of allanite as a geothermometer.Silicic liquids from either basaltic or sedimentary protolithswill be saturated in allanite except for Ca-poor protolithsor at very high temperatures. For conventional subduction geothermsthe low solubility of LREE (+ Th) in liquids raises questionsabout the mechanism of LREE + Th transport from slab to wedge.It is suggested either that, locally, temperatures experiencedby the slab are high enough to eliminate allanite in the residueor that substantial volumes of H₂O-rich fluids must pass throughthe mantle wedge prior to melting. The solubility of accessoryphases in fluids derived from subducted rocks can provide importantconstraints on subduction zone thermal structure. KEY WORDS: subduction; experimental petrology; allanite; solubility; supercritical liquid; eclogite     

Application of a multimetric fish index to assess the environmental condition of South African estuaries

A multimetric fish index, the Estuarine Fish Community Index (EFCI) of Harrison and Whitfield (2004), was applied to data collected for 190 South African estuaries. Estuaries spanned three biogeographic regions and included three distinct estuarine typologies. The EFCI is based on 14 metrics or measures that represent four broad fish community attributes: species diversity and composition, species abundance, nursery function, and trophic integrity. Metric reference conditions and scoring criteria were developed separately for each estuary type within each zoogeographic region. The final EFCI was applied to each estuary by comparing its fish community with the appropriate reference. Index values ranged between 18 (very poor) and 66 (very good). A comparison of the EFCI with independent measures of estuarine condition revealed that the index was able to effectively differentiate between poor and good quality sites. Applying the EFCI to estuaries in which multiple samples were taken also showed that the index is reproducible. The EFCI is both a robust and sensitive method for assessing the ecological condition of estuarine systems; it is also an effective communication tool for converting ecological information into an easily understood format for managers, policy makers, and the general public.