Truncation errors in finite difference models for solute transport equation with first-order reaction

The truncation errors associated with finite difference solutions of the advection-dispersion equation with first-order reaction are formulated from a Taylor analysis. The error expressions are based on a general form of the corresponding difference equation and a temporally and spatially weighted parametric approach is used for differentiating among the various finite difference schemes. The numerical truncation errors are defined using Peclet and Courant numbers and a new Sink/Source dimensionless number. It is shown that all of the finite difference schemes suffer from truncation errors. Tn particular it is shown that the Crank-Nicolson approximation scheme does not have second order accuracy for this case. The effects of these truncation errors on the solution of an advection-dispersion equation with a first order reaction term are demonstrated by comparison with an analytical solution. The results show that these errors are not negligible and that correcting the finite difference scheme for them results in a more accurate solution. (C) 1999 Elsevier Science B.V. All rights reserved.

El Nino and arboviral disease prediction

Recent El Nino events have stimulated interest in the development of modeling techniques to forecast extremes of climate and related health events. Previous studies have documented associations between specific climate variables (particularly temperature and rainfall) and outbreaks of arboviral disease. In some countries, such diseases are sensitive to Fl Nino. Here we describe a climate-based model for the prediction of Ross River virus epidemics in Australia. From a literature search and data on case notifications, we determined in which years there were epidemics of Ross River virus in southern Australia between 1928 and 1998. Predictor variables were monthly Southern Oscillation index values for the year of an epidemic or lagged by 1 year. We found that in southeastern states, epidemic years were well predicted by monthly Southern Oscillation index values in January and September in the previous year. The model forecasts that there is a high probability of epidemic Ross River virus in the southern states of Australia in 1999. We conclude that epidemics of arboviral disease can, at least in principle, be predicted on the basis of climate relationships.

CXTANNEAL: an improved program for estimating solute transport parameters

CXTANNEAL is a program for analysing contaminant transport in soils. The code, written in Fortran 77, is a modified version of CXTFIT, a commonly used package for estimating solute transport parameters in soils. The improvement of the present code is that it includes simulated annealing as the optimization technique for curve fitting. Tests with hypothetical data show that CXTANNEAL performs better than the original code in searching for optimal parameter estimates. To reduce the computational time, a parallel version of CXTANNEAL (CXTANNEAL_P) was also developed. (C) 1999 Elsevier Science Ltd. All rights reserved.

Submarine groundwater discharge and associated chemical input to a coastal sea

This paper presents a theoretical model of flow and chemical transport processes in subterranean estuaries (unconfined brackish groundwater aquifers at the ocean-land interface). The model shows that groundwater circulation and oscillating flow, caused by wave setup and tide, may constitute up to 96% of submarine groundwater discharge (SGWD) compared with 4% due to the net groundwater discharge. While these local flow processes do not change the total amount of land-derived chemical input to the ocean over a long period (e.g., yearly), they induce fluctuations of the chemical transfer rate as the aquifer undergoes saltwater intrusion. This may result in a substantial increase in chemical fluxes to the ocean over a short period (e.g., monthly and by a factor of 20 above the averaged level), imposing a possible threat to the marine environment. These results are essentially consistent with the experimental findings of Moore [1996] and have important implications for coastal resources management.

PAM chlorophyll fluorometry: a new in situ technique for stress assessment in scleractinian corals, used to examine the effects of cyanide from cyanide fishing

Sodium cyanide is being used on reefs in the Asia-Pacific region to capture live fish for the aquarium industry, and to supply a rapidly growing, restaurant-based demand, The effects of cyanide on reef biota have not been fully explored. To investigate its effect on hard corals, we exposed small branch lips of Stylophora pistillata and Acropora aspera to cyanide concentrations estimated to occur during cyanide fishing. Pulse amplitude modulation (PAM) chlorophyll fluorescence techniques were used to examine photoinhibition and photosynthetic electron transport in the symbiotic algae (zooxanthellae) in the tissues of the corals, These measurements were made in situ and in real time using a recently developed submersible PAM fluorometer. In S. pistillata. exposure to cyanide resulted in an almost complete cessation in photosynthetic electron transport rate. Both species displayed marked decreases in the ratio of variable fluorescence (F-v) to maximal fluorescence (F-m) (dark-adapted F-v/F-m), following exposure to cyanide, signifying a decrease in photochemical efficiency. Dark-adapted F-v/F-m recovered to normal levels in similar to 6 d, although intense tissue discolouration, a phenomenon well-recognised as coral 'bleaching' was observed during this period, Bleaching was caused by loss of zooxanthellae from the coral tissues, a well-recognised sub-lethal stress response of corals. Using the technique of chlorophyll fluorescence quenching analysis, corals exposed to cyanide did not show light activation of Calvin cycle enzymes and developed high levels of non-photochemical quenching (q(N)), signifying the photoprotective dissipation of excess light as heat, These features are symptomatic of the known properties of cyanide as an inhibitor of enzymes of the Calvin cycle. The results of this in situ study show that an impairment of zooxanthellar photosynthesis is; the site of cyanide-mediated toxicity, and is the cue that causes corals to release their symbiotic zooxanthellac following cyanide exposure. This study demonstrates the efficacy of PBM fluorometry as a new tool for in situ stress assessment in zooxanthellate scleractinian corals. (C) 1999 Elsevier Science Ltd. All rights reserved.

