Time effects on dynamic soil properties under random excitation conditions

Two types of clayey soils, a kaolinite and a bentonite, were tested using a resonant column apparatus under random excitation conditions. The concept of root mean square (rms) strain was utilized for the purpose of strain calculations during random loading. The conventional estimator of the transfer function was used for random vibration analysis. The effect of confinement duration (at a constant pressure) on dynamic soil properties, namely damping and shear modulus, was evaluated. The results indicate that for both cohesive soils, the effect of time was less pronounced during random vibration than sinusoidal loading at the same rms strain. This effect is however more pronounced when the peak shearing strain of sinusoidal loading is considered. Furthermore, time effects were more pronounced at low strain levels than at high strain levels.     

River channel adjustments downstream from channelization works in England and Wales

Relatively little attention has been given to river channel adjustments that occur downstream from channelization works. This study is concerned with the nature of channel adjustments downstream from a total of 46 channelization works located in low and high energy environments in England and Wales. Channel changes are identified principally by the method of field survey and by reconstructing the original positions of eroded beds and banks. Use is also made of maps, aerial photographs, and engineering drawings of different dates and the technique of space-for-time substitution is applied. Enlargement of channel cross-sections through erosion had occurred downstream from a variety of types, sizes, and dates of channelization works. The maximum increase of channel size was 153 per cent. Out of a total of 14 sites with enlarged channel cross-sections, seven had undergone a change of width only, at a further three width increased rather than depth, and at the remaining four sites depth increases were dominant. These sites all have relatively high stream powers. Factors causing spatial variation of erosion included tree roots locally binding bank sediments and the occurrence of bends. Planform change had taken place at only one site. A further three high stream power sites had downstream reaches incised into bedrock and therefore did not exhibit adjustment. Channel enlargement is explained in terms of increased flood flows downstream from channelization works causing higher stream velocities, which in turn cause erosion, thereby increasing channel width and/or depth. Examination of flow records for 35 stations revealed flood events which would formerly have spread overbank but are now confined by the channelization works and are therefore likely to alter downstream flows. At sites with downstream change it is proposed that the energy of increased flows was sufficient to exceed a threshold required for erosion of perimeter sediments. By contrast the absence of change at a majority of sites in low energy lowland areas could be a reflection of both the incompetence of increased flows to erode and resistance provided by perimeter sediments. Sites with erosion features appear not to have yet attained new equilibrium conditions.     

Inferred ground motions on Guadeloupe during the 2004 Les Saintes earthquake

Accurate estimates of the ground motions that occurred during damaging earthquakes are a vital part of many aspects of earthquake engineering, such as the study of the size and cause of the uncertainties within earthquake risk assessments. This article compares a number of methods to estimate the ground shaking that occurred on Guadeloupe (French Antilles) during the 21st November 2004 (M _w 6.3) Les Saintes earthquake, with the aim of providing more accurate shaking estimates for the investigation of the sources of uncertainties within loss evaluations, based on damage data from this event. The various techniques make differing use of the available ground-motion recordings of this earthquake and by consequence the estimates obtained by the different approaches are associated with differing uncertainties. Ground motions on the French Antilles are affected by strong local site effects, which have been extensively investigated in previous studies. In this article, use is made of these studies in order to improve the shaking estimates. It is shown that the simple methods neglecting the spatial correlation of earthquake shaking lead to uncertainties similar to those predicted by empirical ground-motion models and that these are uniform across the whole of Guadeloupe. In contrast, methods (such as the ShakeMap approach) that take account of the spatial correlation in motions demonstrate that shaking within roughly 10 km of a recording station (covering a significant portion of the investigated area) can be defined with reasonable accuracy but that motions at more distant points are not well constrained.     

Simulation of the dissolution of weathered versus unweathered limestone in carbonic acid solutions of varying strength

