The northern Patagonia Somuncura plateau basalts: a product of slab‐induced,shallow asthenospheric upwelling?

The Meseta de Somuncura forms the largest basaltic plateau (20 000 km²) of southern Argentina (extra-Andean domain). Most of these tholeiitic to alkaline rocks were extruded at ˜ 25 Ma (late Oligocene). The absence of rifting–thinning processes, plume activity, or slab-window phenomena leaves only one major possibility for the generation of Somuncura: asthenospheric ('OIB-like') corner flow leading to a transient thermal anomaly above the subducting plate. It is suggested herein that the intake of hot asthenosphere was forced into a favourable topography (concave-up) of the subducting plate, when a major plate reorganization event (Farallon to Nazca) was taking place in late Oligocene to early Miocene time. The fast and vigorous intake of asthenosphere would have been induced by slab roll-back, leading to decoupling of the subducting plate. The Somuncura volcanic episode can be regarded as a marker of the passage from the extremely oblique subduction of Farallon, to the birth of the Nazca plate and roughly perpendicular convergence between South America and Nazca.     

Reliability assessment of granular filters in embankment dams

Empirical criteria have been used successfully to design filters of most embankment large dam projects throughout the world. However, these empirical rules are only applicable to a particular range of soils tested in laboratory and do not take into account the variability of the base material and filter particle sizes. In addition, it is widely accepted that the safety of fill dams is mainly dependent on the reliability of their filter performance. The work herein presented consists in a new general method for assessing the probability of fulfilling any empirical filter design criteria accounting for base and filter heterogeneity by means of first‐order reliability methods (FORM), so that reliability indexes and probabilities of fulfilling any particular criteria are obtained. This method will allow engineers to estimate the safety of existing filters in terms of probability of fulfilling their design criteria and might also be used as a decision tool on sampling needs and material size tolerances during construction. In addition, sensitivity analysis makes possible to analyse how reliabilities are influenced by different sources of input data. Finally, in case of a portfolio risk assessment, this method will allow engineers to compare the safety of several existing dams in order to prioritize safety investments and it is expected to be a very useful tool to evaluate probabilities of failure due to internal erosion. Copyright © 2006 John Wiley & Sons, Ltd.     

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     

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.