Groundwater waves in a coastal aquifer: A new governing equation including vertical effects and capillarity

Groundwater waves, that is, water table fluctuations, are a natural phenomenon in coastal aquifers. They represent an important part of the interaction between the ocean and aquifer and affect the mass exchange between them. This paper presents a new groundwater wave equation. Because it includes the effects of vertical flows and capillarity, the new equation is applicable to both intermediate-depth aquifers and high-frequency waves. Compared with the wave equation derived by Nielsen ed al. [1997], the present equation provides a closer representation of groundwater waves. In particular, it predicts high-frequency water table fluctuations as observed in the field. A validation of the new equation has been carried out by comparing the analytical solutions to it with predictions from direct simulations using the numerical model SUTRA. The effects of various physical parameters and their relative importance are also discussed.

Determining forest structural attributes using an inverted geometric-optical model in mixed eucalypt forests, Southeast Queensland, Australia

The Montreal Process indicators are intended to provide a common framework for assessing and reviewing progress toward sustainable forest management. The potential of a combined geometrical-optical/spectral mixture analysis model was assessed for mapping the Montreal Process age class and successional age indicators at a regional scale using Landsat Thematic data. The project location is an area of eucalyptus forest in Emu Creek State Forest, Southeast Queensland, Australia. A quantitative model relating the spectral reflectance of a forest to the illumination geometry, slope, and aspect of the terrain surface and the size, shape, and density, and canopy size. Inversion of this model necessitated the use of spectral mixture analysis to recover subpixel information on the fractional extent of ground scene elements (such as sunlit canopy, shaded canopy, sunlit background, and shaded background). Results obtained fron a sensitivity analysis allowed improved allocation of resources to maximize the predictive accuracy of the model. It was found that modeled estimates of crown cover projection, canopy size, and tree densities had significant agreement with field and air photo-interpreted estimates. However, the accuracy of the successional stage classification was limited. The results obtained highlight the potential for future integration of high and moderate spatial resolution-imaging sensors for monitoring forest structure and condition. (C) Elsevier Science Inc., 2000.

Increasing the efficiency of solute leaching: impacts of flow interruption with drainage of the "preferential flow paths"

Most soils contain preferential flow paths that can impact on solute mobility. Solutes can move rapidly down the preferential flow paths with high pore-water velocities, but can be held in the less permeable region of the soil matrix with low pore-water velocities, thereby reducing the efficiency of leaching. In this study, we conducted leaching experiments with interruption of the flow and drainage of the main flow paths to assess the efficiency of this type of leaching. We compared our experimental results to a simple analytical model, which predicts the influence of the variations in concentration gradients within a single spherical aggregate (SSA) surrounded by preferential flow paths on leaching. We used large (length: 300 mm, diameter: 216 mm) undisturbed field soil cores from two contrasting soil types. To carry out intermittent leaching experiments, the field soil cores were first saturated with tracer solution (CaBr2), and background solution (CaCl2) was applied to mimic a leaching event. The cores were then drained at 25- to 30-cm suction to empty the main flow paths to mimic a dry period during which solutes could redistribute within the undrained region. We also conducted continuous leaching experiments to assess the impact of the dry periods on the efficiency of leaching. The flow interruptions with drainage enhanced leaching by 10-20% for our soils, which was consistent with the model's prediction, given an optimised equivalent aggregate radius for each soil. This parameter quantifies the time scales that characterise diffusion within the undrained region of the soil, and allows us to calculate the duration of the leaching events and interruption periods that would lead to more efficient leaching. Application of these methodologies will aid development of strategies for improving management of chemicals in soils, needed in managing salts in soils, in improving fertiliser efficiency, and in reclaiming contaminated soils. (C) 2000 Elsevier Science B.V. All rights reserved.

Overview, critical assessment, and conservation implications of koala distribution and abundance

