NUMERICAL ANALYSIS OF EMBANKMENT FOUNDATION LIQUEFACTION COUNTERMEASURES

Computational simulations are presented for a unique series of centrifuge tests conducted to assess the performance of liquefaction countermeasure techniques. In these centrifuge tests, the dynamic response of an embankment supported on a liquefiable foundation (medium sand) is investigated. The experimental series included: (i) a benchmark test without a liquefaction countermeasure, (ii) foundation densiflcation below the embankment toe, and (iii) use of a sheet-pile containment enclosure below the embankment. This series of experiments documents a wide range of practical liquefaction response mechanisms (including countermeasure implementation). In order to numerically simulate the above centrifuge tests, a new calibrated soil stress-strain constitutive model is incorporated into a two-phase (solid-fluid) fully coupled Finite Element formulation. Comparison of the computational and experimental results demonstrates: (i) importance of post-liquefaction dilative soil behavior in dictating the dynamic response and deformation characteristics of the embankment-foundation system, and (ii) capabilities and limitations of the numerical modeling procedure.     

SEISMIC AMPLIFICATION AT A SOFT SOIL SITE WITH LIQUEFIABLE LAYER

It is well known that local soil conditions play a key role in the amplification of earthquake waves. In particular, a liquefiable shallow soil layer may produce a significant influence on ground motion during strong earthquakes. In this paper, the response of a liquefiable site during the 1995 Kobe earthquake is studied using vertical array records, with particular attention on the effects of nonlinear soil behaviour and liquefaction on the ground motion. Variations of the characteristics of the recorded ground motions are analysed using the spectral ratio technique, and the nonlinearity occurring in the shallow liquefied layer during earthquake is identified. A fully coupled, inelastic finite element analysis of the response of the array site is performed. The calculated stress-strain histories of soils and excess pore water pressures at different depths are presented, and their relations to the characteristics of the ground motions are addressed.     

Numerical Modeling of a Church Nave Wall Subjected to Differential Settlements: Soil-Structure Interaction,Time-Dependence and Sensitivity Analysis

ABSTRACT

Historic masonry structures are particularly sensitive to differential soil settlements. These settlements may be caused by deformable soil, shallow or inadequate foundation, structural additions in the building and changes in the underground water table due to the large-scale land use change in urban areas.

This paper deals with the numerical modeling of a church nave wall subjected to differential settlement caused by a combination of the above factors. The building in question, the church of Saint Jacob in Leuven, has suffered extensive damage caused by centuries-long settlement. A numerical simulation campaign is carried out in order to reproduce and interpret the cracking damage observed in the building.

The numerical analyses are based on material and soil property determination, the monitoring of settlement in the church over an extended period of time and soil-structure interaction. A sensitivity study is carried out, focused on the effect of material parameters on the response in terms of settlement magnitude and crack width and extent. Soil consolidation over time is considered through an analytical approach. The numerical results are compared with the in-situ observed damage and with an analytical damage prediction model.     

Subsoil Interventions Effect on Structural Seismic Response. Part I: Validation of Numerical Simulations

Subsoil interventions to enhance the static soil resistance and reduce deformations may alter significantly the seismic response of the complex soil-foundation-structure system. The aim of this article is to have an insight in the physics of the problems encountered and to validate an adequate numerical modeling procedure to study these effects of the intervention in the global response of the system. Validation concerns wave propagation, site effects, and dynamic soil-structure-interaction issues including the intervention beneath the foundation. Theoretical models-expressions and experimental results from centrifuge tests have been used. The proposed numerical model is proven very efficient to describe the complex dynamic phenomenon and anticipate the seismic response after the employment of subsoil interventions.     

Effects of Near-Fault Pulses on Nonlinear Soil-Structure Systems

ABSTRACT

This paper is focused on effects of near-fault pulse characteristics on seismic performance of soil-structure systems considering foundation uplifting and soil yielding. To this end, an extensive parametric study is conducted. Mid-to-high-rise buildings of different aspect ratios (SR) resting on shallow mat foundations are investigated. Different vertical load-bearing safety factors (F_S) of foundation as well as different soil types (i.e. soft to very dense) are considered in this study. Finite element method is used for numerical modeling. The underlying soil is simply modeled with a set of nonlinear springs and dashpots beneath the foundation. Mathematical near-fault pulse models of fling step and forward directivity are used as input ground motions. The results show that reduction in structural drift demands due to nonlinear soil-structure interaction (SSI) is more considerable in the case of short-period pulses compared to long-period ones. In more precise words, significant reduction occurs when pulse-to-fixed-base period ratio falls within 0.7–1.5 in the case of directivity pulses and 0.5–1.4 in the case of fling pulses. It is also demonstrated that the beneficial effects of nonlinear SSI reduce when the number of stories increases.     

