Evaluation of liquefaction potential using neural-networks and CPT results

In this research, a reliable Cone Penetration Test data set was gathered with a wide range of parameters. This data was incorporated in a Neural-Networks computer software called STATISTICA Neural-Networks. The back propagation algorithm with a multilayer perceptron network is utilized to analyze the liquefaction occurrence in different sites. In this study, different sets of effective parameters for the neural-network analyses are selected such that to reduce the noise and to obtain more accurate results.Considering the relative importance of effective parameters in liquefaction assessment, it is indicated that σ₀, σ′₀ together play a more important role than what previously was assumed and hence the relative importance of the q_c and seismic parameters are decreased compared with the previous works. The results presented here have more accuracy than previous works while at the same time, the range of the parameters used in this study is much wider than what was previously used. This range of parameters makes the proposed method applicable for practical purposes.     

Modelling of oxygen and nitrogen cycling as a function of macrophyte community in the Thau lagoon

A three-dimensional model coupling physical and biological processes for the whole Thau lagoon (Mediterranean coast of France) was developed in order to assess the relationships between macrophytes and the oxygen and nitrogen cycles. Ten species have been inserted as forcing variables in the model. Plankton dynamics, shellfish cultivation impact and mineralization of organic matter are also considered, as well as nutrient and oxygen exchanges between the sediment and the water column.Simulations with and without the macrophytes have shown that the system can be characterized as having a highly structured pattern involving lagoon nitrogen and oxygen cycles. This pattern is created by the combined influence of macrophytes, watershed and oyster farming.The model has been also used to assess the total annual macrophyte production at the whole lagoon scale. Comparisons with phytoplankton production and with results from other temperate lagoons have underlined the high productivity of the Thau lagoon supported by active nutrient regeneration.     

Drainage in finite-sized unsaturated zones

After the initiation of gravity drainage, water is often assumed to be either (a) draining under unit gradient, or (b) at capillary/gravity equilibrium. Both of these simplifications can be useful, but the regimes of validity of each assumption must be delineated. Water pressures are measured versus time and distance as water drains out of a 1.6 m long sand column to determine the relative effects of capillary and gravitational forces during drainage. For medium sized sands (0.15–0.3 mm in diameter), the capillary pressure is constant in space in a large region of the column for over 12 days, and the water continues to flow under unit gradient for relatively long time scales. Similar results are seen for finer sands, but with a much faster approach to equilibrium. Numerical simulations and analytical estimates are presented and compare favorably to the measurements. Together, the experimental, theoretical and analytical results are used to calculate when capillary/gravity equilibrium is reached as a function of porous media properties and length of the unsaturated zone. The ratio of the length of the unsaturated zone to the bubbling pressure is a key parameter in determining the drainage regime, and that even for relatively short unsaturated zones the equilibrium time scale can be on the order of years.     

A unified formulation of the piecewise exact method for inelastic seismic demand analysis including the P‐delta effect

The non‐linear analysis of single‐degree‐of‐freedom (SDOF) systems provides the essential background information for both strength‐based design and displacement‐based evaluation/design methodologies through the development of the inelastic response spectra. The recursive solution procedure called the piecewise exact method, which is efficiently used for the response analysis of linear SDOF systems, is re‐formulated in this paper in a unified format to analyse the non‐linear SDOF systems with multi‐linear hysteresis models. The unified formulation is also capable of handling the P‐delta effect, which generally involves the negative post‐yield stiffness of the hysteresis loops. The attractiveness of the method lies in the fact that it provides the exact solution when the loading time history is composed of piecewise linear segments, a condition that is perfectly satisfied for the earthquake excitation. Based on simple recursive relationships given for positive, negative and zero effective stiffnesses, the unified form of the piecewise exact method proves to be an extremely powerful and probably the best tool for the SDOF inelastic time‐history and response spectrum analysis including the P‐delta effect. A number of examples are presented to demonstrate the implementation of the method. Copyright © 2003 John Wiley & Sons, Ltd.     

System identification of linear structures based on Hilbert–Huang spectral analysis. Part 2: Complex modes