A simulation was undertaken within a climatic chamber to investigate limestone dissolution under varied carbonic acid (H₂CO₃) strengths as a possible analogue for future increases in atmospheric CO₂ arising from global warming. Twenty‐eight samples cut from a block of Bath (Box Hill) limestone from Somerville College, Oxford, which had been removed during restoration after 150 years in an urban environment, were weighed and placed in closed bottles of thin plastic containing varying concentrations of H₂CO₃. Half of the stone samples were derived from exposed surfaces of the stone block (weathered) while the others were obtained from the centre of the block on unexposed surfaces (unweathered). The purpose of this was to compare dissolution of previously weathered versus unweathered surfaces in strong (pH 4·73) versus weak (pH 6·43) solutions of H₂CO₃. A temperature of c. 19 °C was maintained within the chamber representing a plausible future temperature in Oxford for the year 2200 given current warming scenarios. The simulation lasted 25 days with a few stone samples being removed midway. Stone samples show reduced weight in all cases but one. There was greater dissolution of stone samples in a strong H₂CO₃ solution as conveyed by higher concentrations of total hardness and Ca²⁺ in the water samples as well as enhanced microscopic dissolution features identified using SEM. The simulation confirms that enhanced atmospheric CO₂ under global warming, given adequate moisture, will accelerate dissolution rates particularly of newly replaced limestone building stones. However, previously weathered surfaces, such as those on historical stone exposed for a century or more, appear to be less susceptible to the effects of such increased rainfall acidity. Conservation techniques which remove weathered surfaces, such as stone cleaning, may accelerate future decay of historical limestone structures by increasing their susceptibility to dissolution. Copyright © 2006 John Wiley & Sons, Ltd.     

A hyperspectral experiment over tropical forests based on the EO-1 orbit change and PROSAIL simulation

ABSTRACT

We designed a unique hyperspectral experiment from the Earth Observing One (EO-1) orbit change to evaluate solar illumination effects over tropical forests in Brazil. Ten nadir-viewing Hyperion images collected over a fixed site and period of the year (July to August) were selected for analysis. We evaluated variations in reflectance and in 16 narrowband vegetation indices (VIs) with increasing solar zenith angle (SZA) from the pre-drift (2004–2008) to the EO-1 drift period (2011–2016). To detect changes in reflectance and shadows, we applied spectral mixture analysis (SMA) and principal component analysis (PCA) and calculated the similarity spectral angle (θ) between the vegetation spectra measured with variable SZA. The magnitude of the illumination effects was also evaluated from change-point analysis and nonparametric Mann-Whitney U tests applied over the time series. Finally, we complemented our experiment using the PROSAIL model to simulate the VIs variation with increasing SZA resultant from satellite drift. The results showed significant changes in Hyperion reflectance and VIs, especially when the EO-1 crossed the study area at earlier times and larger SZA in 2015 (9:05 a.m.; SZA = 59°) and 2016 (8:30 a.m.; SZA = 67°). Compared to the pre-drift period (10:30 a.m.; SZA = 45°), the SZA differences of 14° (2015) and 22° (2016) increased the shade fractions and decreased the vegetation brightness. PCA separated the pre-drift and drift reflectance datasets, showing shifts in scores due to changes in brightness. θ increased with SZA, indicating changes in the shape of the vegetation spectra with drift. For most VIs, the change-point analysis indicated 2015 (SZA = 59°) as the predominant year of detected changes. Compared to the EO-1 original orbit, the Plant Senescence Reflectance Index (PSRI), Anthocyanin Reflectance Index (ARI) and Structure Insensitive Pigment Index (SIPI) presented the largest positive changes during drift, while the Photochemical Reflectance Index (PRI), Visible Atmospherically Resistant Index (VARI) and Enhanced Vegetation Index (EVI) had the largest negative changes. The effect size of the illumination geometry on these VIs was large, as indicated by increasing values of the Cohen’s r metric toward 2016. The anisotropy of the Hyperion VIs was generally consistent with that from PROSAIL in the simulated pre-drift and drift periods. Focusing on structural indices, it affected the relationships between VIs and simulated leaf area index (LAI) at large SZA.     

Transport of water in frozen soil IV. Analysis of experimental results on the effects of ice content

Effects of ice content on the transport of water in frozen soil are studied experimentally and theoretically under isothermal conditions. A physical law, that the flux of water in unsaturated frozen soil is proportional to the gradient of total water content is proposed. Theoretical justification is made by the use of the two-phase flow theory. The experimental results are shown to support the proposed physical law. The results of this study are presented in two parts and this is the second paper describing the theoretical aspects of the study.     