Development and Verification of a New Site Classification System and Site Coefficients for Regions of Shallow Bedrock in Korea

In this article, site response analyses for 124 sites collected in Korea were performed to evaluate earthquake ground motions in regions of shallow bedrock in Korea. Based on the results, a new two-parameter site classification system was developed for use in regions of shallow bedrock. The system incorporates depth to bedrock and mean V_S of soil above bedrock as parameters for site classification instead of V_S30. Soil sites were divided into seven site classes and the corresponding site coefficients were proposed for each site category. Some verification work demonstrates the superiority of the proposed system over the current seismic codes.     

The Ben Lawers Historic Landscape Project: Simultaneous Multi-element Analysis of Former Settlement and Arable Soils by X-ray Fluorescence Spectrometry

X-ray fluorescence spectrometry (XRFS) was used in the analysis of A horizon soil samples collected from a former farming settlement and its associated area of infield (i.e. arable) located in the Central Highlands of Scotland. To date, XRFS has not been extensively used in geoarchaeological research, but in our study the simultaneous multi-element capabilities of this instrumental technique allowed the total concentrations of 25 major, minor, and trace elements to be fully quantified with acceptable levels of accuracy and precision. Included within this group of chemical elements are a number (e.g. Ba, Ca, P, Pb, Sr and Zn) that have proved to be of value to archaeological interpretation in earlier investigations undertaken in Scotland. In our preliminary work documented here, significant differences were found between the A horizon soils of former settlement and infield areas for 18 chemical elements. Subjecting the XRFS data—and three other measured variables: soil organic carbon (SOC), pH and A horizon depth—to discriminant analysis indicates that soils of former settlement and arable farming can be effectively classified according to their pH, SOC content and Ca, Cu, Mg, Rb, and Zn concentrations. The inference is that areas of former infield and settlement elsewhere at this study location in the Central Highlands may be able to be identified according to their soil chemical composition and use of discriminant function, even though the surface remains of pre-eighteenth century settlement sites are not readily evident today because they were constructed of perishable materials.     

Calibration of an Advanced Constitutive Model for Babolsar Sand Accompanied by Liquefaction Analysis

In our research, we apply numerical modeling for prediction of liquefaction of sands during and after dynamic loading. In numerical modeling, to properly simulate the generation, redistribution, and dissipation of excess pore water pressure during and after dynamic loading, it is important to use a suitable constitutive model for soil. In this article, Dafalias and Manzari’s model [2004] (a critical state bounding surface plasticity model) was used to model the behavior of saturated sand due to relatively simple of formulations and a unique set of input parameters for a wide range of initial stress and void ratio. The attention in this article is on Babolsar sand. After calibration model parameters for Babolsar sand, the analysis of liquefaction using the modeling of a centrifuge test and predictions of model was carried out. The results indicate a reasonable performance of the model for prediction of behavior of types of sands. Also, Babolsar sand has more prone to dilatancy than Nevada and Toyoura sands.     

Rocking Soil-Structure Systems Subjected to Near-Fault Pulses

Seismic performance of rocking soil-structure systems subjected to near-fault pulses is investigated considering foundation uplifting and soil plasticity. An extensive parametric study is conducted including medium-to-high-rise buildings with different aspect ratios based on shallow raft foundation at stiff-to-rock sites. Mathematical directivity and fling pulses are used as input ground motion. The superstructure is assumed to have three different boundary conditions: (a) fixed-base, (b) linear soil-structure interaction (SSI), and (c) nonlinear SSI. Evidently, the prevailing pulse period T_p is a key parameter governing nonlinear SSI effects. The normalized acceleration response spectra reveal that despite beneficial effects of foundation uplifting and soil yielding in most cases, there are some minor regions in which the response accelerations are amplified. In addition, more slender buildings significantly benefit from uplifting and soil yielding when subjected to short- and medium-period directivity pulses compared to squat structures. However, response amplifications with respect to fixed-base structures are considerable in case of slender structures subjected to medium- or long-period directivity pulses. So that neglecting the SSI effects on seismic performance of rocking structures with shallow foundations, as mostly assumed in common practice, may give rise to inaccurate estimations of force demands against near-fault pulselike ground motions. Furthermore, the envelope of residual foundation tilting θ_r is limited to 0.015 rad, in case of directivity pulses.     

A quartic dispersion equation for internal gravity waves in the thermosphere

A new quartic dispersion equation in the square of the complex vertical wave number is derived by employing the ‘shallow atmosphere’ approximation and an ion drag approximation. These approximations allow the coefficients of the quartic equation to be given in terms of the corresponding cubic equation, which neglects the Coriolis force and the zonal ion drag component, but modified to take into account these neglected effects. Coupling between the extraordinary viscosity wave mode and the other three wave modes is highlighted and numerical solutions are compared for this quartic equation, an exact eighth order equation and the cubic equation. For the first time the validity of using the ‘shallow atmosphere’ approximation to describe internal gravity wave motions is demonstrated.     