A method, based on the Hilbert–Huang spectral analysis, has been proposed by the authors to identify linear structures in which normal modes exist (i.e., real eigenvalues and eigenvectors). Frequently, all the eigenvalues and eigenvectors of linear structures are complex. In this paper, the method is extended further to identify general linear structures with complex modes using the free vibration response data polluted by noise. Measured response signals are first decomposed into modal responses using the method of Empirical Mode Decomposition with intermittency criteria. Each modal response contains the contribution of a complex conjugate pair of modes with a unique frequency and a damping ratio. Then, each modal response is decomposed in the frequency–time domain to yield instantaneous phase angle and amplitude using the Hilbert transform. Based on a single measurement of the impulse response time history at one appropriate location, the complex eigenvalues of the linear structure can be identified using a simple analysis procedure. When the response time histories are measured at all locations, the proposed methodology is capable of identifying the complex mode shapes as well as the mass, damping and stiffness matrices of the structure. The effectiveness and accuracy of the method presented are illustrated through numerical simulations. It is demonstrated that dynamic characteristics of linear structures with complex modes can be identified effectively using the proposed method. Copyright © 2003 John Wiley & Sons, Ltd.     

Continuous wavelet-like filter for a spherical surface and its application to localized admittance function on Mars

We have developed a 2D isotropic continuous wavelet-like transform for a spherical surface. The transform is simply defined as the surface convolution between the original field and a kernel, based on the zeroth-order Bessel function with a spherical correction. This spherical correction violates the geometric similarity for the various scales of the kernels, which becomes more apparent at longer wavelengths. We found numerically that this transform is practically equivalent to a Gaussian bandpass filter in the spherical harmonic domain. We have applied this wavelet-like transform on the recently acquired Martian gravity and topography fields. Using a ratio constructed locally from these two fields, we have constructed a map describing the lateral variations of the localized admittance function on Mars.     

Organic litter: dominance over stones as a source of interrill flow roughness on low‐gradient desert slopes at Fowlers Gap,arid western NSW,Australia

Thirty‐six runoff plot experiments provide data on flow depths, speeds, and Darcy–Weisbach friction coefficients (f) on bare soil surfaces, and surfaces to which were added sufficient extra plant litter or surface stones to provide projected cover of 5, 10 and 20 per cent. Precision flow depth data were derived with a computer‐controlled gantry and needle gauge for two different discharges for each plot treatment. Taking a fixed flow intensity (Reynolds number, Re = 150) for purposes of comparison shows means of f = 17·7 for bare soil surfaces, f = 11·4 for added stone treatments, and f = 23·8 for added litter treatments. Many individual values of f for stone treatments are lower than for the bare soil surface, but all litter treatments show increases in fcompared to bare soil. The lowering of f in stone treatments relates to the submerged volume that the stones occupied, and the associated concentration of flow onto a smaller part of the plot surface. This leads to locally higher flow intensities and lower frictional drag along threads of flow that the obstacles create. Litter causes higher frictional drag because the particles are smaller, and, for the same cover fraction, are 100 times more numerous and provide 20 times the edge or perimeter length. Along these edges, which in total exceed 2·5 m g^?1 (equivalent to 500 m m^?2 for a loading of 2 t ha^?1), surface tension draws up water from between the litter particles. This reduces flow depth there, and as a consequence of the lower flow intensity, frictional drag rises. Furthermore, no clear passage remains for the establishment of flow threads. These findings apply to shallow interrill flows in which litter is largely immobile. The key new result from these experiments is that under these conditions, a 20 per cent cover of organic litter can generate interrill frictional retardation that exceeds by nearly 41 per cent that of a bare soil surface, and twice that contributed by the same cover fraction of surface stones. Even greater dominance by litter can be anticipated at the many dryland sites where litter covers exceed those tested here. Copyright © 2002 John Wiley & Sons, Ltd.     

The role of sediment specific gravity and availability on cluster evolution

