Analyses of trace elements on quartz surfaces in sulfidic mine tailings from Kristineberg (Sweden) a few years after remediation

Laser ablation inductively coupled plasma mass spectrometry (LA-ICP-MS) has been applied to determine the elemental composition of the surface layer, as well as of the first interior layer, of quartz grains from the mine tailings from Kristineberg (northern Sweden) in order to determine concentration gradients between these two layers. The quartz grains were collected from the oxidized and unoxidized zones within the tailings. The aim of this study is to assess the role of quartz surfaces as sites for the attenuation of solutes from the mine-tailings porewater. Concentrations of Cu, Ag, Sb, Pb and Bi are highest near the surface of each grain and decrease towards the interior. The surface concentration of Cu, Zn and Pb is more pronounced within the unoxidized than within the oxidized zone of the tailings. Cu exhibits a distinct concentration peak at the surface of the quartz grains below the pre-remediation oxidation front. For Zn and Ce the trend of high surface concentration is less pronounced than for Cu or Pb. Silver, Bi and As are preferably adsorbed within the uppermost layers of the oxidized zone where the pH is as high as 6.2. The conversion of intensity signals of the elements to concentration values in ppm was done by using external standards (NIST silicate glass).     

Using UCODE_2005 and PHREEQC to determine thermodynamic constants from experimental data

This article presents a method for estimating chemical thermodynamic constants from experimental data using the two computer programs UCODE_2005 and PHREEQC. As an example, the conditional stability constant for lead (Pb) complexation by a remediation agent (carboxymethyl-beta-cyclodextrin) is estimated, but the method can be applied to estimate other thermodynamic parameters such as sorption constants and degradation rate constants. Advantages of this technique include estimation of uncertainties associated with estimated parameters, evaluation of information content of observations, statistical evaluation of the appropriateness of the conceptual model, and statistical-based comparison of different models.     

Remotely sensed assessment of water quality levels in the Pearl River Estuary, China

In this paper, a method of assessing water quality from satellite data is introduced. The composite pollution index (CPI) was calculated from measured chemical oxygen demand (COD) and nutrient concentration. The relationships between CPI and 240 band combinations of SeaWiFS water-leaving radiance were analyzed and the optimal band combination for estimating CPI was chosen from the 240 band combinations. An algorithm for retrieval of CPI was developed using the optimal band combination, (L₄₄₃ × L₅₁₀)/(L₄₁₂ + L₄₉₀). The CPI was estimated from atmospherically corrected SeaWiFS data by employing the algorithm. Furthermore, the CPI value range for each water quality level was determined based on data obtained from 850 samples taken in the Pearl River Estuary. The remotely sensed CPIs were then transferred to water quality levels and appropriate maps were derived. The remotely sensed water quality level maps displayed a similar distribution of levels based on in situ investigation issued by the State Ocean Administration, China. This study demonstrates that remote sensing can play an important role in water quality assessment.     

The use of press data in the development of a database for rainfall-induced landslides in Tegucigalpa,Honduras, 1980–2005

The capital city of Honduras, Tegucigalpa, suffers from the occurrence of destructive landslides on a yearly basis. During the rainy season, damages to infrastructure as well as injuries, casualties and homeless individuals resulting from landslides are reported in the press. This paper presents the development of a database for rainfall-induced landslides for the period 1980–2005, based on the news reported by two local newspapers. The editions comprehended during the study period have been scrutinized, and articles focusing on landslides, tropical storms, hurricanes, floods and vulnerability of the city have been collected. The interpretation of these archives has allowed the compilation of valuable data of approximately 400 landslides. The analyses of monthly and annual precipitation during the study period show how extreme rainfall events like Hurricane Mitch in October of 1998 have significantly contributed to the initiation of landslides. In addition, the assessment of the slums and neighborhoods affected by landslides during the study period reveals an evident link between the social and physical vulnerability of Tegucigalpa. In order to estimate the reliability of this press-based database, the set of landslides that have been reported as a result of Hurricane Mitch in the press archives has been compared with two inventories based on the interpretation of aerial photographs taken in 1999 and 2001. It is shown that the analysis of the landslide damage left after the hurricane can be enriched with the detailed temporal data provided in the archives and the precise location of these events determined by the aerial photographs. Despite the difficulties faced in the compilation of this database, a good comprehension of the temporal and spatial distribution of landslides in Tegucigalpa has been achieved.     

Optimization routine for identification of model parameters in soil plasticity

The paper presents an optimization routine especially developed for the identification of model parameters in soil plasticity on the basis of different soil tests. Main focus is put on the mathematical aspects and the experience from application of this optimization routine. Mathematically, for the optimization, an objective function and a search strategy are needed. Some alternative expressions for the objective function are formulated. They capture the overall soil behaviour and can be used in a simultaneous optimization against several laboratory tests. Two different search strategies, Rosenbrock's method and the Simplex method, both belonging to the category of direct search methods, are utilized in the routine. Direct search methods have generally proved to be reliable and their relative simplicity make them quite easy to program into workable codes. The Rosenbrock and simplex methods are modified to make the search strategies as efficient and user‐friendly as possible for the type of optimization problem addressed here. Since these search strategies are of a heuristic nature, which makes it difficult (or even impossible) to analyse their performance in a theoretical way, representative optimization examples against both simulated experimental results as well as performed triaxial tests are presented to show the efficiency of the optimization routine. From these examples, it has been concluded that the optimization routine is able to locate a minimum with a good accuracy, fast enough to be a very useful tool for identification of model parameters in soil plasticity. Copyright © 2001 John Wiley & Sons, Ltd.     