Probabilistic Seismic Design of Pile Foundations in Non Liquefiable Soil by Response Spectrum Approach

The behavior of pile foundations in non liquefiable soil under seismic loading is considerably influenced by the variability in the soil and seismic design parameters. Hence, probabilistic models for the assessment of seismic pile design are necessary. Deformation of pile foundation in non liquefiable soil is dominated by inertial force from superstructure. The present study considers a pseudo-static approach based on code specified design response spectra. The response of the pile is determined by equivalent cantilever approach. The soil medium is modeled as a one-dimensional random field along the depth. The variability associated with undrained shear strength, design response spectrum ordinate, and superstructure mass is taken into consideration. Monte Carlo simulation technique is adopted to determine the probability of failure and reliability indices based on pile failure modes, namely exceedance of lateral displacement limit and moment capacity. A reliability-based design approach for the free head pile under seismic force is suggested that enables a rational choice of pile design parameters.     

Testing settlement models in the early Roman colonial landscapes of Venusia (291 B.C.), Cosa (273 B.C.) and Aesernia (263 B.C.)

This paper examines settlement density and settlement patterns in the Roman colonial territories of Venusia, Cosa and Aesernia, located in three different landscapes of central southern Italy (modern Basilicata, Tuscany and Molise). Using a series of GIS tools, we conducted a comparative analysis of the density and spatial distribution of sites dating to the Hellenistic period (ca. 350–50 b.c.). We used the legacy settlement data collected by previous large-scale, intensive, site-oriented field surveys to test the validity of two competing rural settlement models of early Roman colonization: the conventional model of neatly organized settlements regularly dispersed across the landscape and the recently proposed theory that colonists adopted a polynuclear settlement strategy. After calculating the extent to which the archaeological datasets conform to the regular or polynuclear model, we conclude that only a very small portion of the colonized areas actually meets traditional expectations regarding the organization of early colonial settlements. Our analyses show that the legacy survey data is more consistent with the polynuclear settlement theory, but the data also reveals some completely unexpected patterns, suggesting that early Roman colonial landscapes were more diverse than previously thought.     

Miliacin in palaeosols from an Early Iron Age in Ukraine reveal in situ cultivation of broomcorn millet

During the archaeobotanical investigation of Scythian–Sarmatian period (Early Iron Age), pits with crop processing waste, discovered in the floodplain of Donets River, eastern Ukraine, and charred remains of cereal grains, dominated by broomcorn millet, were recorded. The grains from the pits were radiocarbon dated to the fifth to first century BC. Those pits are distant from any known contemporaneous settlement. The apparent disconnection of these pits from any local settlement suggests that (1) millet was brought from other locations by mobile groups, or (2) millet was cultivated locally by populations whose settlements have left no discernible archaeological trace. The analysis of molecular biomarkers preserved in palaeosols that are stratigraphically connected to the pits revealed high levels of miliacin, a molecule that can be preserved in ancient soils and sediments, and that is consistent with broomcorn millet (Panicum miliaceum). High levels in miliacin in soils stratigraphically connected to the pits are interpreted as the result of a large biomass of P. miliaceum produced at time of soil formation. Our biogeochemical results applied to a palaeosol thus attest to the in situ cultivation of crops dominated by the broomcorn millet during the early Iron Age in the floodplain of Donets River. Biochemical examination of soils and palaeosols can thus provide useful information on past dynamics of land-use by ancient population, especially when settlements or macrobotanical remains are absent.     

Multiphase flow simulation through porous media with explicitly resolved fractures

Accurate simulation of multiphase flow in fractured porous media remains a challenge. An important problem is the representation of the discontinuous or near discontinuous behaviour of saturation in real geological formations. In the classical continuum approach, a refined mesh is required at the interface between fracture and porous media to capture the steep gradients in saturation and saturation‐dependent transport properties. This dramatically increases the computational load when large numbers of fractures are present in the numerical model. A discontinuous finite element method is reported here to model flow in fractured porous media. The governing multiphase porous media flow equations are solved in the adaptive mesh computational fluid dynamics code IC‐FERST on unstructured meshes. The method is based on a mixed control volume – discontinuous finite element formulation. This is combined with the P_N+1DG‐P_NDG element pair, which has discontinuous (order N+1) representation for velocity and discontinuous (order N) representation for pressure. A number of test cases are used to evaluate the method's ability to model fracture flow. The first is used to verify the performance of the element pair on structured and unstructured meshes of different resolution. Multiphase flow is then modelled in a range of idealised and simple fracture patterns. Solutions with sharp saturation fronts and computational economy in terms of mesh size are illustrated.     

