A one-dimensional model for the interaction between continental-scale ice sheets and atmospheric stationary waves

 The great continental ice sheets of the Pleistocene represented a significant topographic obstacle to the westerly winds in northern midlatitudes. This work explores how consequent changes in the atmospheric stationary wave pattern might have affected the shape and growth of the ice sheets themselves. A one dimensional (1-D) model is developed which permits an examination of the types and magnitudes of the feedbacks that might be expected. When plausible temperature perturbations are introduced at the ice-sheet margin which are proportional to the stationary wave amplitude, the equilibrium shape of the ice sheet is significantly altered, and depends on the sign of the perturbation. The proposed feedback also affects the response of the ice sheet to time-varying climate forcing. The results suggest that the evolution of a continental-scale ice sheet with a land-based margin may be significantly determined by the changes it induces in the atmospheric circulation. Received: 1 October 1999 / Accepted: 17 July 2000     

Methods of array element localization for a towed underwateracoustic array

The location of the hydrophones on a towed underwater acoustic array as a function of time (array element localization) is needed for signal processing. Methods to perform this localization using least squares polynomial fitting to data from depth sensors, heading sensors, and sensors detecting a ping from a single source are discussed. Arc distance along the array is used as the independent parameter so that all solutions are constrained to be space curves. Examples of application to real data are presented, and techniques to discriminate against bad sensor data are discussed     

Metamorphism in Archaean greenstone belts: calculated fluid compositions and implications for gold mineralization

An inescapable consequence of the metamorphism of greenstone belt sequences is the release of a large volume of metamorphic fluid of low salinity with chemical characteristics controlled by the mineral assemblages involved in the devolatilization reactions. For mafic and ultramafic sequences, the composition of fluids released at upper greenschist to lower amphibolite facies conditions for the necessary relatively hot geotherm corresponds to those inferred for greenstone gold deposits (X_CO2= 0.2–0.3). This result follows from the calculation of mineral equilibria in the model system CaO–MgO–FeO–Al₂O₃–SiO₂–H₂O–CO₂, using a new, expanded, internally consistent dataset. Greenstone metamorphism cannot have involved much crustal over-thickening, because very shallow levels of greenstone belts are preserved. Such orogeny can be accounted for if compressive deformation of the crust is accompanied by thinning of the mantle lithosphere. In this case, the observed metamorphism, which was contemporaneous with deformation, is of the low-P high-T type. For this type of metamorphism, the metamorphic peak should have occurred earlier at deeper levels in the crust; i.e. the piezothermal array should be of the ‘deeper-earlier’type. However, at shallow crustal levels, the piezothermal array is likely to have been of ‘deeper-later’type, as a consequence of erosion. Thus, while the lower crust reached maximum temperatures, and partially melted to produce the observed granites, mid-crustal levels were releasing fluids prograde into shallow crustal levels that were already retrograde. We propose that these fluids are responsible for the gold mineralization. Thus, the contemporaneity of igneous activity and gold mineralization is a natural consequence of the thermal evolution, and does not mean that the mineralization has to be a consequence of igneous processes. Upward migration of metamorphic fluid, via appropriate structurally controlled pathways, will bring the fluid into contact with mineral assemblages that have equilibrated with a fluid with significantly lower X_CO2. These assemblages are therefore grossly out of equilibrium with the fluid. In the case of infiltrated metabasic rocks, intense carbonation and sulphidation is predicted. If, as seems reasonable, gold is mobilized by the fluid generated by devolatilization, then the combination of processes proposed, most of which are an inevitable consequence of the metamorphism, leads to the formation of greenstone gold deposits predominantly from metamorphic fluids.     

Response of circular footings and anchor plates in non-homogeneous elastic soils

