Environmental issues with Australian ports

Ports and shipping operations are under increased environmental scrutiny in Australia, particularly as so many of Australia's ports and shipping routes are located in or near major environmentally sensitive areas. The environmental risks involved have been highlighted by the publicity given to shipping accidents. Environmental regulations over aspects of port operations, including ballast water management, dredge spoil and waste reception, are also becoming stricter. While port authorities have adopted stringent environmental guidelines, they remain concerned that the full impact of proposed new environmental regulations on maritime industry may not be fully understood. This paper discusses some of the issues involved with achieving a balance between environmental protection and the need to recognise the importance of seaborne trade to Australia and to maintain Australia's competitiveness in international trade.     

Determination of palladium and platinum in seaweed

Palladium and platinum concentrations were determined for twenty two species of seaweed from the California coast, U.S.A.. A new sensitive analytical technique using atomic absorption spectrometry was developed for palladium in seaweed. For the first time, palladium contents in seaweed were obtained. The contents of palladium and platinum varied from 0.09 to 0.61 ng g^–1 and 0.25 to 1.75 ng g^–1 in dried material, respectively.The average ratio of Pt/Pd found in seaweed was 3.5, while the ratio of Pt/Pd was 4.5 in seawater. Therefore, it is postulated that there is no discrimination for seaweed to accumulate palladium and platinum from the ambient seawater due to the chemical similarity of the metals.     

On the CO2O2 system in the Northeastern Pacific

The TCO₂, O₂, TA and δ¹³C data of the 1969 Geosecs Intercalibration Cruise was analyzed and found to be consistent with a vertical mixing model which assumes that each point along a vertical profile is a mixture of the upper and lower boundaries. Calculated regression coefficients are in agreement with the model of Redfield et al. (1963) and with the assumption that TA variation is due to carbonate reaction. Oxygen consumption and TCO₂ production decrease exponentially with depth and approximately 80% of ΔCO₂ can be accounted for, on average, by O₂ consumption. The remaining 20% are probably due to carbonate solution which seems to take place at depths below 2,500 m. The present study suggests that the isotopic composition (δ¹³C) of the carbon source, required to account for most of the oxygen consumed, may be heavier than the value of −23%. assigned to dissolved organic carbon and particulate organic carbon.     

Capturing strain localization in reinforced soils

Lade’s single hardening soil model with Cosserat rotation embodied in the finite element method is employed to investigate the behavior of geosynthetic reinforced soils with special attention to the development of shear banding. The ability of the finite element model to detect shear banding in a reinforced soil is examined against three high quality small-scale laboratory plane strain tests on Toyoura sand with and without reinforcement. These three tests were chosen because of the clear failure surfaces that developed in the soil during loading. The FEM analyses were able to reasonably simulate the plane strain laboratory tests including both unreinforced and reinforced cases. The FEM analyses gave reasonably good agreement with the experimental results in terms of global stress–strain relationships and shear band occurrences. Furthermore, and based on FE analyses of a hypothetical geosynthetic reinforced soil (GRS) retaining wall, it is shown that the geosynthetic reinforcements are very effective in hindering the formation of shear bands in GRS retaining walls when small spacing between the reinforcement layers was used. When used properly, the geosynthetic reinforcements made the soil behave as a truly reinforced mass of considerable stiffness and strength.     

Compressibility of titanosilicate melts

The effect of composition on the relaxed adiabatic bulk modulus (K₀) of a range of alkali- and alkaline earth-titanosilicate [X ₂ ^n/n+ TiSiO₅ (X=Li, Na, K, Rb, Cs, Ca, Sr, Ba)] melts has been investigated. The relaxed bulk moduli of these melts have been measured using ultrasonic interferometric methods at frequencies of 3, 5 and 7 MHz in the temperature range of 950 to 1600°C (0.02 Pa s < _s < 5 Pa s). The bulk moduli of these melts decrease with increasing cation size from Li to Cs and Ca to Ba, and with increasing temperature. The bulk moduli of the Li-, Na-, Ca- and Ba-bearing metasilicate melts decrease with the addition of both TiO₂ and SiO₂ whereas those of the K-, Rb- and Cs-bearing melts increase. Linear fits to the bulk modulus versus volume fraction of TiO₂ do not converge to a common compressibility of the TiO₂ component, indicating that the structural role of TiO₂ in these melts is dependent on the identity of the cation. This proposition is supported by a number of other property data for these and related melt compositions including heat capacity and density, as well as structural inferences from X-ray absorption spectroscopy (XANES). The compositional dependence of the compressibility of the TiO₂ component in these melts explains the difficulty incurred in previous attempts to incorporate TiO₂ in calculation schemes for melt compressibility. The empirical relationship KV^-4/3 for isostructural materials has been used to evaluate the compressibility-related structural changes occurring in these melts. The alkali metasilicate and disilicate melts are isostructural, independent of the cation. The addition of Ti to the metasilicate composition (i.e. X₂TiSiO₅), however, results in a series of melts which are not isostructural. The alkaline-earth metasilicate and disilicate compositions are not isostructural, but the addition of Ti to the metasilicate compositions (i.e. XTiSiO₅) would appear, on the basis of modulus-volume systematics, to result in the melts becoming isostructural with respect to compressibility.     

Green's theorem and Green's functions for the steady-state cosmic-ray equation of transport

