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     

Fracture energy release and size effect in borehole breakout

The paper presents a simple approximate analytical solution of the remote stresses that cause the collapse of a borehole or other circular cylindrical cavity in an infinite elastic space. Regions of parallel equidistant splitting cracks are assumed to form on the sides of the cavity. Their boundary is assumed to be an ellipse of a growing horizontal axis, the other axis remaining equal to the borehole diameter. The slabs of rock between the splitting cracks are assumed to buckle as slender columns, and their post-critical stress is considered as the residual stress in the cracked rock. The buckling of these slab columns is assumed to be resisted not only by their elastic bending stiffness but also shear stresses produced on rough crack faces by relative shear displacements. The energy release from the infinite medium caused by the growth of the elliptical cracking region is evaluated according to Eschelby's theorem. This release is set equal to the energy dissipated by the formation of all the splitting cracks, which is calculated under the assumption of constant fracture energy. This yields the collapse stress as a function of the elastic moduli, fracture energy, ratio of the remote principal stresses, crack shear resistance characteristic and borehole diameter. The collapse stress as a function of crack spacing is found to have a minimum, and the correct crack spacing is determined from this minimum. For small enough diameters, the crack spacing increases as the (4/5)-power of the borehole diameter, while for large enough diameters a constant spacing is approached. In contrast to plastic solutions, the breakout stress exhibits a size effect, such that for small enough diameters the breakout stress decreases as the (? 2/5)-power of the borehole diameter, while for large enough diameters a constant limiting value is approached. Finally, some numerical estimates are given and the validity of various simplifying assumptions made is discussed.     

Ice‐cover variability on shallow lakes at high latitudes: model simulations and observations

A one‐dimensional thermodynamic model for simulating lake‐ice phenology is presented and evaluated. The model can be driven with observed daily or hourly atmospheric forcing of air temperature, relative humidity, wind speed, cloud amount and snowfall. In addition to computing the energy balance components, key model output includes the temperature profile at an arbitrary number of levels within the ice/snow (or the water temperature if there is no ice) and ice thickness (clear ice and snow‐ice) on a daily basis, as well as freeze‐up and break‐up dates. The lake‐ice model is used to simulate ice‐growth processes on shallow lakes in arctic, sub‐arctic, and high‐boreal forest environments. Model output is compared with field and remote sensing observations gathered over several ice seasons. Simulated ice thickness, including snow‐ice formation, compares favourably with field measurements. Ice‐on and ice‐off dates are also well simulated when compared with field and satellite observations, with a mean absolute difference of 2 days. Model simulations and observations illustrate the key role that snow cover plays on the seasonal evolution of ice thickness and the timing of spring break‐up. It is also shown that lake morphometry, depth in particular, is a determinant of ice‐off dates for shallow lakes at high latitudes. Copyright © 2003 John Wiley & Sons, Ltd.     

A micromechanical finite element model for linear and damage‐coupled viscoelastic behaviour of asphalt mixture

This study presents a finite element (FE) micromechanical modelling approach for the simulation of linear and damage‐coupled viscoelastic behaviour of asphalt mixture. Asphalt mixture is a composite material of graded aggregates bound with mastic (asphalt and fine aggregates). The microstructural model of asphalt mixture incorporates an equivalent lattice network structure whereby intergranular load transfer is simulated through an effective asphalt mastic zone. The finite element model integrates the ABAQUS user material subroutine with continuum elements for the effective asphalt mastic and rigid body elements for each aggregate. A unified approach is proposed using Schapery non‐linear viscoelastic model for the rate‐independent and rate‐dependent damage behaviour. A finite element incremental algorithm with a recursive relationship for three‐dimensional (3D) linear and damage‐coupled viscoelastic behaviour is developed. This algorithm is used in a 3D user‐defined material model for the asphalt mastic to predict global linear and damage‐coupled viscoelastic behaviour of asphalt mixture. For linear viscoelastic study, the creep stiffnesses of mastic and asphalt mixture at different temperatures are measured in laboratory. A regression‐fitting method is employed to calibrate generalized Maxwell models with Prony series and generate master stiffness curves for mastic and asphalt mixture. A computational model is developed with image analysis of sectioned surface of a test specimen. The viscoelastic prediction of mixture creep stiffness with the calibrated mastic material parameters is compared with mixture master stiffness curve over a reduced time period. In regard to damage‐coupled viscoelastic behaviour, cyclic loading responses of linear and rate‐independent damage‐coupled viscoelastic materials are compared. Effects of particular microstructure parameters on the rate‐independent damage‐coupled viscoelastic behaviour are also investigated with finite element simulations of asphalt numerical samples. Further study describes loading rate effects on the asphalt viscoelastic properties and rate‐dependent damage behaviour. Copyright © 2006 John Wiley & Sons, Ltd.     

Wavelet analysis of inter‐annual variability in the runoff regimes of glacial and nival stream catchments,Bow Lake,Alberta