The axisymmetric elastic response of circular footings and anchor plates in a linearly non-homogeneous elastic soil is analysed. It is assumed that footings/anchors are flexible and subjected to axisymmetric vertical loads. The response of the footing/anchor is modelled by using the classical Poisson–Kirchhoff thin plate theory. A variational technique is used to analyse the interaction problem. A representation for the contact stress is established by using a fundamental solution corresponding to a unit vertical pressure acting over an annular region in the interior of the non-homogeneous soil. The fundamental solution can be derived by using rigorous analytical procedures. The influence of the footing flexibility and the degree of soil non-homogeneity on the displacements, bending moments and contact stresses of a surface footing is examined over a wide range of governing parameters. In the case of anchor plates the influence of depth of embedment, degree of soil non-homogeneity and anchor flexibility on the anchor displacement is investigated.     

Comparison of results from several AOGCMs for global and regional sea-level change 1900–2100

 Sea-level rise is an important aspect of climate change because of its impact on society and ecosystems. Here we present an intercomparison of results from ten coupled atmosphere-ocean general circulation models (AOGCMs) for sea-level changes simulated for the twentieth century and projected to occur during the twenty first century in experiments following scenario IS92a for greenhouse gases and sulphate aerosols. The model results suggest that the rate of sea-level rise due to thermal expansion of sea water has increased during the twentieth century, but the small set of tide gauges with long records might not be adequate to detect this acceleration. The rate of sea-level rise due to thermal expansion continues to increase throughout the twenty first century, and the projected total is consequently larger than in the twentieth century; for 1990–2090 it amounts to 0.20–0.37 m. This wide range results from systematic uncertainty in modelling of climate change and of heat uptake by the ocean. The AOGCMs agree that sea-level rise is expected to be geographically non-uniform, with some regions experiencing as much as twice the global average, and others practically zero, but they do not agree about the geographical pattern. The lack of agreement indicates that we cannot currently have confidence in projections of local sea-level changes, and reveals a need for detailed analysis and intercomparison in order to understand and reduce the disagreements. Received: 1 September 2000 / Accepted: 20 April 2001     

Metamorphic evolution of exhumed middle to lower crustal rocks in the Mojave Extensional Belt, southern California, USA

The Waterman Metamorphic Complex of the central Mojave Desert was exposed as a consequence of early Miocene detachment-dominated extension. However, it has evidence consistent with a more extensive geological history that involves collision of a crustal fragment(s), tectonic thickening by overthrusting and two periods of extension. The metamorphic complex contains granitoid intrusives and felsic mylonitic gneisses as well as polymetamorphic rocks that include marble, calc-silicate, quartzite. mafic granulite, pyribolite, amphibolite, migmatite and biotite schist. The latter group of rocks was affected by an initial series of high-grade metamorphic events (M1 and M2) and a localized lower grade overprint (M3). The initial metamorphism (M1) can be separated into two stages along its high-grade P–T path: M1a, a granulite facies metamorphism at 800–850° C and 7.5–9 kbar and Mlb, an upper amphibolite facies overprint at 750–800° C and 10–12 kbar. M1a developed mineral assemblages and textures consistent with granulite facies conditions at a reduced activity of H₂O and is associated with intense ductile deformation (D1) and minor local partial melting. M1b overprinted the granulite assemblages with a series of hydrous phases under conditions of increasing pressure and H₂O activity and is accompanied by little or no deformation. M2 developed at lower pressures and temperatures (650–750° C, 4.5–5.5 kbar) and is distinguished by a second local overprint of hydrous phases that reflects an input of aqueous fluids probably associated with the intrusion of a series of granitic dykes and veins. Effects of M3 are confined to the Mitchel detachment zone, an anastomosing early Miocene detachment fault, and are characterized by local ductile/brittle deformation (D2) of the pre-existing high-grade rocks and granitoid intrusives and by the production of mylonites and mylonitic gneisses under greenschist facies conditions (300–350° C, 3–5 kbar). The initial overprint (M1a) represents metamorphism, devolatilization and minor partial melting of supracrustal rocks under granulite facies conditions as a consequence of tectonic and, possibly, magmatic thickening. The increasing pressure transition of M1a to M1b reflects a period of continued compressional tectonism, thrusting and influx of H₂O, in part, locally related to crystallization of partial melts. The near isothermal decompression between M1b and M2 probably represents a pre-112-Ma extensional episode that may have been the result of a decompressional readjustment of a thickened crust. Following the initial extensional event, the metamorphic complex remained at depths of 10–17 km for at least 90 Ma until it was uplifted following Miocene extension. M3 develops locally in response to this second extensional period resulting from the early Miocene detachment faulting.     