Green's Theorem is developed for the spherically-symmetric steady-state cosmic-ray equation of transport in interplanetary space. By means of it the momentum distribution functionF _o(r,p), (r=heliocentric distance,p=momentum) can be determined in a regionr _arr_bwhen a source is specified throughout the region and the momentum spectrum is specified on the boundaries atr _a andr _b . Evaluation requires a knowledge of the Green's function which corresponds to the solution for monoenergetic particles released at heliocentric radiusr _o , Examples of Green's functions are given for the caser _a =0,r _b = and derived for the cases of finiter _a andr _b . The diffusion coefficient is assumed of the form = _o(p)r ^b . The treatment systematizes the development of all analytic solutions for steady-state solar and galactic cosmic-ray propagation and previous solutions form a subset of the present solutions.     

Thermal expansion and spontaneous strain associated with the normal-incommensurate phase transition in melilites

The linear thermal expansions of åkermanite (Ca₂MgSi₂O₇) and hardystonite (Ca₂ZnSi₂O₇) have been measured across the normal-incommensurate phase transition for both materials. Least-squares fitting of the high temperature (normal phase) data yields expressions linear in T for the coefficients of instantaneous linear thermal expansion, $$\alpha _1 = \frac{1}{l}\frac{{dl}}{{dT}}$$ for åkermanite: $$\begin{gathered} \alpha _{[100]} = 6.901(2) \times 10^{ - 6} + 1.834(2) \times 10^{ - 8} T \hfill \\ \alpha _{[100]} = - 2.856(1) \times 10^{ - 6} + 11.280(1) \times 10^{ - 8} T \hfill \\ \end{gathered} $$ for hardystonite: $$\begin{gathered} \alpha _{[100]} = 15.562(5) \times 10^{ - 6} - 1.478(3) \times 10^{ - 8} T \hfill \\ \alpha _{[100]} = - 11.115(5) \times 10^{ - 6} + 11.326(3) \times 10^{ - 8} T \hfill \\ \end{gathered} $$ Although there is considerable strain for temperatures within 10° C of the phase transition, suggestive of a high-order phase transition, there appears to be a finite ΔV of transition, and the phase transition is classed as “weakly first order”.     

Vertical motions in an intense magnetic flux tube

The recent discovery of localised intense magnetic fields in the solar photosphere is one of the major surprises of the past few years. Here we consider the theoretical nature of small amplitude motions in such an intense magnetic flux tube, within which the field strength may reach 2 kG. We give a systematic derivation of the governing expansion equations for a vertical, slender tube, taking into account the dependence upon height of the buoyancy, compressibility and magnetic forces. Several special cases (e.g., the isothermal atmosphere) are considered as well as a more realistic, non-isothermal, solar atmosphere. The expansion procedure is shown to give good results in the special case of a uniform basic-state (in which gravity is negligible) and for which a more exact treatment is possible.The form of both pressure and velocity perturbations within the tube is discussed. The nature of pressure perturbations depends upon a critical transition frequency, _p , which in turn is dependent upon depth, field strength, pressure and density in the basic (unperturbed) state of the tube. At a given depth in the tube pressure oscillations are possible only for frequencies greater than _p for frequencies below _p exponentially decaying (evanescent) pressure modes occur. In a similar fashion the nature of motions within the flux tube depends upon a transition frequency, _v . At a given depth within the tube vertically propagating waves are possible only for frequencies greater than _v ; for frequencies below _v exponentially decaying (evanscent) motions occur.The dependence of both _v and _p on depth is determined for each of the special cases, and for a realistic solar atmosphere. It is found that the use of an isothermal atmosphere, instead of a more realistic temperature profile, may well give misleading results.For the solar atmosphere it is found that _v is zero at about 12 km above optical depth ₅₀₀₀= 1, thereafter rising to a maximum of 0.04 s^–1 at some 600 km above ₅₀₀₀ = 1. Below ₅₀₀₀ = 1, in the convection zone, _v has a maximum of 0.013 s^–1. The transition frequency, _p , for the pressure perturbations, is peaked at 0.1 s^–1 just below ₅₀₀₀ = 1, falling to a minimum of 0.02 s^–1 at about one scale-height deeper in the tube     

Evidence linking coronal transients to the evolution of coronal holes

The positions of X-ray coronal transients outside of active regions observed during Skylab were superposed on H synoptic charts and coronal hole boundaries for seven solar rotations. We confirmed a detailed spatial association between the transients and neutral lines. We found that most of the transients were related to large-scale changes in coronal hole area and tended to occur on the borders of evolving equatorial holes.Skylab Solar Workshop Post-Doctoral Appointee, 1975–1977.     

Reservoir description using dynamic parameterisation selection with a combined stochastic and gradient search

There is a correspondence between flow in a reservoir and large scale permeability trends. This correspondence can be derived by constraining reservoir models using observed production data. One of the challenges in deriving the permeability distribution of a field using production data involves determination of the scale of resolution of the permeability. The Adaptive Multiscale Estimation (AME) seeks to overcome the problems related to choosing the resolution of the permeability field by a dynamic parameterisation selection. The standard AME uses a gradient algorithm in solving several optimisation problems with increasing permeability resolution. This paper presents a hybrid algorithm which combines a gradient search and a stochastic algorithm to improve the robustness of the dynamic parameterisation selection. At low dimension, we use the stochastic algorithm to generate several optimised models. We use information from all these produced models to find new optimal refinements, and start out new optimisations with several unequally suggested parameterisations. At higher dimensions we change to a gradient-type optimiser, where the initial solution is chosen from the ensemble of models suggested by the stochastic algorithm. The selection is based on a predefined criterion. We demonstrate the robustness of the hybrid algorithm on sample synthetic cases, which most of them were considered insolvable using the standard AME algorithm.