Numerical viscoelastic modelling by the spectral Laguerre method

Seismic wave propagation in a viscoelastic media can be described by a system of integro-differential equations. The solution of such equations requires special methods when using finite-difference techniques in the time domain. In the frequency domain, the integral terms are represented by complex elastic parameters. This paper presents an efficient algorithm for viscoelastic modelling based on the integral Laguerre transform for the approximation of temporal derivatives and for the calculation of convolution integrals. For the calculation of spatial derivatives, it is possible to use various methods: finite-difference and finite-element techniques, spectral and pseudo-spectral methods. We then obtain a system of algebraic equations with a matrix independent of the parameter m , i.e. the degree of the Laguerre polynomials. In this case, only the right-hand side of the system has recurrent dependence on the parameter m , which is an analogue of the temporal frequency in the frequency domain. The obtained system with a large number of right-hand sides can be solved using fast methods, where the matrix is transformed only once, as opposed to the frequency-domain approach, when the matrix is transformed for each temporal frequency.     

Smart optics in ast

Smart optics means much more than adaptive optics on telescopes: these new technologies are changing the way space instruments are built and operated and are bringing new technologies into everyday life in the form of cheaper, lighter and more robust optical systems. Steve Welch, Peter Doel, Alan Greenaway and Gordon Love summarize the work of the Smart Optics Faraday Partnership in the UK.     

Computer Simulation of Shaped Charge Detonation Using Meshless Particle Method

Shaped charges are widely applied in the defense industry and civil engineering for penetrating hardened targets. Simulations of shaped charges are generally based on traditional grid-based numerical methods. Different from the existing literature, the SPH method is applied in this paper to simulate the detonation process of shaped charge explosive. The combination of the meshless, Lagrangian and particle nature inherent in the SPH method is very attractive in dealing with large deformations and large inhomogeneities in extremely transient situations such as detonation process. The SPH method and some related numerical aspects are briefly discussed first, followed by applications in simulating shaped charge detonations. The effects of different detonation cavity angles and different charge head lengths are investigated. The results demonstrate that major shaped charge detonation phenomena are captured with revealing observations.     

The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

A coupled,dynamical and chemical model for the prestellar core L1544: Comparison of modeled and observed C<Superscript>18</Superscript>O,HCO<Superscript>+</Superscript>, and CS emission spectra

We consider radiative transfer in C¹⁸O, HCO⁺, and CS molecular lines in a spherically symmetrical, coupled, dynamical and chemical model of a prestellar core whose evolution is determined by ambipolar diffusion. Theoretical and observed line profiles are compared for the well-studied core L1544, which may be a collapsing protostellar cloud. We study the relationship between the line shapes and model parameters. The structure of the envelope and kinematic parameters of the cloud are the most important factors determining the shape of the lines. Varying the input model parameters for the radiative transfer—the kinetic temperature and microturbulent velocity—within the limits imposed by observations does not result in any substantial variations of the line profiles. The comparison between the model and observed spectra indicates that L1544 displays a flattened structure, and is viewed at an oblique angle. A two-dimensional model is needed to reproduce this structure.     

Possible refugia for reefs in times of environmental stress

This paper investigates the refuge potential of (1) upwelling areas, (2) coral areas at medium depth, and (3) offshore bank and island reefs in a scenario of increased global warming, and thus increased sea surface temperature (SST) and increased solar UV radiation. (1) Observations on coral health and water temperature in the subtropical Atlantic (Eleuthera and Cat Island, Bahamas) and Indian Ocean (Sodwana Bay, South Africa) suggest a link between cool water delivered by upwelling and coral health. After the 1998 bleaching event, caused by strong SST anomalies, coral health and recovery from the previous year's bleaching was significantly better on the narrow southern Cat Island shelf (70% of corals healthy) where the presence of cold water was observed, which was attributed to small-scale upwelling, than on the wide northern Eleuthera shelf (44% of corals healthy), where downwelling of hot bank waters was believed to have damaged corals. In South Africa, regular, short-term upwelling events in five summers reduced SST to well below bleaching level. (2) In the northern Red Sea (Safaga Bay) and in South Africa (Sodwana Bay), wide areas with either coral frameworks or non-framework communities exist. Calculations show that if the top 10 m (20 m) of the ocean became inhospitable to corals, still 50.4% (17.5%) of the coral area would remain intact in the Red Sea and 99% (40%) in South Africa. (3) Offshore bank and island reefs investigated in the Turks, Caicos, and Mouchoir Banks and Grand and Little Cayman showed high rates of mortality and coral diseases. The most remote sites (Mouchoir Bank) were not the healthiest. Refuge areas appear to exist in (1) and (2), but in (3) only if vigorous water-circulation is encountered.     

On Rhythmical Layering of Rocks Formed from Basaltic Magma

Formation of rhythmical layering in intrusive basic and ultrabasic rock bodies is explained in different ways, in particular, by the movement of microparticles (mineral clusters) of plagioclase and pyroxene within the basalt melt under the influence of thermal and gravitational forces. A model of cluster movement is proposed as a consequence of the forces appearing when ultrasonic elastic fluctuations pass through the melt. The model is based on fundamental dynamic equations. Depending on cluster density, wave parameters, and magmatic chamber size, in the melt there can form different combinations of rhythmically alternating layers of different composition.     

Changes in physical and chemical properties of three soil types in India as a result of amendment with fly ash and sewage sludge

The use of coal fly ash and domestic sewage sludge in agriculture is being considered as one of the methods for recycling of these wastes in an environmental beneficial manner. Mixtures with soil were prepared at different proportions of fly ash and sludge, either alone or in combination at a maximum application rate of 52 t ha^-1. The changes in the selected properties and heavy metal contents of three soil types in India were studied after incubating the respective mixtures for 90 days at near field capacity moisture level. Sewage sludge, due to its acidic and saline nature, high organic matter and heavy metals content, had more impact on soil properties than the fly ash. Sludge application produced several changes including an increase in available nitrogen, organic carbon, salinity and water-holding capacity of the soils. The concentrations of major cations and heavy metals also increased because of the sludge application and the pH was decreased. However, the levels of individual metal concentrations in all the mixture types were below the allowable limits prescribed by several environmental agencies. Using fly ash either alone or in equal quantity with sewage sludge had little influence on soil properties and heavy metal content. The relative availability (RA) of heavy metals in three soils amended with 52 t ha^-1 of sewage sludge was observed to be highest in oxisol, followed by alfisol and vertisol.