Multi‐variable relations for soil effects on seismic ground motion

Soil effects on peak ground acceleration, velocity and elastic response spectra (5% damping) are expressed by simple approximate relations in terms of five key parameters: (a) the fundamental vibration period of the non‐linear soil, T_S, (b) the period of a bedrock site of equal thickness, T_b, (c) the predominant excitation period, Te, (d) the peak seismic acceleration at outcropping bedrock, a, and (e) the number of significant excitation cycles, n. Furthermore, another relation is proposed for the estimation of T_S in terms of the soil thickness H, the average shear wave velocity of the soil V?_S,o and a. The aforementioned parameters were first identified through a simplified analytical simulation of the site excitation. The multivariable approximate relations were then formulated via a statistical analysis of relevant data from more than 700 one‐dimensional equivalent‐linear seismic ground response analyses, for actual seismic excitations and natural soil conditions. Use of these relations to back‐calculate the numerical results in the database gives an estimate of their error margin, which is found to be relatively small and unbiased. The proposed relations are also independently verified through a detailed comparison with strong motion recordings from seven well‐documented case studies: (a) two sites in the San Fernando valley during the Northridge earthquake, and (b) five different seismic events recorded at the SMART‐1 accelerometer array in Taiwan. It is deduced that the accuracy of the relations is comparable to that of the equivalent‐linear method. Hence, they can be readily used as a quick alternative for routine applications, as well as for spreadsheet computations (e.g. GIS‐aided seismic microzonation studies) where numerical methods are cumbersome to implement. Copyright © 2003 John Wiley & Sons, Ltd.     

Saturation of the Large Aperture Scintillometer

The saturation aspects of a large aperture (0.3 m) scintillometer operating over a 10-km path were investigated. Measurements were made over mainly forested, hilly terrain with typical maximum sensible heat fluxes of 300–400 W m ⁻², and over flat terrain with mainly grass, and typical maximum heat fluxes of 100–150 W m⁻². Scintillometer-based fluxes were compared with eddy-correlation observations. Two different schemes for calculating the reduction of scintillation caused by saturation were applied: one based on the work of Hill and Clifford, the other based on Frehlich and Ochs. Without saturation correction, the scintillation fluxes were lower than the eddy-correlation fluxes; the saturation correction according to Frehlich and Ochs increased the scintillometer fluxes to an unrealistic level. Correcting the fluxes after the theory of the Hill and Clifford gave satisfying results     

The topographic bias by analytical continuation in physical geodesy

This study emphasizes that the harmonic downward continuation of an external representation of the Earth’s gravity potential to sea level through the topographic masses implies a topographic bias. It is shown that the bias is only dependent on the topographic density along the geocentric radius at the computation point. The bias corresponds to the combined topographic geoid effect, i.e., the sum of the direct and indirect topographic effects. For a laterally variable topographic density function, the combined geoid effect is proportional to terms of powers two and three of the topographic height, while all higher order terms vanish. The result is useful in geoid determination by analytical continuation, e.g., from an Earth gravity model, Stokes’s formula or a combination thereof.     

Comparison of Galerkin-type discretization techniques for two-phase flow in heterogeneous porous media

A comparison of Standard Galerkin, Petrov-Galerkin, and Fully-Upwind Galerkin methods for the simulation of two-phase flow in heterogeneous porous media is presented. On the basis of the coupled pressure-saturation equations, a generalized formulation for all three finite element methods is derived and analysed. For flow in homogeneous media, the Petrov-Galerkin method gives excellent results. But this method fails miserably for problems with heterogeneous media. This is because it is not able to capture correctly processes that take place at interfaces when, for instance, the capillary pressure-saturation relationship after Brooks and Corey is assumed. The Fully-Upwind Galerkin method is superior to the Petrov-Galerkin approach because it is able to give correct results for flow in homogeneous and heterogeneous media for the two models of van Genuchten and Brooks-Corey. The widely used formulation which is correct for the homogeneous case cannot be used for heterogeneous media. Instead the straightforward approach of gradp_c in combination with a chord-slope technique must be utilized.