The aperiodic broad-band X-ray variability of Cygnus X-3

We study the soft X-ray variability of Cygnus X-3. By combining data from the All-Sky Monitor and Proportional Counter Array instruments on the RXTE satellite with EXOSAT /Medium Energy (ME) detector observations, we are able to analyse the power density spectrum (PDS) of the source from 10⁻⁹ to 0.1 Hz, thus covering time-scales from seconds to years. As the data on the longer time-scales are unevenly sampled, we combine traditional power spectral techniques with simulations to analyse the variability in this range. The PDS at higher frequencies  (≳10⁻³ Hz)  are for the first time compared for all states of this source. We find that it is for all states well described by a power law, with index  ∼−2  in the soft states and a tendency for a less steep power law in the hard state. At longer time-scales, we study the effect of the state transitions on the PDS, and find that the variability below  ∼10⁻⁷ Hz  is dominated by the transitions. Furthermore, we find no correlation between the length of a high/soft-state episode and the time since the previous high/soft state. On intermediate time-scales, we find evidence for a break in the PDS at time-scales of the order of the orbital period. This may be interpreted as evidence for the existence of a tidal resonance in the accretion disc around the compact object, and constraining the mass ratio to   M ₂/ M ₁≲ 0.3  .     

Analysis of digital elevation data for palaeosurfaces in south-western Sweden

A 3873-km² large rectangular area of the Precambrian basement in south-western Sweden has been investigated regarding gross morphology of palaeosurfaces. GIS and digital elevation data were used to construct maps of relative and absolute relief and E–W trending height profiles. Three different palaeosurfaces were recognised and delimited; the very even sub-Cambrian peneplain in the east at an intermediate altitudinal position, an uplifted and broken part of the sub-Cambrian peneplain in the centre at the highest present elevation, and in the western coastal areas a joint-aligned valley landscape with high relative relief, probably an exhumed Mesozoic etch-surface. Phanerozoic faulting has dissected the different palaeosurfaces into tectonic blocks, which are shown from trend surface analysis to be tilted in different directions and to different degrees, characteristic and decisive for the morphological development.     

TWO-PHASE OIL MIGRATION IN COMPACTING SEDIMENTARY BASINS MODELLED BY THE FINITE ELEMENT METHOD

The upstream-weighted finite element method with lumped mass matrix is applied to the modelling of oil migration in compacting sedimentary basins. An implicit formulation is made in Lagrangian co-ordinates of a pressure, saturation and a temperature equation, which is based on immiscible two-phase flow of oil and water. The formulation accounts for the compaction of the sediments, the generation of oil from solid organic material (kerogen), the eventual pore space generated by kerogen breakdown, and the density variations of the fluids which may set up thermal convection. The model is validated by comparison with results from a one-dimensional (1D) fractional flow-based migration model. A 2D case example showing oil expulsion from source rocks, and the filling of a trap is presented. The mass balance of the model is easily checked because all oil in the basin originates from breakdown of kerogen. Compared with other alternatives, the simple upstream-weighted finite element method is suggested as a possible first choice for a numerical method for the modelling of oil migration in compacting sedimentary basins. It easily deals with the complex geometry of a basin, it yields reasonably good results, is simple to implement, and the same implementation applies to all spatial dimensions. © 1997 by John Wiley & Sons, Ltd.     

A METHOD TO CORRECT YIELD SURFACE DRIFT IN SOIL PLASTICITY UNDER MIXED CONTROL AND EXPLICIT INTEGRATION

When applying an explicit integration algorithm in e.g. soil plasticity, the predicted stress point at the end of an elastoplastic increment of loading might not be situated on the updated current yield surface. This so-called yield surface drift could generally be held under control by using small integration steps. Another possibility, when circumstances might demand larger steps, is to adopt a drift correction method. In this paper, a drift correction method for mixed control in soil plasticity, under drained as well as undrained conditions, is proposed. By simulating triaxial tests in a Constitutive Driver, the capability and efficiency of this correction method, under different choices of implementation, have been analysed. It was concluded that the proposed drift correction method, for quite marginal additional computational cost, was able to correct successfully for yield surface drift giving results in close agreement to those obtained with a very large number of integration steps. © 1997 by John Wiley & Sons, Ltd.     

Modeling undertow due to random waves

A numerical model of undertow due to random waves is developed. The model includes three sub-models: (i) a model for multi-directional and multi-frequency random wave transformation, (ii) a surface roller evolution model, and (iii) a model for calculating the vertical distribution and the mean value of the undertow velocity. The calculation of wave trough level is performed based on a theory for the wave asymmetry. The model was successfully validated against small- and large-scale laboratory experiments. Thus, the model is expected to provide reliable input for the modeling of sediment transport and morphological change due to waves and currents.