MODIFICATION OF SITE RESPONSE IN PRESENCE OF LOCALISED SOFT LAYER

Experimental and numerical simulations are performed to evaluate the modification of ground response resulting from either the presence of soft layers or occurrence of partial liquefaction. Results from two densely instrumented dynamic centrifuge tests are presented to show the ambiguous role played by the presence of a soft layer. It was found that the lateral extent of the soft layer has significant influence on the overall response of the layered strata and any structure founded on it. The experimental observations are supported by simplified numerical analysis. The amplification or deamplification of the input motion is found to be a function of the ratio of the width of soft layer to the wave length. Based on the numerical analysis, a general function describing the site amplification is presented which may be used as a guide in seismic design of foundations in such layered strata.     

Energy use by Eem Neanderthals

An analysis of energy use by Neanderthals in Northern Europe during the mild Eem interglacial period is carried out with consideration of the metabolic energy production required for compensating energy losses during sleep, at daily settlement activities and during hunting expeditions, including transport of food from slain animals back to the settlement. Additional energy sources for heat, security and cooking are derived from fireplaces in the open or within shelters such as caves or huts. The analysis leads to insights not available from archaeological findings that are mostly limited to durable items such as those made of stone: Even during the benign Eem period, Neanderthals faced a considerable heat-loss problem. Wearing tailored clothes or some similar measure was necessary for survival. An animal skin across the shoulder would not have sufficed to survive even average cold winter temperatures and body cooling by convection caused by wind. Clothes and particularly footwear had to be sewn together tightly in order to prevent intrusion of water or snow. The analysis of hunting activity evolvement in real time further shows that during summer warmth, transport of meat back to the base settlement would not be possible without some technique to avoid that the meat rots. The only likely technique is meat drying at the killing site, which indicates further skills in Neanderthal societies that have not been identified by other routes of investigation.     

Classic Lowland Maya Political Organization: A Review

This paper reviews recent archaeological research concerning Classic Maya lowland political systems (ca. A.D. 250–1000). It focuses specifically on (1) subsistence practices revealed through the analysis of prehistoric climate, available resources, agricultural technologies, and diet; (2) population distribution, density, and size revealed through the analysis of settlement practices and architectural function; (3) social differentiation and interaction revealed through the analysis of burial practices, diet and health, architecture, and production, consumption, and exchange patterns; and (4) ancient Maya political economy (how it was funded) revealed through the analysis of community organization, ritual activities, the Classic Maya collapse, and warfare. It finally ends with a brief discussion of the future of Maya archaeology. A key factor that recurs throughout this review is the noticeable amount of variability that existed—varied resources, subsistence strategies, settlement practices, and social and political systems. An understanding of this variability is the key to appreciate fully the Classic Maya.     

Measurements of gradients in ionospheric parameters with a new nine-position experiment at Millstone Hill

A new nine-position experiment is now routinely carried out with the Millstone Hill incoherent scatter radars which allows estimation of spatial gradients in the measured ionospheric scalar parameters N_e, T_e, and T_i, and in the components of the ion velocity vector v_i. Use of this technique results in improved estimates of basic and derived parameters from incoherent scatter data at times of significant gradients. We detail the data analysis method and present the first results from this new experiment. The gradients in N_e and in the components of vi are used to compute the motion term in the ionospheric F region continuity equation ▿ · (Nv), which is then combined with ∂N/∂t to estimate the O⁺ recombination rate β at night. Meridional neutral winds U_mer are computed from the field-aligned ion velocity v_∥ and a calculation of the O⁺ diffusion velocity v_d, and it is found that horizontal gradients in the ion velocity field at times significantly affect the calculation of the neutral winds.     

Subsistence and mobility strategies in the Epipalaeolithic: a stable isotope analysis of human and faunal remains at 'Uyun al-Hammam,northern Jordan

This paper investigates evidence for subsistence and settlement activities in the Levantine Middle Epipalaeolithic through the application of stable isotope analysis to human and non-human remains from the site of 'Uyun al-Hammam, northern Jordan. In general, bone from the site suffers a high degree of diagenesis and collagen could not be extracted for analysis here. Carbon and oxygen isotopic values from human tooth enamel samples were variable, but within expected values for the Middle Epipalaeolithic, whereas animal carbon and oxygen isotopic values varied widely, most likely due to hydrological and climatological factors. Carbon and isotopic values for the human samples indicated a predominantly C₃ plant dietary input, while animal samples appeared to have varying amounts of C₄ inputs into their diet. This is the first isotopic analysis conducted on material from the Middle Epipalaeolithic of the southern Levant and, as such, even results constrained by temporal and climatological variants within the region contribute to the overall knowledge of settlement and subsistence strategies during this cultural period.