Multiple-source localization using GPS technology and receivedarrival time structure analysis in an air-deployed system

An efficient and robust method has been developed to locate multiple impulsive sources in an ocean environment. Global position system (GPS) receivers were installed on sonobuoys to obtain their locations within a few meters of accuracy. A sonobuoy field was deployed in a ring-type pattern. Charges were then set off at arbitrary locations within the ring, High-resolution plots were used to obtain direct path and/or first bottom bounce arrivals on each buoy. A model grid of arrival times was constructed, corresponding to the dimensions of the buoy field. A ray model previously developed here at the Applied Research Laboratories at the University of Texas at Austin (ARL:UT) was used to obtain model travel times. The minimum value of the least-square-type error between the real arrival times and the modeled travel times resulted in an unambiguous location of the source, within the limits of the grid spacing chosen. This value was calculated by picking one receiver as the reference and then summing the timing errors of the remaining receivers relative to the reference. Successive iterations with finer grid spacings result in source localization within the accuracy of the buoy locations. The localization routine was extended by allowing permutations of the pulse arrivals on each buoy to account for multiple sources closely separated in time and/or space. An automated correlation technique is presented as an alternative to the leading edge-detection method used here for obtaining relative arrival times. Two proof-of-concept experiments were performed and some results of data obtained at Lake Travis and the Gulf of Mexico are presented     

Kinetic studies of lead(II) uptake by the seagrass Zostera muelleri in water by radiotracing,atomic absorption spectrometry and electrochemical techniques

Evidence was found that uptake of lead from seawater in both model and natural systems by the leaves of the seagrass Zostera muelleri does occur for live, dead and scraped leaves at all the lead concentrations tested. Positive uptake of lead was measured using the three analytical techniques of radio-tracer, differential pulse anodic stripping voltammetry and atomic absorption spectrometry; similar uptake profiles were obtained for each technique. Profiles always showed an initial rapid uptake phase followed by a gradual transition to a plateau phase at which an approximate equilibrium between lead attached to seagrass and free ionic lead in seawater was reached. The presence of active lead uptake processes was indicated in experiments which attempted to remove lead from leaves by a chelating agent.     

Four-element CODAR beam forming

An HF radar called the Coastal Ocean Dynamics Applications Radar (CODAR) is presently being used in several forms to measure ocean surface parameters. The original version was developed by the National Oceanic and Atmospheric Administration (NOAA) and utilizes a four-element receive array. The array consists of four equally spaced elements arranged on a circle with a radius of 0.2151 wavelengths (at 25.4 MHz). It was designed to measure ocean currents using a direction-finding technique based on an extension to a simple two-element interferometer. The problem of determining the bearing of a radiating source can be readily shown to be equivalent to that incurred in spectral estimation. In an attempt to improve upon the processing of existing data, modern nonlinear spectral estimation techniques are applied in a beam-forming bearing estimation procedure and compared against several direction-finding algorithms. Enhancement of bearing estimators via analysis of the eigenstructure of a spatial correlation matrix is included. Antenna response patterns are calculated and used to investigate properties of direction-finding algorithms. Simulated data are used for a comparison of direction finding and beam forming. The asymmetrical bias of each method is investigated to determine its effect on the error in estimating the angle of arrival of a radar target.