Regional and national surveys provide a broadscale description of the koala's present distribution in Australia. A detailed understanding of its distribution is precluded, however, by past and continuing land clearing across large parts of the koala's range. Koala population density increased in some regions during the late 1800s and then declined dramatically in the early 1900s. The decline was associated with habitat loss, hunting, disease, fire, and drought. Declines are continuing in Queensland and New South Wales. In contrast, dense koala populations in habitat isolates in Victoria and South Australia are managed to reduce population size and browse damage. Current understanding of koala distribution and abundance suggests that the species does not meet Australian criteria as endangered or vulnerable fauna. Its conservation status needs to be reviewed, however, in light of the extensive land clearing in New South Wales and Queensland since the last (1980s) broadscale surveys. Consequently, we recommend that broadacre clearing by curtailed in New South Wales and Queensland and that regular, comprehensive, standardized, national koala surveys be undertaken. Given the fragmentation of koala habitat and regional differences in the status of the koala, we recommended that studies on regional variation in the koala be intensified and that koala ecology in fragmented and naturally restricted habitats be developed. More generally, the National Koala Conservation Strategy should be implemented.

Cultural revolution in whale songs

Humpbacks have picked up a catchy tune sung by immigrants from a distant ocean. The song patterns of humpback whales (Megaptera novaeangliae) depend on where they live, with populations inhabiting different ocean basins normally singing quite distinct songs. Here we record a unique and radical song change in the song of humpback whales in the Pacific Ocean off the Australian east coast. Their song was replaced rapidly and completely by the song of the Australian west coast population from the Indian Ocean, apparently as a result of the introduction of only a small number of 'foreign' singers. Such a revolutionary change is unprecedented in animal cultural vocal traditions and suggests that novelty may stimulate change in humpback whale songs.

Tidal fluctuations in a leaky confined aquifer: Dynamic effects of an overlying phreatic aquifer

Tidal fluctuations in a leaky confined coastal aquifer are damped significantly due to leakage into an overlying phreatic aquifer. Jiao and Tang [1999] presented an analytical solution to a simple model describing this phenomenon. Their solution assumes that the tidal fluctuations in the overlying phreatic aquifer are negligible (i.e,, a static phreatic aquifer), Here we examine dynamic effects of the overlying aquifer based on a new approximate analytical solution. The numerical results indicate that the dynamic effects can be significant for a relatively large leakage and a high transmissivity of the phreatic aquifer.

A method for setting the size of plant conservation target areas

Realistic time frames in which management decisions are made often preclude the completion of the detailed analyses necessary for conservation planning. Under these circumstances, efficient alternatives may assist in approximating the results of more thorough studies that require extensive resources and time. We outline a set of concepts and formulas that may be used in lieu of detailed population viability analyses and habitat modeling exercises to estimate the protected areas required to provide desirable conservation outcomes for a suite of threatened plant species. We used expert judgment of parameters and assessment of a population size that results in a specified quasiextinction risk based on simple dynamic models The area required to support a population of this size is adjusted to take into account deterministic and stochastic human influences, including small-scale disturbance deterministic trends such as habitat loss, and changes in population density through processes such as predation and competition. We set targets for different disturbance regimes and geographic regions. We applied our methods to Banksia cuneata, Boronia keysii, and Parsonsia dorrigoensis, resulting in target areas for conservation of 1102, 733, and 1084 ha, respectively. These results provide guidance on target areas and priorities for conservation strategies.

Sudden drawdown and drainage of a horizontal aquifer

Drainage of a saturated horizontal aquifer following a sudden drawdown is reanalyzed using the Boussinesq equation. The effect of the finite length of the aquifer is considered in detail. An analytical approximation based on a superposition principle yields a very good estimate of the outflow when compared to accurate numerical solutions. An illustration of the new analytical approach to analyze basin-scale field data is used to demonstrate possible field applications of the new solution.

Similitude applied to centrifugal scaling of unsaturated flow

Centrifuge experiments modeling single-phase flow in prototype porous media typically use the same porous medium and permeant. Then, well-known scaling laws are used to transfer the results to the prototype. More general scaling laws that relax these restrictions are presented. For permeants that are immiscible with an accompanying gas phase, model-prototype (i.e., centrifuge model experiment-target system) scaling is demonstrated. Scaling is shown to be feasible for Miller-similar (or geometrically similar) media. Scalings are presented for a more, general class, Lisle-similar media, based on the equivalence mapping of Richards' equation onto itself. Whereas model-prototype scaling of Miller-similar media can be realized easily for arbitrary boundary conditions, Lisle-similarity in a finite length medium generally, but not always, involves a mapping to a moving boundary problem. An exception occurs for redistribution in Lisle-similar porous media, which is shown to map to spatially fixed boundary conditions. Complete model-prototype scalings for this example are derived.