Determination of the regional deformation rates of Shanghai and Kashima VLBI stations based on ITRF97

The vertical deformation rates (VDRs) and horizontal deformation rates (HDRs) of Shanghai VLBI station in China and Kashima and Kashima34 VLBI stations in Japan were re-analysed using the baseline length change rates from Shanghai to 13 global VLBI stations, and from Kashima to 27 stations and from Kashima34 to 12 stations, based on the NASA VLBI global solution glb1123 (Ma, 1999). The velocity vectors of the global VLBI stations were referred to the ITRF97 reference frame, and the Eulerian vectors of different models of plate motion were used for comparative solutions. The VDR of Shanghai station is estimated to be −1.91±0.56 mm/yr, and those of Kashima and Kashima34 stations, −3.72±0.74mm/yr and −8.81±0.84mm/yr, respectively. The difference between the last two was verified by further analysis. Similar estimates were also made for the Kokee, Kauai and MK_VLBA VLBI stations in mid-Pacific.     

Isotope hydrological study of mean transit time in the granitic Strengbach catchment (Vosges massif,France): application of the FlowPC model with modified input function

Measurements of ¹⁸O concentrations in precipitation, soil solution, spring and runoff are used to determine water transit time in the small granitic Strengbach catchment (0·8 km²; 883–1146 m above sea level) located in the Vosges Mountains of northeastern France. Water transit times were calculated by applying the exponential, exponential piston and dispersion models of the FlowPC program to isotopic input (rainfall) and output (spring and stream water) data sets during the period 1989–95. The input function of the model was modified compared with the former version of the model and estimated by a deterministic approach based on a simplified hydrological balance. The fit between observed and calculated output data showed marked improvements compared with results obtained using the initial version of the model. An exponential piston version of the model applied to spring water indicates a 38·5 month mean transit time, which suggests that the volume in the aquifer, expressed in water depth, is 2·4 m. A considerable thickness (>45 m) of fractured bedrock may be involved for such a volume of water to be stored in the aquifer. Copyright © 2005 John Wiley & Sons, Ltd.     

A critical state model for sands dependent on stress and density

An elastoplastic model for sands is presented in this paper, which can describe stress–strain behaviour dependent on mean effective stress level and void ratio. The main features of the proposed model are: (a) a new state parameter, which is dependent on the initial void ratio and initial mean stress, is proposed and applied to the yield function in order to predict the plastic deformation for very loose sands; and (b) another new state parameter, which is used to determine the peak strength and describe the critical state behaviour of sands during shearing, is proposed in order to predict simply negative/positive dilatancy and the hardening/softening behaviour of medium or dense sands. In addition, the proposed model can also predict the stress–strain behaviour of sands under three-dimensional stress conditions by using a transformed stress tensor instead of ordinary stress tensor. Copyright © 2004 John Wiley & Sons, Ltd.     

Plane of elastic non-resisting tension material under foundation structures

In the frame of 2D-static problems one approaches the problem of elastic-NRT (not-resisting tension) semi-plane loaded on its limit line. This problem is intended to model the stress situation induced in the soil by a foundation structure. The solution, in terms of activated stress field, is searched for in the class of stress fields satisfying equilibrium and admissibility conditions, by applying an energy approach. Copyright © 2004 John Wiley & Sons, Ltd.     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Use of multi‐platform,multi‐temporal remote‐sensing data for calibration of a distributed hydrological model: an application in the Arno basin,Italy

Images from satellite platforms are a valid aid in order to obtain distributed information about hydrological surface states and parameters needed in calibration and validation of the water balance and flood forecasting. Remotely sensed data are easily available on large areas and with a frequency compatible with land cover changes. In this paper, remotely sensed images from different types of sensor have been utilized as a support to the calibration of the distributed hydrological model MOBIDIC, currently used in the experimental system of flood forecasting of the Arno River Basin Authority. Six radar images from ERS‐2 synthetic aperture radar (SAR) sensors (three for summer 2002 and three for spring–summer 2003) have been utilized and a relationship between soil saturation indexes and backscatter coefficient from SAR images has been investigated. Analysis has been performed only on pixels with meagre or no vegetation cover, in order to legitimize the assumption that water content of the soil is the main variable that influences the backscatter coefficient. Such pixels have been obtained by considering vegetation indexes (NDVI) and land cover maps produced by optical sensors (Landsat‐ETM). In order to calibrate the soil moisture model based on information provided by SAR images, an optimization algorithm has been utilized to minimize the regression error between saturation indexes from model and SAR data and error between measured and modelled discharge flows. Utilizing this procedure, model parameters that rule soil moisture fluxes have been calibrated, obtaining not only a good match with remotely sensed data, but also an enhancement of model performance in flow prediction with respect to a previous calibration with river discharge data only. Copyright © 2006 John Wiley & Sons, Ltd.     

A qualitative and quantitative evaluation of experimental model catchment evolution

Due to the geological time scales required for observation of catchment evolution, surrogates or analogues of field data are necessary to understand long‐term processes. To investigate long‐term catchment behaviour, two experimental model catchments that developed without rigid boundaries under controlled conditions are examined and a qualitative and quantitative analysis of their evolution is presented. Qualitatively, the experimental catchments have the visual appearance of field scale data. Observation demonstrates that changes in catchment shape and network form are conservative. Quantitative analysis suggests that the catchments reach an equilibrium form while a reduction in the channel network occurs. While the catchments are laboratory scale models, the results provide insights into field scale behaviour. Copyright © 2003 John Wiley & Sons, Ltd.     

Weathering process at the natural fission reactor of Bangombé

 The uranium deposits in the basin of Franceville (Gabon) host the only natural fission reactors known in the world. Unique geological conditions favoured a natural fission reaction 2 Ga ago. This was detected by anomalous isotopic compositions of uranium and rare earth elements (REE), which are produced by the fission reaction. In total, 16 reactor zones were found. Most of them are mined out. The reactor zone of Bangombé, is only 10–11 m below the surface. This site has been influenced by surface weathering processes. Six drill cores have been sampled at the site of the reactor zone of Bangombé during the course of the study and only one drill core (BAX 08) hit the core of the reactor. From these data and previous drilling campaigns, the reactor size is estimated to be 10 cm thick, 2–3 m wide and 4–6 m long. The migration of fission products can be traced by the anomalous isotope ratios of REE because of the fission process. The ¹⁴⁹Sm/¹⁴⁷Sm ratio close to the reactor zone is only 0.28 (normal: 0.92) because of the intense neutron capture of ¹⁴⁹Sm and subsequent transmutation, whereas ¹⁴⁷Sm is enriched by the fission reaction. Similar changes in isotopic patterns are detectable on other REE. The isotope ratios of Sm and Nd of whole rock and fracture samples surrounding the reactor indicate that fission-genic REE migrated only a few decimetres above and mainly below the reactor zone. Organic matter (bitumen) seems to act as a trap for fission-genic REE. Additional REE-patterns show less intense weathering with increasing depth in the log profile and support a simple weathering model. Received: 26 November 1999 · Accepted: 2 May 2000     

A stepwise method for determining the number of component distributions in a mixture

This paper considers the problem of estimatingm, the number of components in a finite mixture of distributions from a parametric family. A step-up procedure using the bootstrap method is proposed. Some properties of the procedure are illustrated with simulation studies. An example of the method, applied to orientation of beach clasts, is given.