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.     

The ellipsoidal variables

We perform a Monte Carlo analysis of the eclipse probabilities of short-period binary systems in an attempt to estimate the distribution of ellipsoidal variable system with Main-Sequence components. Our results suggest that possibly as many as one in four of the non-eclipsing, spectroscopic binaries withP(d)10, and Main-Sequence components withM _v 10 could be ellipsoidal variable systems. This result, while an upper limit, would seem to hold irrespective of primary spectral type (B5 to G5).     

Primordial magnetic fields and symmetry restoration

In the standard cosmological model symmetry breaking in grand unified theories will occur at times 10^–39 s after the initial singularity when the Universe has cooled to a temperature 10¹⁶ GeV. We investigate here whether it is possible for a uniform, large-scale, magnetic field present in the early universe to delay significantly the time at which symmetry breaking occurs. Given the present magnitude of the intergalactic B-field (10^–11–10^–9 G) it is found that no significant effects are introduced.     

Basement–cover relationships of the Kalak Nappe Complex, Arctic Norwegian Caledonides and constraints on Neoproterozoic terrane assembly in the North Atlantic region

The Kalak Nappe Complex (KNC) has been regarded as Baltica passive margin metasediments telescoped eastwards onto the Baltic (Fennoscandian) Shield during the Caledonian Orogeny. Recent studies have questioned this interpretation, instead pointing to a Neoproterozoic exotic origin. In an effort to resolve this controversy we present a Sm–Nd and U–Th–Pb study of gnessic units, traditionally considered as the depositional basement, along with cover rock sediments and intrusives. Late Palaeoproterozoic gneisses now beneath the KNC were deposited after 1948 ± 33 Ma, before intrusion of the Tjukkfjellet Granite at 1796 ± 3 Ma, and were affected by later melting events at 1765 ± 9 and 1727 ± 9 Ma. These gneisses are interpreted as part of the Baltic Shield and underlie the KNC across a tectonic contact. An unconformity between psammites of the KNC and other paragneisses previously considered as its Precambrian basement is reinterpreted as a modified sedimentary contact between Neoproterozoic metasediments. These metasediments have statistically very similar detrital zircon populations with grains as young as 1034 ± 22, 1025 ± 32 and 1014 ± 14 Ma. The results indicate that the KNC sediments were deposited during the Neoproterozoic in basins along the Laurentian margin of eastern Rodinia and were not connected to Baltica via a depositional basement. Dating of the 851 ± 5 Ma Eidvågvatnet and 853 ± 4 Ma Nordneset granites shows that intrusive material associated with the Porsanger Orogeny (c. 850 Ma) affected a considerable region of the upper KNC terrane. Later Neoproterozoic events at 711 ± 6, 687 ± 12 and 617 ± 6 Ma are also recognised the latest of which may be an expression of rifting. Since early Neoproterozoic magmatism (c. 840–690 Ma) is unknown in Baltica, these results support an exotic origin for the KNC terranes.     

Travel time based approach for the assessment of vulnerability of karst groundwater: the Transit Time Method

Karst aquifers represent important water resources in many parts of the world. Unfortunately, karst aquifers are characterised by high contamination risks. This paper presents a travel time based method for the estimation of karst groundwater vulnerability. It considers (1) physics-based lateral flow within the uppermost weathered zone (epikarst) in a limestone-dominated region and (2) high velocities of vertical infiltration at discrete infiltration points (e.g. sinkholes) or lines (e.g. dry valleys, faults). Consequently, the Transit Time Method honours the actual flow path within the unsaturated zone of a karst aquifer system. A test site in Northern Jordan was chosen for the demonstration of the assessment technique, i.e. the catchment area of the Qunayyah Spring north of the capital Amman. The results demonstrate that zones of highest vulnerability lie within valleys and nearby main fault zones. It also reveals that regions, categorised as protected areas by other methods due to thick unsaturated zones, contribute to a major degree to the total risk.     

A design methodology for rock slopes susceptible to wedge failure using fracture system modelling

This paper demonstrates how the use of fracture system modelling can be linked to limit equilibrium analysis of rock slopes susceptible to wedge failure. The use of fracture systems highlights some of the limitations inherent in traditional structural data analysis and representation. Consequently it allows for more comprehensive input data that can be used for stability analysis of rock slopes. In particular the developed methodology addresses important issues such as spatial variability and wedge size distributions. The paper introduces a series of guidelines for interpretation of the results of rock slopes. The proposed techniques arguably result in an improved level of confidence in the design of rock slopes susceptible to wedge failure.