Continuous wavelet analyses of hourly time series of air temperature, stream discharge, and precipitation are used to compare the seasonal and inter‐annual variability in hydrological regimes of the two principal streams feeding Bow Lake, Banff National Park, Alberta: the glacial stream draining the Wapta Icefields, and the snowmelt‐fed Bow River. The goal is to understand how water sources and flow routing differ between the two catchments. Wavelet spectra and cross‐wavelet spectra were determined for air temperature and discharge from the two streams for summers (June–September) 1997–2000, and for rainfall and discharge for the summers of 1999 and 2000. The diurnal signal of the glacial runoff was orders of magnitude higher in 1998 than in other years, indicating that significant ice exposure and the development of channelized glacial drainage occurred as a result of the 1997–98 El Niño conditions. Early retreat of the snowpack in 1997 and 1998 led to a significant summer‐long input of melt runoff from a small area of ice cover in the Bow River catchment; but such inputs were not apparent in 1999 and 2000, when snow cover was more extensive. Rainfall had a stronger influence on runoff and followed quicker flow paths in the Bow River catchment than in the glacial catchment. Snowpack thickness and catchment size were the primary controls on the phase relationship between temperature and discharge at diurnal time scales. Wavelet analysis is a fast and effective means to characterize runoff, temperature, and precipitation regimes and their interrelationships and inter‐annual variability. The technique is effective at identifying inter‐annual and seasonal changes in the relative contributions of different water sources to runoff, and changes in the time required for routing of diurnal meltwater pulses through a catchment. However, it is less effective at identifying changes/differences in the type of the flow routing (e.g. overland flow versus through flow) between or within catchments. Copyright © 2003 John Wiley & Sons, Ltd.     

Delamerian‐aged metamorphism in the southern Curnamona Province,Australia: implications for the evolution of the Mesoproterozoic Olarian Orogeny

Monazite electron microprobe U–Th–Pb and garnet Sm–Nd isotopic data from metapelitic assemblages in the Willyama Supergroup in the southern Curnamona Province, south‐central Australia, indicate that the terrain underwent regional greenschist to amphibolite‐grade metamorphism during the c. 500 Ma Delamerian Orogeny. The Delamerian‐aged mineral assemblages include prograde garnet–staurolite and kyanite‐bearing associations that overprint andalusite‐ and sillimanite‐bearing assemblages that are interpreted to have developed during the c. 1600 Ma Olarian Orogeny. Importantly, the development of secondary kyanite‐bearing assemblages in the southern Curnamona Province has been used previously to suggest that the Olarian Orogeny followed an anticlockwise P–T evolution. If such assemblages are the product of c. 500 Ma metamorphism, then the anticlockwise P–T path is an apparent path, due to the overprint of a distinct metamorphic cycle c. 1100 Ma later. Making such distinctions is therefore extremely important when using the textural and metamorphic evolution of polycyclic terrains to model the thermal behaviour of the crust during orogeny. This study highlights the utility of in situ geochronology, linking age data to petrologically important phases and assemblages.     

Radon variations in a Hungarian village

 A steady radon exhalation is assumed in most publications. In a village of North-East Hungary, however, high radon concentrations have been measured, differing strongly in neighbouring houses and varying in time, due to the interplay of geochemical phenomena. Received: 20 November 1995 · Accepted: 18 June 1996     

Rock warming and drying under simulated intertidal conditions,part II: weathering and biological influences on evaporative cooling and near‐surface micro‐climatic conditions as an example of biogeomorphic ecosystem engineering

The way in which rocks and engineering materials heat‐up and dry‐out in the intertidal zone is of relevance to both weathering and ecology. These behaviours can be measured in the laboratory under controlled conditions designed to replicate those occurring in the field. Previous studies have demonstrated differences in thermal behaviours between rock types and through time as a result of soiling in terrestrial environments, but the influence of weathering and colonization on rock behaviours in the intertidal zone has not been previously assessed. We measured the warming and drying of blocks of rock (limestone and granite) and marine concrete during ‘low‐tide’ events simulated in the laboratory, before and after a period of exposure (eight months) on rock platforms in Cornwall, UK. As well as differences between the material types, temperatures of control (unexposed) and field‐exposed blocks differed in the order of 1 to 2 °C. Drying behaviours were also different after field exposure. Differences during the first few hours of exposure to air and heat were attributed to discolouration and albedo effects. Over longer periods of time, changes in the availability of near‐surface pore water as a result of micro‐scale bioerosion of limestone and the development of bio‐chemical crusts on marine concrete [observed using scanning electron microscopy (SEM)] are suggested as mechanisms enhancing and reducing, respectively, the efficiency of evaporative cooling. The retention of moisture by epilithic biofilms may also influence thermal and drying behaviours of granite. These observations represent one of the first examples of cross‐scalar biogeomorphic linkages in the intertidal zone. The significance of the results for the subsequent efficiency of weathering, and near‐surface micro‐climatic conditions experienced by colonizing organisms is discussed. The involvement of microorganisms in the creation of more (or less) ecologically stressful conditions through the alteration of substratum geomorphic properties and behaviours is suggested as an example of ‘biogeomorphic ecosystem engineering’. Copyright © 2011 John Wiley & Sons, Ltd.