The objective of this experimental study was to account for the role of sediment availability and specific gravity on cluster formation and cluster geometric characteristics (spacing and size). To isolate the effects of sediment availability and specific gravity on cluster evolution, mono‐sized spheres were used to simulate the cluster evolutionary cycle. Overall, twelve experimental runs were carried out in the laboratory flume. Six of these tests were performed by using glass spheres (specific gravity, SG = 2·58) and the other six by employing an equal combination of glass and Teflon spheres (SG = 2·12) of the same diameter to evaluate the role of specific gravity on cluster evolution. The three sediment availability conditions that were investigated here simulated isolated gravel elements, pool–riffle sequences and densely packed gravel‐bed. An advanced image analysis technique was employed to track the evolution of cluster microforms and provide quantitative information about the size and shape of clusters and the number of clusters per unit area. The results of this study showed that: (1) sediment availability affects the architecture and size of cluster microforms; and (2) clusters consisting of mono‐sized sediments start disintegrating at twice the incipient conditions. By performing complementary tests for the isolated gravel elements case, it was found that the evolutionary cycle of individual clusters could be described as follows, in order of increasing stress: no cluster→two particle cluster→comet→triangle→rhomboid→break up. Copyright © 2003 John Wiley & Sons, Ltd.     

Testing a model for predicting the timing and location of shallow landslide initiation in soil‐mantled landscapes

The growing availability of digital topographic data and the increased reliability of precipitation forecasts invite modelling efforts to predict the timing and location of shallow landslides in hilly and mountainous areas in order to reduce risk to an ever‐expanding human population. Here, we exploit a rare data set to develop and test such a model. In a 1·7 km² catchment a near‐annual aerial photographic coverage records just three single storm events over a 45 year period that produced multiple landslides. Such data enable us to test model performance by running the entire rainfall time series and determine whether just those three storms are correctly detected. To do this, we link a dynamic and spatially distributed shallow subsurface runoff model (similar to TOPMODEL) to an in?nite slope model to predict the spatial distribution of shallow landsliding. The spatial distribution of soil depth, a strong control on local landsliding, is predicted from a process‐based model. Because of its common availability, daily rainfall data were used to drive the model. Topographic data were derived from digitized 1 : 24 000 US Geological Survey contour maps. Analysis of the landslides shows that 97 occurred in 1955, 37 in 1982 and ?ve in 1998, although the heaviest rainfall was in 1982. Furthermore, intensity–duration analysis of available daily and hourly rainfall from the closest raingauges does not discriminate those three storms from others that did not generate failures. We explore the question of whether a mechanistic modelling approach is better able to identify landslide‐producing storms. Landslide and soil production parameters were ?xed from studies elsewhere. Four hydrologic parameters characterizing the saturated hydraulic conductivity of the soil and underlying bedrock and its decline with depth were ?rst calibrated on the 1955 landslide record. Success was characterized as the most number of actual landslides predicted with the least amount of total area predicted to be unstable. Because landslide area was consistently overpredicted, a threshold catchment area of predicted slope instability was used to de?ne whether a rainstorm was a signi?cant landslide producer. Many combinations of the four hydrological parameters performed equally well for the 1955 event, but only one combination successfully identi?ed the 1982 storm as the only landslide‐producing storm during the period 1980–86. Application of this parameter combination to the entire 45 year record successfully identi?ed the three events, but also predicted that two other landslide‐producing events should have occurred. This performance is signi?cantly better than the empirical intensity–duration threshold approach, but requires considerable calibration effort. Overprediction of instability, both for storms that produced landslides and for non‐producing storms, appears to arise from at least four causes: (1) coarse rainfall data time scale and inability to document short rainfall bursts and predict pressure wave response; (2) absence of local rainfall data; (3) legacy effect of previous landslides; and (4) inaccurate topographic and soil property data. Greater resolution of spatial and rainfall data, as well as topographic data, coupled with systematic documentation of landslides to create time series to test models, should lead to signi?cant improvements in shallow landslides forecasting. Copyright © 2003 John Wiley & Sons, Ltd.     

dC Variations in Late Jurassic Carbonates, Jaisalmer Formation, Western India

The carbon isotope measurements, carried out on subsurface carbonate samples from Oxfordian Jaisalmer Formation, western India, yield positive d¹³C values up to +3.17%. The most positive Oxfordian C-isotope value corresponds to the carbon isotope excursion measured in samples from from other late Jurassic basins of world. The latest Oxfordian C-isotope values of Jaisalmer Basin fluctuate around 2% while the C-isotope values of 1.50% mark the base of Kimmeridgian. The Oxfordian C-isotope excursion appears to correspond to a time of overall increased organic carbon burial triggered by increased nutrient transfer from continents to oceans during a time of rising global sea level.